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Abstract 

A HIGH-LEUEL COMPUTER GRAPHICS IMPLEMENTATION OF 
THREE-DIMENSIONAL B-SPLINE SURFACE GENERATION 

by Dennis Michael Bryant 

Chairperson of the Supervisory Committee: 

Professor D. E. Calkins 
Department of Mechanical Engineering 

The difficulty of accurately defining, representing 
and visualizing three-dimensional sculptured surfaces has 
concerned designers and engineers throughout history, but 
never more so than today. Technology has given man the 
ability, and the necessity, to build automobile bodies, 
ship hulls and aircraft fuselages with complex curvature 
and extremely smooth lines. The ability to efficiently 
develop and define the form of such surfaces has not kept 
pace with the requirement to manufacture them. 

This study explores the application of a surface B- 
spline for surface definition coupled with high-level 
computer graphics to enhance visualization. A computer 
program, SPLINE, is introduced which facilitates 

interactive construction of the B-spline control graph and 
real-time manipulation of the model on an Evans and 
Sutherland PS 300 distributed graphics system. A uniform 
nonrational surface B-spline is utilized for interactive 
modeling, with the capability preprogrammed of 

implementing nonuniform rational surface B-splines upon 



later addition of an inversion algorithm. 

While the program is not fully functional, due to 
apparent shortcomings of the host computer operating 
system, the viability and desirability of such a system is 
fully demonstrated. Heuristic techniques are discussed 
for approaching model development, based on properties of 
the B-spline and characteristics of SPLINE. Step-by-step 
illustrated case studies further develop the modeling 
techniques. Examples include the generation of B-spline 
representations of a sailing yacht hull and of a Porsche 
automobile body, as well as the freeform development of a 
sculptured surface ab initio. 

Structured and heavily commented program listings are 
included, with flow diagrams, function networks, and 
operating instructions. Program expansions and hardware 
improvements are recommended and existing provisions 
within the program for these improvements are identified. 
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CHAPTER 1 



SCULPTURED SURFACE REPRESENTATION 
1.1. Sculptured Surfaces 

The computer program described herein uas developed 
to demonstrate the feasibility of designing and defining 
sculptured surfaces through the use of a B-spline 
algorithm coupled with highly interactive computer 
graphics. The accurate definition, representation, and 
visualization of three-dimensional surfaces is a problem 
which has concerned designers and engineers throughout 
history. Even relatively simple three-dimensional shapes 
can often be difficult to depict on a two-dimensional 
drawing such that the viewer can readily generate a mental 
image of the object or the craftsman can construct the 
object true to the designer’s intentions. 

Particularly troublesome, however, are sculptured 
surfaces, those having curvature in more than one 
direction. Sculptured surfaces form a vital part of 
aircraft, ship, and automobile designs, being necessary in 
aircraft and ship design for aerodynamic and hydrodynamic 
streamlining and used extensively in automobile design 
primarily for styling reasons. Nor is the use of these 
troublesome surfaces for the sake of appealing to the eye 
restricted to such major market items as automobiles; as 
Rogers and Satterfield C1J point out, today even the 
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shampoo bottle is likely to sport a complex sculptured 
surface . 

Given the difficulty of representing these surfaces 
in a tuo-dimensional medium, it is not surprising that the 
historical solution has often been simply to avoid 
producing such drawings at all. Letcher CHI presents a 
brief but interesting history of man’s attempts to skirt 
the problem in the context of ship hull design. The form 
of the earliest hulls, and of many small craft in low- 
technology cultures today, was undoubtedly preserved by 
carving a solid model from wood. Later, the use of the 
master-section-and-batten method was developed and has 
continued to some extent until today. A few transverse 
sections are erected with the desired form, and the 
remainder of the hull is defined by the installation of 
the hull planking as it is sprung about the master 
sections . 

Even as late as the early 1300’s, the lines drawing 
was not the principal tool in developing and visualizing 
the lines of a ship’s hull. While the lines were 
ultimately preserved with such a drawing, they were 
usually taken from a carved model. Often the model was 
even carved in a material with separable layers so that 
the lines could be directly traced onto the paper. 

1.2. The Ship’s Lines Drawing 



The problem of defining and representing a sculptured 



3 



surface is perhaps best illustrated by a description of 
the lines drawing, that instrument which is used to 
develop and define the surface of ship and boat hulls. As 
is standard with the representation of virtually all 
three-dimensional objects, three views of the hull are 
presented on the lines drawing: a top view Cknown as the 
half-breadth plan, because only one half of the 
symmetrical hull is shown } ; a front view Cknown as the 
body plan, and which actually is a stern view as well}, 
and a side or profile view Clabeled the sheer plan}. 

Intersections of the vessel’s hull with evenly spaced 
cutting planes are projected onto these three reference 
planes. On the body plan, the form of the hull at 
transverse sections is shown. Depending on the size and 
complexity of the vessel’s hull, the sections are shown at 
ten, twenty or forty foot intervals, with a few 
intermediate stations shown in areas of greatest 
curvature. For the sheer plan, vertical cutting planes 
are taken at one to four foot intervals, yielding a set of 
buttocks lines. Waterlines are similarly generated on the 
half-breadth plan by intersecting the hull with vertical 
longitudinal cutting planes. Nonorthogonal cutting planes 
further serve to define the shape through the production 
of diagonals or cant lines. A sample lines drawing is 
shown in figure 1. 

The task of the naval architect in developing a 
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completed set of lines is not an easy one. Having 
visualized the desired Form, the architect must accurately 
translate it into a set of lines on paper which satisfies 
quantitative, qualitative, and aesthetic considerations. 
Working with irregularly shaped curves, the designer must 
nevertheless maintain specific areas within the curves to 
yield proper ship stability and handling characteristics. 
In order to obtain fair (smooth} lines, the lines are 
drawn with a drafting spline, an elastic lath constrained 
by weights called ducks. The lines must not only be 
smooth and contain the appropriate areas, but the 
projection of all points on the lines onto all three 
reference planes must be consistent. These requirements 
often appear to be nearly mutually exclusive. 

In short, because of the difficulties of Cl} visual- 
izing the sculptured surface, C2} obtaining smooth lines, 
C3} maintaining agreement of all views, and C4} achieving 
and measuring the required form characteristics such as 
section areas, the design of a ship’s lines is a tedious 
and difficult iterative task for even the most experienced 
architect. Letcher C23 cites Kinney as saying, ”’It takes 
between five and ten working days for an experienced man 
to produce a good set of lines, ’ regardless of the size of 
the vessel . ” He also points out that anything more than a 
slight change requires starting over almost from the 
beginning . 
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When the naval architect has finished the lines 
drawing, however, the task of hull definition is not yet 
complete. The lines are next translated into a tabulated 
form called a table of offsets. This table lists the 
distances, at each station, of waterlines from the 
centerline and of buttocks lines above the baseline. An 
example is shown in table 1 . 

Greater accuracy is needed for production purposes 
than is obtainable from the relatively small scale lines 
drawings. For this reason the lines are redrawn full size 
on a large wooden floor Cthe mold loft} and faired again. 
In very recent practice, the full scale lofting may be 
replaced by drawing the lines in one-sixteenth scale on a 
marble slab with a very hard pencil or by entering the 
offsets into a computer for generation of a large scale 
drawing on a plotter. In any case, the large scale 
fairing is used to identify necessary revisions to the 
lines and to produce a more accurate table of offsets. 

1.3. Mathematical lu Defined Surfaces 

It is not hard to understand, then, the benefits 
which would be derived if the hull form, or shampoo 
bottle, could be defined by a precise mathematical 
relationship. The need for full scale lofting would be 
eliminated, since accuracy could be obtained initially to 
whatever degree was desired. Precision would no longer 
depend upon the ability to measure small increments on a 



TABLE 1 



Table of Offsets for a Mariner Class Uessel 
Source: Owen C31 



(Dimensions given in feet, inches and eighths) 



Half Breadths at Waterlines and Decks 



Station 


O’-O” 

UJL 


4 ’ -0 ” 
UJL 


a ’ -o ” : 

UJL ! 


; 16 ’ -o ” 

UJL 


CFUD END) 








1 


0 FP 


— 


3-0-5 


2-B-2 


1-0-1 


1/2 


— 


4-0-7 


5-0-7 


4-3-5 


1 


— 


6-10-2 


7-11-6 


0-6-1 


2 


2-4-0 


13-7-6 


16-9-3 


19-9-0 


3 


9-2-0 


24-1-0 


27-7-1 


30-10-3 


4 


21-4-0 


33-3-0 


35-0-1 


37-2-2 


5 


2B-5-6 


36-2-6 


37-10-3 


30-0-0 


S 


24-10-4 


34-B-l 


37-0-7 


37-11-7 


7 


13-6-4 


27-1-4 


31-2-6 


35-0-1 


B 


4-5-2 


15-0-3 


19-3-3 


26-2-2 


9 


— 


6-0-7 


7-9-4 


11-2-4 


9 1/2 


— 


2-B-7 


3-4-1 


4-3-0 


AP 10 
CAFT End) 


— “ 


““ “ 


“ ““ 







Heiahts Above Molded Baseline 


Station 


4-ft 

buttock 


0-ft 

buttock 


16-ft ! 

buttock ! 


I 24-ft 
buttock 


CFUJD End) 
0 FP 


47-5-6 


57-0-7 


1 


i 


1/2 


2-3-0 


40-0-3 


55-4-4 


— 


1 


0-11-7 


0-0-7 


43-4-1 


56-0-1 


2 


0-0-4 


0-0-7 


6-9-2 


31-5-7 


3 


— 


0-0-0 


0-6-7 


3-11-1 


4 


— 


— 


0-0-0 


0-1-1 


5 


— 


— 


— 


0-0-0 


6 


— 


— 


— 


0-0-0 


7 


— 


— 


0-0-6 


2-1-0 


0 


o-o-o 


0-3-5 


4-9-6 


13-3-3 


9 


1-2-1 


0-6-6 


21-11-3 


30-0-4 


9 1/2 


14-3-4 


23-4-1 


30-3-6 


30-4-5 


AP 10 


20-5-0 


31-6-5 


30-11-3 


— 


CAft End) 


37-1-4 


41-7-0 


44-1-6 


— 



a 



drawing, but could be achieved by simply outputting a few 
more digits from the computer’s memory. 

Preparation of a table of offsets would not be a 
manual task, but simply a matter of solving the 
mathematical relationship at any desired point on the 
surface. Thus the representation would not only be more 
accurate, but could be more complete as well. A lines 
drawing only accurately defines the surface at paints 
which fall exactly an the lines drawn. Determination of 
intermediate paints requires drawing and fairing 
additional lines or performing some manner of 
interpolation. Interpolation will yield only a rough 
approximation, given the complex curvature of the farm. A 
mathematical representation can be equally accurate at any 
paint which is defined by the mathematical relationship. 

A mathematical representation could also provide 
direct input for the calculation of stability data, as 
well as for fabrication. The need for scale model testing 
would be greatly reduced with the substitution of computer 
simulations, which would no longer need to rely upon mere 
approximations of the form. 

Perhaps most exciting are the prospects for using a 
mathematical representation of a sculptured surface as 
input to a numerically controlled machining CNCM3 process. 
Rogers, Rodriguez, and Satterfield have already reported 
substantial progress in this application C4,5,SD. A 
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system of computer programs was developed on a 
microprocessor-based graphics system interfaced with a CNC 
controller to drive a milling machine for the production 
of wooden towing tank models. The numerical control 
milling data base is actually generated directly by a 
computer program from a mathematical definition of the 
hull. This mathematical definition is created by the 
program in response to manipulation of a graphical image 
on the computer screen by the user. Using this system, 
towing tank models havB been produced in only two days. 

Considering the beleaguered state of the shipbuilding 
and automotive industries in the United States today, it 
is impossible not to become excited with the possibilities 
of larger scale implementations of this concept. 

1.4. High-Level Computer Graphics 

The implementation of mathematical definition of 
sculptured surfaces does not require enormous computing 
power. The FAIRLINE series of computer programs developed 
by Letcher C2,7,83 for the design of sailing yacht hulls 
was designed for use on a Texas Instruments pocket 
calculator. On the other hand, defining the surface is 
only part of the problem. Closely related is the problem 
of representing the surface during the surface development 
process and subsequent to surface definition. It may be 
difficult to visualize a three-dimensional surface while 
looking at a two-dimensional drawing, but it is even more 
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difficult to visualize a complex sculptured surface by 
reviewing a computer-generated table of coordinates. 

It is when coupled with high-level computer graphics 
that the mathematical representation begins to be as 
useful as it is powerful. High-level computer graphics 
have been defined by Calkins C93 as being characterized 
by: 

C1D the ability to display three-dimensional 
objects , 

C2) real-time image manipulation of the three- 
dimensional objects displayed, and 

C3) color-shaded object surfaces. 

Calkins summarizes high-level computer graphics as the 
’’capability of real time manipulation of realistically 
rendered 3-dimensional objects." CIt will be shown later 
that while the computer system utilized for the research 
described in this paper does not meet the third criterion, 
it certainly fulfills the intention of this summation.) 
With high-level computer graphics, the designer can rotate 
the object about to obtain another perspective or zoom in 
on a troublesome detail. The results of changes made in 
one view can be immediately reviewed from other angles to 
determine their suitability. 

It is this merging of visualization of sculptured 
surfaces through high-level computer graphics and 
definition by mathematical representation which has been 
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explored in this research. 



CHAPTER 2 



B-SPLINE SURFACES 

The form of mathematical surface representation 
selected for implementation in the computer program under 
consideration is the surface B-spline. In particular, a 
nonuniform rational surface B-spline algorithm has been 
utilized, though in its present form the program only 
allows the user direct control of a subset of this, the 
uniform rational surface B-spline. 

The B-spline concept was developed by Schoenberg [101 
in 1346 and introduced to computer-aided geometric design 
by Riesenfeld Clll in 1973. Several properties of the B- 
spline curve or surface generation algorithms make them 
very suitable for sculptured surface modeling using 
interactive computer graphics. Briefly, the B-spline 
offers inherent curve smoothing, considerable data 
compression [complex surfaces can be fully defined with a 
small amount of data), local control Ca portion of a curve 
or surface can be modified without affecting the entire 
surface!, the flexibility of including sharp 
discontinuities in a curve or surface, and a convenient 
means of interactively manipulating the curve. 

This last property will be shown to be a particularly 
useful characteristic. Barsky and Greenberg C1B3 discuss 
the difficulty of performing shape design using 
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mathematical methods which require the designer to 
directly specify values for algorithm variables. It can 
require an intimate knowledge of complex mathematics for a 
designer to be able to readily perceive the necessary 
adjustment to a mathematical algorithm which will effect 
the desired change in a portion of a surface. In 
contrast, modeling with B-splines requires virtually no 
knowledge of mathematics by the program user. Fish C133 
observed, ’’The nice thing about spline curves is that they 
allow you to treat a collection of pieces of various 
polynomials as a single entity, maintaining continuity of 
tangent directions, curvatures, and so on without working 
at it.” Undoubtedly the same considerations led Rogers 
C43, working with a number of mathematical curve 
representations, to conclude that the B-spline was ’’most 
generally useful.” 

Uhile a detailed understanding of the mathematics 
behind B-splines is not necessary in order to use them, 
some understanding will nevertheless aid the user in 
developing the ability to more quickly model a desired 
surface. A brief introduction follows. flare extensive 
treatments of the subject can be found in Rogers and Adams 
C14J, who provide useful examples, and in Bartels, Beatty, 
and Barsky C15J, who provide a particularly methodical and 
understandable approach. 
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5.1. Parametric Curves 

The Cartesian coordinate system is essential to most 
surface design. It is in this coordinate system that most 
of the formulations are presented with which the naval 
architect or aircraft engineer will work to determine the 
performance characteristics of the craft. It is also 
Cartesian coordinates in which the graphics computer will 
require data to be presented for display and manipulation 
on the screen. At the same time, the Cartesian coordinate 
system has significant liabilities when working with 
sculptured surfaces. 

Host significantly, many sculptured surfaces or 
curves cannot be defined by an explicit function Ca 
function having only one value of / for each value of x) . 
The transverse section of a ship with tumble-home is a 
ready example. Also, it will be shown that calculating 
one Cartesian coordinate for a curve at even increments of 
another Cartesian coordinate can result in very irregular 
definition resolution. For these reasons, a parametric 
representation is preferable. The B-spline can be 
parametric. In a parametric representation, the Cartesian 
coordinate for a point on a curve is not defined in terms 
of the other Cartesian coordinate Ci.e. y - fCx)), but 
coordinates are instead defined as functions of a 
parametric variable: thus, 



x 



f Cu) 
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and 



y " gfu) . 

where a is the parametric variable. 

The parametric representation addresses both of the 
objections to the purely Cartesian coordinates. Implicit 
functions in x and y can be handled without ambiguity and 
the resolution of the curve representation is generally 
much more uniform with equal parameter spacing than with 
equal Cartesian coordinate spacing. 

Rogers and Adams 111411 provide an illustration similar 
to figure B, which compares representations of a circular 
arc using a nonparametric representation with data points 
at equal x intervals, using uniform parametric spacing, 
and using equal arc lengths. The nonparametric 
representation uses the relationship, 

y - 4Cl-x=J . Cl) 

The uniform parametric spacing uses a parameter, a, which 
varies from zero to one. The points on the curve are 
defined by 

x - Cl-u=)/Cl+u=) 

and 



y - Bu/Cl+u=) . CB) 
The equal arc length representation defines the points on 
the curve as 

x “ COS0 



and 
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y - sin0, C3) 
In each case, the arc is represented by a straight-line 
connection between six points. 




Fig. 2. Comparison of parametric and nanparametric 
representations . 



While the uniform parametric spacing does not produce 
a resolution quite as uniform as that of equal arc 
lengths, it is clearly superior to the equal x coordinate 
spacing. Sharpe and Thorne C161 discuss an approach for 
using the equal arc length spacing with a parametric 
representation, but more straightforward equal parameter 
spacing was felt to be sufficient for the purposes of this 
program . 

2.2. Two-Dimensional B-Splines 

The B-spline is most easily understood initially by 
looking at a two-dimensional uniform nonrational case, the 
simplest form of the B-spline. The B-spline curve is a 
piecewisB polynomial curve. That is, the curve is broken 
into segments which are defined by separate polynomials. 
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This curve is everywhere continuous, both in position and 
in up to the CK-H) derivative, K being the order of the 
spl ine . 

The general B-spline curve is defined by: 
n 

QCu) - E BiNi.kCu) C4) 

1-1 

Ni.xCu) - 1 if Xi<_u<Xi +l 
- 0 otherwise 

(u-Xi)Ni,k-i(u!) C x i +k-u )Ni*i . k— i (u } 

Ni, k (u) - + 

Xi*k-Xi + i 



where 



QCu} - CxCu} , yCu}D , the curve, 
n - number of control vertices, 

B i - control vertices’ coordinates, 

K ” order of the B-spline curve C implies 
degree CK-1 } } , 

Ni.kCu) - B-spline basis functions, 

and 



Xi - elements of the knot vector. 

While it can be seen that the B-spline algorithm can 
be somewhat confusing mathematically, the use of it is 
actually quite intuitive. The curve being formed is 
shaped by controlling a relatively small number of points 
known as control points or vertices . The curve will not 
interpolate Cpass through} the control vertices, but will 
mimic the overall shape of the array of vertices. For a 
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better visualization of this, the vertices are frequently 
Joined by straight line segments, forming a control 
polygon. An example control polygon with its associated 
B-spline curve is shown in figure 3. The vertices are 
represented by the squares, the control polygon by the 
connecting straight line segments, and the B-spline curve 
by the curved line. 




Fig. 3. Third order B-spline curve with control 
polygon . 



An apt analogy to the effect of the control vertices 
upon the curve is that of a set of magnets to a ferrous 
wire. While the curve does not pass directly through the 
control vertices, it is molded by their placement, being 
drawn toward them as if toward magnets. (loving a control 
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vertex pulls the curve in the vertex’s local region in the 
direction of the movement. Placing vertices close 
together creates a stronger effect on the curve in the 
region of those vertices. 

Every point on the curve is a weighted average of the 
position of a number of control vertices. Each vertex 
has a unique basis function, which defines the weight of 
the vertex with respect to the parameter a. These basis 
functions, polynomials in the parametric variable, are 
independent of the shape of the controlling polygon, being 
completely defined by the order C k~) and a knot vector. 

The knot vector divides the curve parameter Ca} into 
a series of discrete values. This subdivides the B-spline 
curve into a series of segments defined by separate 
polynomials. The Joined segments form a piecewise 
polynomial curve. If the divisions in the curve are 
equal Cthe knot vector specifies uniform step width 
between knots}, the B-spline is said to be uniform. 

The parameter a varies between an initial value and 
some final maximum as the curve is traversed. Ualues of a 
that correspond to Joints between the spans are called 
knots. In the case of a uniform B-spline, the knots will 
typically correspond to integer values of a. The knot 
vector gives the sequence of these knots and is 
constrained to be non-decreasing. Thus the knot vector 
defines the range of the parameter a which corresponds to 
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each segment of the curve . 

When the knot vector is specified, a multiplicity K 
at each end of the knot vector ensures that the ends of 
the curve interpolate the first and last vertices in the 
control polygon. This is a very convenient constraint 
uihen performing modeling with a B-spline. The knot vector 
for a uniform third order B-spline with four vertices, 
which varies from zero to Cn-K+11, is CO, 0, 0, 1 , 2, 2, 21 . 

The basis functions, then, indicate the degree and 
extent of the influence of each control vertex. This 
provides the local control which was mentioned above. The 
maximum extent of influence of a single control vertex in 




Fig. . Localized effect of moving a single control 
vertex . 




21 



parameter space is plus or minus half the order of the 
basis function. Conversely, the shape of every span is 

affected by K successive vertices. Figure 4 illustrates 
the localized effect of moving a single control vertex. 

A comparison of figures 5 and 6 will demonstrate the 
effect of curve order upon the basis functions. Figure 5 
shows the values for the basis functions for a third order 
curve with five vertices. Figure 6 shows the basis 
function values for a fourth order curve, also with five 
vertices. In each case, the parameter a has been 
normalized to unity. Notice that the basis functions for 
the fourth order curve are generally lower in value at any 
given value of a, but are nonzero over a greater range. 
Thus a higher order basis function yields a curve which is 
less markedly influenced locally by the movement of a 
single control vertex, but which is smoother and more taut 
than a lower order curve. This effect can be readily 
visualized by examining figures 7 and 8, which show third 
and fourth order B-spline curves created by the same 
control polygons. 

A few characteristics of the B-spline curve will be 
convenient to keep in mind when manipulating a curve: 

Cl } The B-spline curve will have at most as 
many extreme and inflection points as thB 
control polygon. 

C2D A straight line will be produced in the B- 
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Fig. 6. Fourth order basis functions for five vertices. 




Fig. 7. 



Third order B-spline curve. 
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spline curve by K ccllinear central 

vertices . 

C3) K-l callinear ccntrol vertices will yield a 
curve tangent to the control polygon. 

C4) A knuckle, or sharp discontinuity, uiill 
occur at the point where K-l control 

vertices are colocated. An example is 
shown in figure 9. 




Fig. 9. Knuckle produced by colocated vertices. 



C5) Each span of a B-spline curve of order K 
lies within the convex hull of its K 
associated vertices. An illustration of 
this concept is shown in figure 10, where 
the B-spline curve would lie in the shaded 
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area . 

As previously mentioned, an advantage of the B-spline 
is the resultant data compression. Rogers and Satterfield 
C1D present an example wherein they successfully 
represented the forebody of a U.S. Navy ammunition ship 
using only a five by eight polygonal net with two- 
dimensional uniform B-splines. 

£.3. Surface B-Splines 

Nearly all systems developed for computer-aided 
sculptured surface design have defined the surface in 
terms of a net of lines. Typically, a series of B-spline 
curves in parallel planes is intersected with another 
series of curves in planes perpendicular to the first. 
The result is not a true surface definition, but 
definition of a net of lines lying on the surface. To 
define any point not lying on a net line requires an 
interpolative approximation. 

The program which is the subject of this study uses a 
true surface algorithm rather than such a network of two- 
dimensional representations. The surface created is a 
mosaic of surface patches, Just as the two-dimensional 
curve was a piecewise polynomial curve, and the surface, 
like the curve, is everywhere continuous with continuity 
of the derivatives to a degree dependent upon the order of 
the basis functions. Control is achieved through 
manipulation of a net of control vertices, with the 
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pi ^^ ggntrol VgrtBH 





Fig. 10. The convex hull of a third order Cal, 
fourth order Cbl and fifth order Cel polygon. 
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control graph being analogous to the control polygon of 
the two-dimensional case. The control graph, as well as 
the surface, is three-dimensional. 

The surface Cor tensor-product) B-spline is defined 
very similarly to the two-dimensional spline as: 

n m 



QCu, w) 
Nx, xCu) 

Ha. i Cu) 

Nx.kCu) 

hc.^Cw) 



E £ B i . jN i „ « C u) tla , i_ C w ) 
1-1 J“1 



C5) 



- 1 if Xi<u<Xi*i 
“ 0 otherwise 

“ 1 if yj<_u<yj*i 

- 0 otherwise 

^ * ) N x , k — i C u ) C x i +k _ u )N x + i , k— i (u) 



Xi+K-i-Xx X i X x ♦ i 

Cw-y a Cyj+L-w)Nj + i, L -i(w) 



yj+L- i -ya 



yj*L _ yo* i 



where 



Q(u,w) 

u , w 
n 



m 

Bx . j 
K 



CxCu , w) , y Cu , w) , zCu, w) 1 , surface data 
points, 

parametric variables, 

number of control vertices in the a 
direction , 

number of control vertices in the w 
direction , 

control graph points, 

order of the surface in the a 



direction , 
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and 



L 

Ni , k(u)“ 
Ha.i_Cuj3- 

Xx 



order of the surface in the m 
direction , 

basis functions in the a direction, 
basis functions in the w direction, 
elements of the knot vector in the u 
direction, 



yj " elements of the knot vector in the m 
direction . 

The computer program developed for this study 
implements a very direct approach to surface calculation, 
computing all basis functions initially, then performing 
the summation of equation C5) . A more efficient 
algorithm, based on the work of Cox and de Boor, is 
presented in appendix C-27 of Rogers and Adams C143 . 

The surface B-spline is directly analogous to the 
two-dimensional B-spline in virtually every way, having 
very similar properties extended into three dimensions. 
For example, just as the effect of moving a single control 
vertex was seen in two dimensions in figure 4, the effect 
in three dimensions is shown in figure 11. CThe effect 
extends toward the back of the picture the same as toward 
the front, although some of that effect is lost to view in 
this illustration due to depth cueing.) 

One difference between the two and three-dimensional 
B-splines is the effect of a repeated knot vector. A 
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Fig. 11. Localized effect of moving a single control 
vertex in a surface B-spline. 



knuckle can be introduced into a two-dimensional curve by 
use of a repeating knot value in the knot vector. The 
same is not true with a three-dimensional surface, unless 
the knuckle is desired across the entire surface, since 
the same knot vector is used for all sets of vertices in 
the a Cor wl direction. Repeating a knot for one or two 
positions in the control graph will result in repeated 
knot values for the corresponding positions in all sets of 
vertices in the net. For this reason, knuckles will be 
introduced into the surface in this program only by 
colocating vertices. Multiple vertices will not indicate 
repeated knot values. This concept is dealt with in more 
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depth by Rogers and Satterfield Cl, 173. 

5.4. Nonuniform Rational B-Splines 

The B-spline algorithm implemented in this program is 
a still more general case of the surface B-spline, the 
nonuniform rational B-spline CNURBD . Both non-uniformity 
and rationality extend the power of the algorithm, as 
explained by Tiller C1S3. 

A limitation of nonrational B-splines is the 
inability to accurately represent circles, conics, and 
quadric primitives. The rational B-spline can do so. 

Rationality refers to the addition of a fourth coordinate. 
Given a point in three-dimensional Euclidean space, 
P~Cx,y,zD, a corresponding point is notu defined in four- 
dimensional space, P M “CHx , Hy , Hz , HD . The H is the 
homogeneous coordinate and is constrained to be non-zero. 
By adding a homogeneous coordinate to the definition of 
each control vertex, the effect of that vertex is 
effectively scaled. Thus the B-spline algorithm now 

becomes : 

n m 

Q m Cu,uj) - £ £ B H i, J Ni >K Cu)n J , L (w). C6) 

1-1 J-l 

Application of the rational B-spline to the development of 
surfaces is further dealt with by Tiller. 

As alluded to above, the knot vector spacing can also 
be used in some cases to control the surface shape. Thus 
the algorithm implemented in this program mill allow 
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nonuniform knot vectors. 

Because of their tremendous power and flexibility, 
nonuniform rational B-splines have been recognized now as 
an IGES standard for curve and surface definition C19J. 
The Alpha_l project at the University of Utah Computer 
Science Department, working with Riesenfeld, has done much 
work in exploring the potential of NURB’s C13J. 

On the other hand, manipulation of the homogeneous 
coordinates and of the knot vector does not produce 
intuitively predicted results. Since this negates one of 
the greatest advantages of the B-spline for interactive 
modeling, the user of the computer program developed in 



this work 


has not been given direct control over these 


variables . 


Uniform knot spacing, and a homogeneous 


coordinate 


of one, are default values for interactive 


modeling . 


By including these variables in the algorithm, 


though not 


under user control, the potential of utilizing 



nonuniform rational B-splines in an inversion algorithm 
Cdiscussed in chapter 6) is retained. 



CHAPTER 3 



SPLINE: A B-SPLINE SURFACE MODELING PROGRAM 

The end product of this study is a computer program 
for the Evans and Sutherland PS 300 computer graphics 
workstation which Facilitates the design and 
representation of sculptured surfaces using a surface B- 
spline. This program is known as SPLINE. 

3.1. Hardware 

The PS 300 is a distributed graphics system, 
employing its own microprocessor for the processing of 
display-related data and relying upon a host computer for 
other operations. The host used in the development of 
this program is the PDP-11/44 with an RSX-11M+ operating 
system. The computer is the University of Washington’s 
Mechanical Engineering Department computer and is set up 
on a time sharing system. Memory for execution of 
individual programs is limited to 64 kilobytes, with 
extension to IBB kilobytes possible by specifying an ID 
flag during compilation. This flag overrides the 
reservation of memory normally restricted to system 
housekeeping tasks. 

The PS 300 station which was used is the single-user 
version of the Evans and Sutherland graphics workstation. 
Properly designated the PS 330, this unit is the most 
basic system of the PS 300 family and is commonly referred 
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to by the family designation. The unit upon which this 
program was developed was running P5 Version S firmware, a 
preliminary version of the operating system and 
proprietary programming language. The PS 300 
workstation is pictured in figure IE. As can be seen 
there, the workstation is comprised of the control unit, 
display screen, keyboard, data tablet, and control dials. 



Fig. IE. P5 300 workstation. 

The control unit consists of a Graphics Control 
Processor, mass memory, and the Display Processor. The 
Graphics Control Processor is a Motorola G8000 
microprocessor with E4-bit address space and S56 kilobytes 
of local memory for the graphics firmware. A full 
megabyte is available in mass memory, with 



16-bit 
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communication with the Graphics Control Processor and 33- 
bit communication with the Display Processor, for very 
fast display access. 

The PS 300 display screen is a monochrome vector 
refresh display with B1SE by BISS resolution. White lines 
and characters are displayed on a black screen with sixty- 
four levels of intensity. While solids with color shading 
cannot be displayed, the fine scaling of intensity allows 
display of wire frame models which are unusually easy to 
interpret. Through a feature called depth cueing, the 
intensity of lines is gradually decreased as their 
distance behind the plane of the screen is increased. 
This fading with depth helps to give the perception of 
three dimensions while viewing a model on the two- 
dimensional screen. 

An additional benefit of depth cueing is that wire 
frame lines in the background can be clipped from view, 
eliminating confusion with foreground lines. This is 
especially important since the PS 300 does not provide 
hidden line removal. A comparison of figures 13 and 14 
will illustrate this advantage. 

The keyboard doubles as an input device for the 
PS 300 and a terminal for the host computer. Features 
include a full keyboard and twelve programable function 
keys with LED labels. A separate pad with eight 
programable control dials also features LED labels. 
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Fig. 13. Dire frame model without depth cueing. 




Fig. 14. Wire frame model with depth cueing. 
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The coordinate position of a stylus on the data 
tablet is converted to a digital equivalent for the 
positioning of a cursor on the screen. This can be 
programmed for menu selection or picking of points or 
lines on a displayed model. 

The PS 300 performs graphical manipulations and the 
associated mathematical transformations independently from 
the host. The transformations and matrix concatenations 
are accomplished without the user having to program the 
matrix arithmetic. 

While graphics manipulations can be performed 
without host interaction, communication between the host 
computer and the PS 300 is essential for virtually all 
applications. The PS 300 has no local storage capability, 
so its program must be stored an the host system and 
passed to the PS 300 ’ s memory as ASCII data. The 
programming language available on the PS 300, which is 
unique to Evans and Sutherland systems, is highly 
efficient for graphics display and manipulation, but other 
mathematics and data processing which would be 
straightforward in another high level computer language 
can be very difficult. The programmer must make a 
division of the graphics application into graphics 
processing and analysis functions in accordance with each 
computer’s strengths. An intelligent division can 
minimize communication between the computers, but 
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substantial communication will remain necessary. 

Communication between the PS 300 and the host is 
accomplished through a 9600 baud interface. Data transfer 
protocol is established through the PS 300 Host-Resident 
1/0 Subroutines CPSID’s). These are host-resident FORTRAN 
subroutines which provide two-way communication between 
the computers . 

3.3. Programming the PS 300 

The PS 300 has its own high level ASCII command 
language with English-like commands, detailed in the Evans 
and Sutherland manual , PS 300 Computer Graphics Sustem 
C203 . Structured objects are created and accessed by name 
with commands such as UECTOR_LIST Cto create an object as 
a list of vectors), TRANSLATE, ROTATE, and SCALE Cto 
perform the named operations) . 

Models are created as hierarchical groupings of 
graphical data, transformations applied to the data, and 
attributes such as intensity and level of detail. 
Primitive objects are defined and used to build more 
complex objects. For example, an automobile wheel might 
be defined by means of a vector list, then instances of 
that wheel translated and rotated to create four wheels 
and position them relative to a separately defined body, 
forming a representation of a car. A sample display tree 
is shown in figure 15. The three basic forms of display 
tree nodes can be seen there: data nodes Csuch as vector 
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Fig. 15. Display tree for an automobile. 
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lists defining objects} , operation nodes Csuch as 
translations and rotations), and instance nodes. Instance 
nodes group other elements together under a single name. 
Applying a transformation to an element results in the 
transformation of all subordinate elements as uell. 

flanipulation and control of models is accomplished 
through the creation of function networks. Ualues from 
interactive devices such as function keys or control dials 
are processed through user-designed function networks and 
applied to interaction points in the model’s structure. 
The function networks are combinations of connected black 
box functions. These function blacks accept data, 
manipulate it according to predefined rules for the given 
function, and output new values to another function input 
or to interaction points in the display tree. 

For example, ADDER might be defined as an instance of 
the function, F : ADDC , which sums any number sent to its 
first input and a constant value held on its second input. 
If the function keys are connected to input one of ADDER 
and the integer ”2” is sent to the second input, whenever 
a function key is pressed ADDER will output an integer 
equal to two plus the number of the key . The programming 
statements to accomplish this would be: 

ADDER : **F : ADDC ; CDefines ADDER as an instance of 

the function, ADDC.} 

SEND FIXC2) TO < 2> ADDER ; tPuts the value, 2, on 
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ADDER’S constant input.} 
CONNECT FKEYS< 1> :< 1> ADDER ; ^Connects function key 

output to ADDER’S 
trigger input.} 

CONNECT ADDER< 1> : C input and name of the element 

to which the sum is to be 
sent} 

CRoutes output to desired 
location . } . 

Available functions include mathematical operations, 
Boolean logic operations, routing functions, and data 
conversion . 

3.3. SPLINE 

Figure IB presents an overview of SPLINE, showing the 
interrelationship of major hardware and software 
components. Items marked by dashed lines have not yet 
been implemented, but provision far future addition of 
these features has been made within the program 
architecture and source coding. All programming on the 
PDP-11 was performed in FORTRAN language and all 

programming on the PS 300 was in the language described 
above. A more detailed description of the programs is 

contained in appendix A and the full FORTRAN and PS 300 
program listings are contained in appendices B and C, 
respectively . 

Program execution will vary slightly depending on 
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Fig. 16. Overview of SPLINE. 



42 



whether it is desired to perform a freeform ah initio 
design or to represent a previously known surface with a 
B-spline. In respect to this distinction, Wu, Abel and 
Greenberg C21D distinguish between shape design and shape 
representation . They identify shape design as having a 
primary concern of conceiving, displaying and modifying 
complex forms. Shape representation, on the other hand, 
is characterized by a need to describe or approximate 
existing surfaces for the purpose of graphical 
manipulation to enhance visualization or for geometric 
calculation to permit quantification. SPLINE will allow 
both approaches, shape design (freeform surface design) 
and shape representation (modeling a B-spline surface to 
approximate a previously known surface) . 

If it is desired to represent a known surface, data 
points on that surface may be input through a two or 
three-dimensional digitizer or by direct creation of an 
appropriately formatted data file. The program will allow 
screen display of these known points for direct visual 
comparison with the B-spline surface during modeling. In 
either case, the user must select the number and initial 
relative position of the control vertices (by designating 
the size and shape of a planar matrix of points), as well 
as the order of the surface in respect to each of the 
parametric directions. Provision has been made for future 



addition of an inversion algorithm 



which would calculate 
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initial control vertex locations to approximate known 
surface data . 

The control vertices, the B-spline surface and known 
surface points may be displayed on the screen in any 
combination through function key menu selections, and can 
be rotated, translated, or scaled in real time through use 
of the control dials. The B-spline surface is formed and 
modified by moving the control vertices about the screen. 
Uertex movement is communicated to the FORTRAN program, 
which recalculates the B-spline surface and sends a new 
vector list to the PS 300. 

Uertices are selected for movement by using the 
stylus on the data tablet to locate a cross hair on the 
screen over the desired vertex for picking. Individual 
vertices or entire rows or columns of vertices may be 
selected for movement. Movement is defined by 
manipulation of labeled control dials. 

Through the use of z-clipping, narrow sections of the 
model may be viewed without display of the remainder of 
the structures. This enables very accurate comparisons of 
any combination of display elements to be made very 
easi ly . 

B-spline surfaces are displayed as nets of constant 
parameter lines. Provision has been made for future 
implementation of orthogonal plane intersection mapping, 
which would enable display of more traditional 
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representations, such as waterlines and transverse section 
lines, yet based on the B-spline surface. 

3.4. Operating Sustem Limitations 

At this time the functionality of the program is 
severely limited by uhat appears to be a problem in the 
operating system or FORTRAN compiler of the PDP-11 
computer . 

During programming, it was noted that previously 
trouble-free subroutines would often fail after adding 
seemingly unrelated lines of FORTRAN program. This was 
found to be caused by variables changing value during 
program execution, but not under program control. Efforts 
to eliminate or program around this problem were not 
successful . Discussion with the flechanical Engineering 
Department computer system manager and with Professor 
Ganter of the thesis committee revealed that the PDP-11 
suffers from known operating system and compiler flaws 
which cause addressing problems resulting in similar 
symptoms . 

While the exact manifestation of the problem varied 
with the state of program development, the current 
symptoms include failure of the host-resident input/output 
library subroutines CPSIO’s) after execution of the 
subroutine which calculates the B-spline surface. This 
precludes display of the B-spline surface. A slightly 
earlier version of the program will allow the B-spline 
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surface to be calculated and displayed, but at an 
indeterminate time afterward will refuse to answer further 
calls to the communication subroutines. This usually 
prevents modification of the B-spline surface once 
displayed . 

Troubleshooting efforts have included careful review 
of the program for any definition of array elements 
outside of dimensioned array size, successful operation of 
each subroutine independently, and attempts to reestablish 
communication immediately after failure. The problem was 
further clouded by consultation with the Evans and 
Sutherland Software Technical Support Department CEEJ, 
which indicated that use of the PSIO subroutines has been 
discontinued in favor of more reliable communication 
packages and that there are documented and undocumented 
flaws in the preliminary PS 300 operating systems such as 
the one used for development of SPLINE. Nevertheless, as 
will be shown in chapter 5, the project results are 
sufficiently successful to demonstrate the viability of 
such a computer program . 



CHAPTER 4 



OPERATING INSTRUCTIONS 

4.1. Preparing Known Surface Data 

If it is desired to display known points from a 
surface for purposes of modeling a B-spline representation 
of the known surface, this data must be properly formatted 
before entering SPLINE. The data will automatically be 
correctly formatted by the digitizing program if it is 
prepared through use of the two-dimensional digitizing 
program, DIGIT, or the three-dimensional digitizing 
program, DIGITS. However, the first parameter in the 
global header may need to be adjusted Csee below) . 

The PREPS input format is used for the preparation of 
known surface data. PREPS is a data pre-processor used at 
the University of Washington for conversion of data 
describing the geometry of an object into a PS 300 syntax- 
compatible data file. The input file format for PREPS has 
been adopted as the standard for data input to SPLINE in 
order to facilitate the use of SPLINE with other PREPS- 
compatible programs, such as DIGIT and DIGITS. 

The PREPS input file format is detailed in Calkins’ 
laboratory book CSS) . A summary follows. 

The data file format prescribes a global header 
followed by a series of subheaders, each with an 
associated set of substructure data. SPLINE will allow 
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the display of individual data points represented in the 
data File, as well as lines sequentially connecting the 
data points within a substructure. Division of a 
structure into substructures will provide the most 
meaningful display if the substructures correspond to 
logical divisions in the actual structure. A common 
choice is the representation of profile lines or section 
contours as individual substructures. 

4.1.1 The global header. The global header is the 
first group of data in the PREPS-format data file and 
contains information affecting the display of all 
substructures. SPLINE reads the data with a list-directed 
format, so the exact format of the data in the file is not 
important. However, the correct sequence must be 

maintained . 

The global header contains five parameters: 

1. ZNORH the normalization factor for all 

data point coordinate values. 
This parameter should be set to 
the largest absolute value of any 
x, y, or z coordinate. 

E. INSTRT not used by SPLINE. Should be 

set to 1 . 

3. XEGDY global translation of all data in 

the input file along the x axis. 
SPLINE assumes that the 
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normalized data runs from 0 to 1 . 
If the origin is not actually at 
one end of the structure, use 
this to move the structure 
accordingly for the most 
convenient display . 

4. YBODY the same as XBQDY, for the y 

coordinate . 

5. 2BQ0Y the same as XBQDY, for the z 

coordinate . 

4.1.8. Substructure headers . Each set of 

substructure data is comprised of seven subheader 

parameters immediately followed by the substructure data. 
The subheader parameters are: 

1 . NPTS the number of data points in the 

substructure . 

8. NPTS8 the number of repetitions of this 

substructure to be displayed. 
Normally set to 1 . If several 

identical substructures are to be 
displayed, this can be set to the 
number or repetitions and ICOQRD 
to the coordinate which vanes 
between repetitions. The series 
of varied coordinate values is 



then entered for that coordinate 



in the data, one per substructure 



I COORD 



XTRAN 



YTRAN 



ZTRAN 



I REF 1 



repetition. Examples appear in 

C23H . 

the coordinate which remains 
constant within a substructure, 
or varies between repeated 

substructures . 

ICOORD “1, x coordinate 
- 2, y coordinate 
”3, z coordinate 
“ 0, no repetition, 
translation of the substructure, 
in non-normalized coordinates, 
along the x axis. 

the same as XTRAN, for the y 

axis . 

the same as XTRAN, for the z 

axis . 

the plane about which mirror 
image reflection occurs. 

IREF1 - 1, yz plane 
“ 2, xz plane 

- 3, xy plane 

- 0, no reflection. 

For modeling purposes, it is 
usually best to select no 
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reflection, as discussed in 
chapter 5 . 

4.1.3. Substructure data . Each subheader is 

immediately followed by the corresponding substructure 
data. The data is not paired for each point, but is 
listed with all x coordinates first, all y coordinates 
next, and all z coordinates last. If ICODRD is not set to 
0, only one value is entered for the constant coordinate, 
except as noted under NPTS2, above. 

ft sample data file is presented below, including a 
global header and subheader with substructure data . 
Comments contained in brackets are not part of the data 
f ile . 



3. 

10 0 
6 11 

0 0 0 

0 . 585E+00 
“0 . E55E+00 
-0.157E+00 
-0.962E+00 
-0.4BEE+00 
7 11 

0 0 0 



cznorh: 

0 CINSTRT, no x, y, z translation.] 
C First subheader; 6 points, one 
substructure, x coordinate 

constant for substructure.] 

0 CNo translation, no reflection.] 
CConstant x coordinate.] 



-0.202E+00 -0 . 191E+00 

-0 . B61E-01 -0.0 CG y 

-0.711E+00 -0 . 5B6E+00 

-0.370E+00 -0 . 362E+00 CG z 

CBegin second substructure 
0 



values . ] 



values . ] 
header . ] 
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0.256E+01 CBegin second substructure data ....1 

4.2. Preparing the Computer Sustem 

In order to run SPLINE, both the PS 300 computer and 
the host PDP-11 must be operational. Before logging on to 
tbe host, it is suggested that the PS 300 be readied. 

SPLINE requires P5 Uersion 8 firmware on the PS 300. 
The version can be checked by removing and inspecting the 
disk installed in the lower right-hand corner of the 

PS 300 cabinet. If it is necessary to reload the 
firmware, use the fallowing sequence: 

1 . Turn an the power switch located on the 

display screen. 

2. Turn off the power switch located on the 

upper right-hand corner of the PS 300 
cabinet . 

3. Insert the P5 Uersion 8 disk. 

4. Turn on the cabinet power switch. 

5. Wait for the DUAL LINE READY message on the 
display screen, which should appear in 
about five minutes. 

It will be necessary to log on to the PDP-11 in the 
normal manner in order to run SPLINE. 

4.3. SPLINE Operation Instructions 

After lagging an to the PDP-11 computer, execute the 
program by entering ’’RUN SPLINE.” A series of 
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’’Working...” statements will begin to appear on the screen 
as the PS 300 ’ s program File is read From the disk and 
passed From the host computer to the PS 300 under FORTRAN 
program control . 

4.3.1. Initial Formatting selections. AFter loading 
the PS 300 program, the LED indicators on the PS 300 
keyboard and control dials uill contain illuminated labels 
and a series oF questions uill be presented on the PDP-11 
terminal screen. 

The initial question is, ”Do you uish to display a 
known surface From a PREPS-Format File?” Answer ”Y” or 
”N”, For yes or no. IF a ”Y” is entered, the user uill be 
prompted For the name oF the data File containing the 
PREPS-Formatted data. See ’’Preparing Known Surface Data” 
above . 

The next series oF questions asks the user to speciFy 
the size oF the initial control vertex matrix Ci.e. 8,5) 
and the length-to-uidth ratio oF the matrix. The maximum 
matrix dimensions uhich can be selected are C10,B). The 
minimum dimensions are (5,5). The initial control 
vertices uill be positioned in the xy plane at z=0, uith 
the user’s First response representing the number oF 
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the x direction Clength) and that in the y direction 
Ciuidth) . Figure 17 illustrates the initial control vertex 
matrix resulting From selecting a matrix size of CB,5D 
with a length-to-width ratio of 3. 

The user is finally prompted to select the order of 
the B-spline in the longitudinal (u) and transverse Cw) 
directions. Third or fourth order B-splines may be 
selected. AFter entering the selections, small cubes 
representing the control vertices will begin to appear on 
the PS 300 screen . 




4.3.5. Enabling the PS 300 keuboard . After all 



initial formatting selections have been made Cin response 
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to display screen prompts! , the Following message will 
appear on the screen: 

*»*»*»'l*'i»«»'***»*»***»*«»4t***'l»*>*>*'l»*»**4t*»'l»'»*»»*'**l'* 

** PRESS <SHIFT-LINE/LOCAL> ** 

** ON THE PS-300 KEYBOARD ** 

** FOLLOWED BY <TERN> ** 

The LINE/LOCAL and TERN keys are located in the keypad at 
the left side of the PS 300 keyboard. Press the 
LINE/LOCAL key while holding down the SHIFT key. This 
enables the twelve Function keys at the top of the 
keyboard. Press the TERN key if it is desired to remove 
the text display From the screen. 

4.3.3. Control dials. The control dials are labeled 
by LED indicators and have the Following Functions: 

GLOBAL X moves the model in the x direction in 
the screen coordinate system Cleft and 
right! . 

GLOBAL Y moves the model in the y direction in 
the screen coordinate system Cup and 
down! . 

GLOBAL 2 moves the model in the z direction in 
the screen coordinate system Cm and 
out of the screen! . As the model 
moves into the screen Caway From the 
viewer! past the z=0 plane, depth 
cueing will cause portions of the 
model to Fade, creating the illusion 
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4.3. 

f unction 

functions 

available 



of depth. Depth cueing can be 
disabled by moving the entire model 
forward of the z-0 plane Ctoward the 
viewer) with the GLOBAL 2 dial. 

SCALING scales the entire display, allowing 
the user to examine details of the 
display more closely . 

ROTATE X rotates the model about the screen’s x 
(horizontal) axis. 

ROTATE Y rotates the model about the screen’s y 
(vertical) axis. 

ROTATE 2 rotates the model about the screen’s z 
axis (turning it within the plane of 
the screen) . 

4. Ham menu. The LED labels above the keyboard 

keys present the mam menu, from which SPLINE 
are selected. The following options are 



FI UERTICES is a three-way toggle controlling the 

display of the control vertices. 
Uertices are displayed on the screen 
as small cubes. Pressing the FI key 
selects one of the following options, 
in sequence: 

(1) vertices displayed, 

(2) vertices displayed with straight 
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lines connecting them in matrix 
position order, 

C3) vertices not displayed. 

F S B-SPLINE toggles display of the B-spline 

surface on and off . When the F2 key 
is pressed, the B-SPLINE label above 
the key will blink while the host 
computer calculates a vector list for 
the B-spline surface and sends it to 
the PS 300. When the surface is 
displayed, the label will cease 
blinking . 

F3 KNOWN PT is a three-way toggle controlling the 

display of known surface points 
contained in a PREPS-format file. 
Points are displayed as asterisks 
C*’s). Pressing the F3 key selects 
one of the following options, in 
sequence : 

C1D known points are displayed, 

CB) known points are displayed and 
points within substructures are 
connected in sequence by straight 
lines, 

C3) known points are not displayed. 

F4 



INTERACT toggles SPLINE into and out of the 
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interactive modeling mode, giving the 
user the ability to manipulate control 
point positions. When the F4 key is 
pressed, several function key and 
control dial labels change to indicate 
new assigned functions or availability 
of previously assigned functions. 
Interactive modeling is discussed in 
section *±.3.5 below. 

F5 CLIP toggles z-clipping on and off. When 

z-clipping is enabled, only portions 
of the model lying within a narrow 
band parallel to the plane of the 
screen can be viewed. This allows 
slices or sections of the model to be 
viewed without confusion caused by 
lines in front of and behind the 
section being examined. When the F5 
key is pressed, the screen will be 
blanked and a CLIPPING label will 
appear above the eighth control dial. 
Turning the dial slightly will cause a 
slice of the model to be displayed. 
Adjusting the CLIPPING dial will vary 
the depth of the slice which is 
displayed. By turning the GLOBAL 2 
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dial, it is possible to move the model 
through the clipping window, allowing 
the viewing of any portion of the 
model . Note that it is possible to 
move the model entirely outside of the 
clipping window, in which case it 
cannot be viewed. 

F5 REFLECT has not been implemented. When 

implemented, it will toggle on and off 
the display of a reflection of the 
model about its XZ plane. This label 
will not be displayed during 
interactive modeling, indicating that 
the reflection function is not 
available at that time. CThis 
prevents the updating of the reflected 
surface from slowing down the 
computer’s response to interactive 
inputs during modeling.) 

F7 XY has not been implemented. When 

implemented, this key will toggle on 
and off the display of intersections 
of the XY plane Cbody-fixed 
coordinates) with the B-spline 
surface. On a ship model, this 
corresponds to the display of 
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waterlines. This Function is not 
available during interactive modeling. 

F3 XZ has not been implemented. When 

implemented, it will Function 
similarly to the F7 key, described 

above, and will display lines 
corresponding to a ship’s buttocks 
lines . 

F9 YZ has not been implemented. When 

implemented, it will Function 
similarly to the F7 key, described 

above, and will display lines 
corresponding to transverse sections 
oF a ship hull . 

F10 I/O selects the input/output menu. 

Pressing this key causes the labels 
above the Function keys to display a 
menu oF input and output Functions 
described in section 4.3.6, below. 

Fll RESET is a three-way toggle which resets the 

displayed model to each one oF three 

orthogonal views, in sequence. 
Pressing the Fll key resets the scale 
Factor and global translations to 
those oF the initial display and 

orients the XY, XZ, or YZ plane oF the 
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model parallel to the screen. 
Repeatedly pressing the key will cycle 
the display through the three 
or ientations . This Function key is 
particularly useful when orienting the 
model for display of slices with the 
CLIP function. 

FIE QUIT stops execution of the program and 

initializes the PS 300, clearing the 
screen. If it is desired to use the 
keyboard to communicate with the host 
computer at this point, press the TERN 
key to return the screen to a text 
display mode, then press the 
LINE/LOCAL key, followed by RETURN. 
The host computer’s prompt should now 
be visible. 

4.3.5. Interactive modeling. The interactive 
modeling mode is entered and exited by pressing the 
INTERACT CF43 key. This mode is the heart of the computer 
program, allowing the user to select and relocate control 
vertices in order to mold a B-spline surface. 

When the F4 key is pressed, the labels above the 
three global translation control dials CGLOBAL X, 
GLOBAL Y, and GLOBAL 2) change to UERTEX X, UERTEX Y, and 
UERTEX 2 and begin to blink. These dials can now be used 
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to move the selected vertex in the indicated direction in 
the body-f ixed coordinate system. If it is necessary to 
move the entire body, the global translation functions can 
be restored by toggling the interactive modeling function 
off with the F4 key. The blinking labels will remind the 
user, while in the interactive modeling mode, that dial 
movement will relocate vertices rather than translate the 
entire model, in order to prevent inadvertent movement of 
vertices . 

Individual control vertices are selected for movement 
with the stylus and data tablet. When the tip of the 
stylus is placed near the surface of the data tablet, a 
tracking cross will appear on the screen. By moving the 
stylus lightly over the data tablet surface, the tracking 
cross can be positioned over the desired control vertex. 
The vertex is then selected by pressing the stylus 
momentarily against the data tablet surface while keeping 
the cursor positioned over the control vertex. This 
activates a switch in the tip of the stylus, causing the 
data tablet to report the stylus position, leading to 
identification of the selected vertex. The vertex can 
then be moved into the desired position with the vertex 
translation dials. If it is necessary to reorient the 
model in order to obtain a better view or reduce ambiguity 
for the selection or location of control vertices, this 
may be accomplished with the rotation dials or RESET 
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function key, which remain operational while in the 
interactive modeling mode. 

When a control vertex is relocated, the movement is 
reported in increments to the host computer, which 
recalculates vector lists defining the B-spline surface 
and the net of lines connecting the vertices. If the 
vertex-connecting lines or B-spline surface are displayed 
while in the interactive modeling mode, some delay will be 
noted between the time of vertex movement completion and 
the corresponding change in the model. This delay is due 
to the time required for the host computer to communicate 
the updated vector lists to the PS 300. The delay will be 
about two to three seconds for update of the vertex- 
connecting lines and about ten seconds for update of the 
B-spline surface. 

When the interactive modeling mode is activated, the 
label above the F7 key will read ’’PICK ROW.” Pressing 
this key will alter the control vertex selection function 
such that selections made with the stylus will affect an 
entire row or column of vertices at a time. The column or 
row is selected similarly to the selection of an 
individual vertex, except that the cursor is placed over 
one of the vertex-connecting lines in the row or column, 
instead of over the individual vertex. Turning the vertex 
translation dials will then move the selected row or 



column 



of vertices in the desired direction 



while 
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retaining the orientation of the individual points within 
the row or column to each other . The user can return to 
selection of individual control vertices by pressing the 
’’PICK PT” CF0} key. 

After moving a vertex it will occassional ly be 
noticed that the vertex-connecting lines will not quite 
reflect the new position of the vertex. Since the vertex 
movement is reported to the host computer incrementally, 
it is possible that a slight additional movement will 
trigger the update. At other times, a flaw in the 
communications between the host computer and the PS 300 
causes an incremental change to be reported to the host as 
a very large or indecipherable movement. In that case, 
the distorted movement value is disregarded by the host 
program, which results in loss of the true incremental 
change. While there is no way to reconcile the difference 
between the control vertex position as viewed on the 
screen and that held in the host computer’s memory, it 
should be noted that the intersection of the vertex- 
connecting lines, and not the vertex cube, will indicate 
the true position upon which the B-spline surface will be 
calculated . 

4.3.6 Input/output menu. The input/output menu is 
displayed in the function key labels when the I/O CF10) 
key is pressed. The menu functions have not yet been 
implemented Cwith the exception of the FIE key}, but when 
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implemented will offer 
FI PRT UERT 



F4 PLI SCRN 



FB SAU UERT 



FB LD UERT 



FIB EXIT I/O 



the following options: 

prints a list of the current 
vertex locations and of the 13- 
spline surface order selections, 
plots a drawing of the model, as 
oriented and displayed on the 
screen, on the laser printer, 
saves the current vertex location 
and surface order information to 
a disk file. 

loads vertex location and surface 
order information from a disk 
file created with the SAU UERT 
function . 

returns the user to the main 



menu . 



CHAPTER 5 



MODEL I NG TECHNIQUES AND CASE STUDIES 

The B-spline modeling of sculptured surfaces greatly 
reduces the requirement for the more traditional 
techniques of surface development and representation 
reviewed in chapter 1. As mas shown there, use of the 
traditional techniques requires that engineering and 
drafting skills be accompanied by an artistic touch in 
order to visualize and fair the surface in three 
dimensions. This ability is developed by the naval 
architect, and others who work often uith sculptured 
surfaces, through years of experience and practice. 

□n the other hand, while the B-spline techniques 
deemphasize the artistic requirement, it will still be 
necessary for the designer to develop an intuition for the 
placement of control vertices which cannot be entirely 
reduced to a set of mathematical rules without sacrificing 
the control and flexibility of interactive modeling. The 
heuristic observations and case studies which follow are 
offered in order to facilitate the rapid acquisition of 
that intuition. 

5.1. Modeling Tips and Techniques, 

5.1.1. Known data. If the surface design is not to 
be developed ab initio , the preparation of a data file of 
known surface points will be essential. The selection of 
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such points for inclusion in this data File requires some 
consideration and advance planning. Primarily, the goal 
is to select sufficient points to accurately represent the 
form and boundaries of the surface, yet not so many as to 
create unnecessary clutter when the known data is 
displayed . 

Two other considerations should be kept in mind 
during data point selection. First, the practicalities of 
displaying and visualizing the two surfaces which are to 
be matched Cthe known surface and the B-spline surface or 
control graph) dictate that most of the modeling be 
performed while viewing two-dimensional sections Cor very 
thin three-dimensional sections) . This generally means 
locating control vertices nearly in the planes of known 
substructure data, such as the planes of a ship’s 
transverse sections. Greater flexibility will be 
available in locating the control vertices for maximum 
surface control if the vertices’ general locations are 
considered before finalizing the file of known surface 
data. Sections 5.1.3 and 5.1.4 will address the problem 
of vertex placement. 

Secondly, if a symmetrical body such as a ship hull, 
automobile body, or airplane fuselage is being modeled, it 
will generally be easier if only one half of the structure 
is displayed, both of the known surface data and of the B- 
spline surface and control graph. Modeling two identical 
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halves requires twice the manipulation as modeling a 
single half and results in much more difficulty 
interpreting the wire frame model due to ambiguities 
between lines in the foreground and those in the 
background. It is also virtually impossible to exactly 
duplicate the control graph for the halves of the body 
through a purely visual comparison. The only drawback to 
modeling a single half of a symmetrical body is that it is 
more difficult to ensure fairness at the juncture of the 
two halves. This problem can be remedied through use of 
the reflection display function, when that function is 
implemented . 

5.1.5. Surface order . SPLINE offers the choice of 
two surface orders: third order or fourth order. The 
order may be specified independently in the a and w 
directions. As noted in chapter 5, a fourth order curve 
produces the smoothest lines, but at the expense of some 
local control and the ability to easily effect sudden 
changes in curvature. It should also be remembered that a 
fourth order curve requires three coincident vertices Cor 
rows of vertices J to create a knuckle or chine, while the 
third order curve requires only two. 

The third order surface will normally produce 
satisfactory results with the least amount of effort. The 
fourth order curve should be considered where smoothness 
of lines is an overriding constraint, such as in the 



66 



longitudinal lines of some ships and aircraft. 

5.1.3. The control graph. Selection of the 
dimensions of the control graph Ci.e., the number of 
vertices contained} is best accomplished by considering 
the locations where control points will be required for 
sufficient control of the surface. The most obvious 
considerations include: 

Cl} Uertices will be required to define the 
boundaries Cends and edges) of the surface. 

C2) A vertex or row of vertices will generally 
be required at each maximum or minimum 
point . 

C3) Particularly sharp curvatures or sudden 
changes in curvature will generally require 
two or more control vertices to achieve . 
Inflection points will require sufficient 
control vertices on each side of the point 
of inflection to establish curvature in 
opposite directions. 

C4) Knuckles or hard chines will require two 
coincident rows of vertices in a third 
order surface or three coincident rows of 
vertices in a fourth order surface. 

C5) While there is no constraint imposed by the 
B-spline algorithm nor by the computer 
program which would dictate that control 
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vertices within a row or column of the 
control graph remain in a plane, the 
practicalities of modeling usually make 
such placement the wisest choice, as 
discussed above. In most cases, comparison 
with known surface points will be most 
efficiently accomplished if the rows Cor 
columns) of the control graph are parallel 
to the known-surface substructures. It is 
not practical to constrain both the rows 
and the columns of the control graph to 
planes, but with one or the other so 
constrained, some effort should be made to 
keep the other in as close to a planar 
placement as possible. This will be 
emphasized by the case study in 
section 5.3. 

After the number of control vertices has been 
established, required manipulation of the vertices can be 
minimized by selecting a length-to-width ratio for the 
initial control graph which approximates that of the body 
to be modeled. This will be illustrated in the case 
studies which follow. 

5.1.H. Control vertex placement. The approach to 
control vertex placement which was found to be most useful 



is to first locate the rows or columns of the control 
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graph in the desired plane and then to use the CLIP 
Function to compare an individual row or column of the 
control graph with the nearest substructure of the known 
surface data. In this way , most of the interactive 
modeling is performed in the more readily visualized two- 
dimensional realm. Placing vertices is primarily a matter 
of keeping in mind the general properties of the B-spline, 
as innumerated in chapter 2. Those properties are briefly 
reviewed below. 

Cl) The B-spline curve does not interpret the 
control polygon, but is contained within 
its convex hull, taken K vertices at a 
time. This dictates locating the control 
vertices slightly outside of the curvature 
of the desired surface. 

C2) The B-spline surface interpolates the 
corners of the control graph. 

C3) While modeling is most efficiently 
performed in two dimensions, it must be 
kept in mind that the influence of each 
control vertex is three-dimensional . This 
will occassional ly require moving a vertex 
from its apparently proper position in two 
dimensions to accommodate some sudden 
change in the surface in the third 
dimension, or to avoid distorting the 
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surface by locating a vertex far from the 
plane of the remainder of its row or column 
in the control graph. 

C4) The slope of the curve at the end points of 
the control polygon is determined by a 
straight line connecting the end point with 
the next point in the polygon. Again, the 
effect of vertices in the third dimension 
may cause unexpected deviations. 

C5) A straight line is produced by K collinear 
vertices . 

CB1 A curve tangent to the control polygon is 
produce by K—l collinear vertices. 

C7D A knuckle is produce by K—l coincident 
vertices . 

5.1.5. Menu selections. The selection of display 
options from the menu can have a significant effect on the 
efficiency with which interactive modeling is performed. 
The requirement for displaying the control vertices is 
obvious, but it is useful to display the lines connecting 
the vertices of the control graph as well. Display of the 
connecting lines prevents confusion concerning the 
identity of a vertex in respect to its position in the 
control graph matrix. If two paints in a control polygon 
are inadvertently interchanged, the result can be dramatic 
anomalies in the B-spline surface. 
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The CLIP function should be used as much as possible, 
since accurate depth perception is difficult with a uire 
frame model, even with depth cueing. This function is 
best disabled, however, when rotating the model. The 
RESET function is particularly useful for orienting the 
model for viewing sections in orthogonal planes with 
clipping . 

While the object of the modeling is to mold a 8- 
spline surface, most of the initial modeling is best 
performed without constant display of the surface. 
Display of the B-spline surface adds a profusion of lines 
to the already crowded screen display and slows response 
time considerably, since the host computer must 
communicate a new vector list for the surface to the 
PS 300 each time a vertex is moved. Due to the highly 
intuitive nature of B-spline manipulation, it is passible 
to achieve surprisingly accurate results by making the 
initial placement of control vertices without viewing the 
resulting surface at all. Due to the current limitations 
of the operating system, the images presented in the three 
case studies which follow were generated in this manner. 
Once the B-spline surface was initially displayed, further 
manipulation was not passible, as discussed in 
section 3.4. As will be seen, little further manipulation 
would be required for very accurate surface 
representation . 
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5.I.G. General □rocedure. The basic procedure for 
performing surface modeling with SPLINE can be summed up 
in a few steps. 

Cl} Preplan the initial placement of control 
graph rows and columns and select 
corresponding control graph dimensions. 

C2} Do as much control vertex placement as 
possible by moving entire rows or columns 
of the control graph at a time. 

C3} Use the CLIP function while modeling 

substructures in two dimensions, initially 
without display of the B-spline surface. 

C4} After initial control vertex placement, 
view the entire control graph in three 

dimensions. Ensure that the individual 
control graph rows and columns do not have 
undesired extremes. 

C5} Display the B-spline surface. 

C6) flake minor adjustments to the control graph 
as needed . 

5.2. Case Studu ttl; Modeling the Australia II 

ThB first case study is a model of the hull of the 

sailing yacht, Australia 11, excluding the famous winged 

keel. This case study illustrates the general approach to 
modeling a B-spline surface to an existing known surface 
or to preliminary design lines. The known surface data 
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IB. Lines drawing of the Australia II. Source: 
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was extracted from a lines drawing (figure IB} contained 
in Stannard’s biography of Ben Lexcen C241. 

Before entering SPLINE, the known surface data was 
put into a PREPS-f ormat file by digitizing it from the 
lines drawing with the two-dimensional digitizing program, 
DIGIT. Transverse sections were designated as 
substructures with about six to ten points each. Since 
the body plan presents the fore and aft section profiles 
on opposite sides of the xz plane, it was necessary to 
manually modify the data file created by DIGIT to place 
all half-section profiles on the same side of the plane. 
This was a simple matter of deleting minus signs from the 
data . 

The final step before entering SPLINE was to 
determine the required control graph dimensions. It was 
apparent that modeling would be performed most 
conveniently in two-dimensional sections corresponding to 
substructures of the digitized data. Since the body shape 
was relatively uncomplicated in the body plan, it was 
decided to use only five rows in the control graph. 

The longitudinal lines of the yacht required greater 
consideration. It will be useful to view figure IS, which 
shows the selected control graph as squares connected by 
net lines and known surface points as unconnected 
asterisks, while examining the selection process which led 
to this control graph. The most easily selected locations 
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Fig 19. Initial control graph for Australia II, 
after relocating columns to known hull stations. 



for columns of the control graph were a station at each 
end of the hull to, define those boundaries, and a station 
to define the maximum point of the keel. Since the shape 
of the forward sections is relatively uniform, two 
additional stations were expected to be sufficient to 
define the forebody. 

The after half of the hull presented more interesting 
problems. The stern is not a U-shaped section, as most of 
the hull, but a U-shaped section. The stern also has a 
short portion which is somewhat elevated and very narrow, 
with nearly parallel sides. A second station was selected 
at the stern to help maintain the shape and size in that 
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region . 

Between the keel and the stern, the transition from 
the U-section to the U-section is relatively abrupt, as 
can be seen in figure BO, which shows the two consecutive 
transverse sections. C Asterisks represent the known 
surface points.! A control graph column was selected for 




Fig. BO. Abrupt transition m the after hull of 
Austral 2 a II. 



each of the illustrated sections, to define the 
approximate beginning and end of the transition region. 
It was expected that additional control would be required 
to ensure that the U-shaped station did not unduly 
influence the hull forward of the transition, nor the U- 
shaped section influence the shape aft of the transition. 
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The proximity of the two control graph columns selected 
for the stern was thought to be sufficient to maintain the 
shape between the U-shaped station aft of the transition 
and those at the stern. An additional U— shaped station 
was selected adjacent to the one located Just forward of 
the transition, to minimize influence of the U— shaped 
stations forward of the transition. This made a final 
count of nine columns. 

Upon entering SPLINE, a control graph size of C9,5i 
was selected, based on the above planning. The basis 
functions were defined as third order in both directions, 
and an initial control graph length-to-width ratio of 
eleven was chosen to roughly correspond to that of the 




Fig. 81. Initial control graph for Australia 11. 
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hull. The resulting initial control graph can be seen in 
figure 21, along with the digitized known surface data. 

Interactive modeling of the yacht hull was begun by 
selecting the ROW PICK function and moving the control 
graph columns into the selected positions, as shown in 
figure 19. The RESET function was then used to orient the 
hull for a bow view and the display was scaled to fill the 
screen. The CLIP function was then activated and the 
clipping window adjusted to a depth which allowed viewing 
the control vertices and known points of a single 
transverse station at a time Cfigure 223. 



Individual transverse sections were modeled one at a 




Fig. 22. A station of the Australia II before 
modeling . 
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time by selecting and moving the control vertices 
individually and moving the hull through the clipping 
window with the GLOBAL 2 dial to view each new station. 
All Five control vertices were located coincidentally at 
the bou . Figure 23 shows a more typical station with 
relocated vertices. Notice that the vertices are located 
slightly outside of the curvature of the known surface, 
since the B-spline surface will be contained in the convex 
hull of the control graph, but will not interpolate it. 
This offset is exaggerated in the illustration. The 
small tabs extending from the control vertices are the 
control graph lines which extend beyond the clipping 
window to vertices in other control graph columns. 




Fig. 23. Control polygon for a U-shaped station. 
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Figure 24 illustrates the control polygon For a U- 
shaped section. The control vertices at the right are 
located above the bottom of the hull in order to remain 
outside of the curvature which occurs in the longitudinal 
direction Cperpendicular to the screen) . Three collmear 
vertices are used to define the flat bottom between the 
two right-hand points, as well as the slope of the bottom 
as it begins the upward curvature between the second and 
third collinear vertices. 




Fig. 24. Control polygon for a U-shaped section. 



The final control graph is pictured in figure 25. 
Notice that the net of lines is uniform in both the 
transverse and longitudinal directions. The resulting 13- 
spline surface is shown in figure 26 with the known 
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Fig. 26. 



B-spiine modal of Australia II. 
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surface paints Far comparison. The match is generally 
quite good, but an additional control graph column is 
needed near the boui to hold the keel line downward. A 
closer view of the transition from U-shaped to U-shaped 
sections is shown in figure 27, which depicts the B-splme 
surface and its associated control graph. 




Fig. 27. B-splme model of Australia II stern. 



The next step in modeling the hull would be to again 
view the transverse stations one at a time and make minor 
adjustments in the control graph. Figure 28 compares the 
control vertices (squares!, known surface (asterisks! and 
B-splme surface (solid line! for a U-shaped station. It 
can be seen that only very minor adjustments are 



indicated . 




Fig. 20. 



B-spline surface at a U-shaped station. 
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Figure S3 shows the B-spline surface achieved at a U- 
shaped section and illustrates one cf the pitfalls cf 
considering primarily tujc-dimensional sections while 
modeling a three-dimensional surface. The three collinear 
vertices at the bottom would dictate a straight line in 
the right-hand segment of the B-spline curve in two 
dimensions. However, the surface is pulled downward at 
the right-hand edge by the effects of the vertices at 
stations located further forward. A considerable upward 
adjustment of the right-hand control vertex may be 
necessary to overcome this pull from other stations. 

5.3. Case Studu #5: Modeling a Porsche 344 

The second case study involves very similar 
considerations to the first and will be presented more 
briefly. The body of a Porsche 344 automobile was 
modeled, reemphasizing the techniques discussed previously 
and illustrating the versatility of the B-spline and of 
SPLINE. 

The known surface data for the Porsche was 
approximated by digitizing points from a lines drawing and 
profile contained in advertising literature. Once again, 
DIGIT was used to create the data file. The data for half 
of the symmetrical body was presented to SPLINE without 
reflection . 

The shape of the surface was considered in order to 
estimate the required dimensions of the control graph. 
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Fig . 31 . 



Relocated control graph columns for Porsche 344 . 
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Figures 30 and 31 illustrate the selection of transverse 
stations For modeling, similar to those chosen For the 
previous case study. In this case, the control graph is 
presented in an edge view. Looking First at the 
longitudinal lines, it was determined that ten control 
graph columns would be required. A station was selected 
at each end, another near the Front Cto avoid the problem 
encountered in modeling the bow oF the Australia ID, and 
a station near the middle at the point oF maximum height 
in the proFile. Two stations were determined to be 
necessary to achieve the small-radius curvature at the 
base oF the windshield, and two more at the top. Finally, 
two intermediate stations were positioned in the rear halF 
oF the body. IF an eleventh station had been available, 
it would have been place midway along the hood oF the 
automobile . 

With consideration to the shape oF the transverse 
sections, seven control graph rows were deemed necessary. 
Figure 32 depicts a section across the hood being modeled 
and illustrates the reason For selection oF seven control 
graph rows. Since the sharp curvatures and sudden 
transitions in the body dictated use oF a B-spline which 
was third order in both dimensions, three vertices were 
required to achieve the nearly Flat hood and rooF oF the 
automobile, with an additional vertex at the hard chine 



between the hood and Fender. 



CThe brighter vertex seen 
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there is actually two coincident vertices.) Three 
vertices mold the curve at the side of the automobile. 




Fig. 32. Modeling a forebody section of a Porsche 944. 



As detailed previously in case study #1, interactive 
modeling was performed by first moving control graph 
columns to the selected stations, then modeling individual 
transverse stations with clipping enabled. Examining the 
control vertices of figure 32 from right to left, the 
first three vertices are placed nearly linearly to give 
the nearly flat surface of the hood. The third and fourth 
vertices are placed at the same point to achieve the 
knuckle. The fifth vertex serves two functions. The 
control polygon line connecting this vertex with the 
vertices at the knuckle defines the slope of the curve as 
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it flows downward from the knuckle. The position of the 
fifth and sixth vertices define the maximum for the curve 
in the fender. Finally, the last vertex defines the lower 
boundary of the surface. 
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Fig. 33. Control graph and known surface for a 
Porsche 344. 

Figure 33 presents the final control graph after 
modeling the individual sections, along with the known 
surface data. While each individual station appeared to 
be adequately modeled, displaying the control graph lines 
Cfigure 343 revealed a problem. Near the midsection of 
the automobile, one vertex was placed much lower on the 
body's contour than the corresponding vertices in adjacent 
columns. While the placement had appeared adequate when 
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Fig. 34. Control graph for Porsche 944 with an 
inconsistency in one row at the midsection. 




Fig. 35. B-spline representation of a Porsche 944. 
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viewed in the transverse plane, it was now apparent that 
this vertex was sharply out of position when viewed as 
part of its control graph row. 

The effect of this placement can be seen in 
figure 35. Much of the apparent distortion in the body at 
the midsection is only an illusion caused by the manner in 
which the net of lines representing the B-spline surface 
is displayed. Some real distortion does also occur in 
this region due to the added pull of the misplaced vertex. 
This case study illustrates well the importance of 
considering and reviewing the control graph in three 
dimensions . 

Figures 35 and 37 present comparisons of the B-splme 




Fig. 35. B-splme surface with a hard chine. 
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Fig. 37. B-spline representation of a Porsche surface. 



surface to the control vertices and known surface points. 
Note that the two coincident vertices in figure 36 produce 
a hard chine, while placing the vertices slightly apart, 
as in figure 37, produces a small-radius curvature. 

The effect on this station of control vertices at 
other stations can again be seen in figure 36, where the 
surface of the hood is pulled up from the desired form by 
the control vertices at the next station, which are 
located higher for the windshield. Figure 37 also 
illustrates the remarkably accurate surface which can be 
obtained even without viewing the B-spline surface during 
initial control graph arrangement. 
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5.4. Case Studu #3: Freeform Modeling 

The final case study presents an example cf Freeform 
modeling, in this case of a Stetson hat. In ab initio 
design, the engineer often has a general idea of the form 
which he or she is trying to achieve, but is not 
constrained by a requirement to exactly match previously 
existing lines. For this reason, no attempt was made to 
compare this model to a preexisting design. 

In the absence of such previous form definition, it 
will often be beneficial for the designer to sketch a few 
basic lines to aid in visualizing the shape and the 
required control graph dimensions. Because of the 
difficulty in visualization of the more abstract Form, in 
this particular case it was also found to be useful to 
model both halves of the object, even though it is 
symmetrical . 

Aside from these initial considerations, the modeling 
approach is basically the same as that in the first two 
case studies. A third order B-spline was used with a 
seven by seven control vertex matrix, based on form 
considerations which will become apparent as the model 
development is reviewed. A length to width ratio of 1.3 
was selected. Figure 3B presents the initial control 
graph . 

Uiewing the planar control graph as being in the 
horizontal plane of the hat’s brim, vertices defining the 
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Fig. 39. 



Centre! graph reus and columns relocated. 
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sides of the crown were to require placement nearly 
vertically above one another. For this reason, the 
ROW PICK feature was used to relocate rows and columns of 
the control graph as shown in figure 39. Notice that the 
second column from the right exhibits the disparity 
between control vertex position and control graph line 
position which was discussed at the end of section 4.3.5. 
When moving the column of control vertices, one increment 
of the movement was distorted during communication between 
the two computers and thus ignored by the host computer’s 
program. In all subsequent operations, there will 
continue to be a disparity between the position of these 
vertices as displayed by the PS 300 and as used in the 
preparation of vector lists and calculation of surfaces by 
the host computer. The intersection of the control graph 
lines, and not the individual cubes, marks the vertex 
position upon which the B-spline surface will be based. 

After positioning individual vertices to define the 
rounded shape of the surface edges C figure 40), the 
vertices controlling the surface of the hat’s crown were 
pulled up from the surface. To achieve uniformity of the 
similar halves, the paints were first selected with the 
stylus while the control graph was in an orientation 
similar to that shown in figure 41, then the display was 
rotated to present an edge view Cfigure 42) before the 
selected vertices were moved. Also in the interest of 



Fig. 40. Surface edge boundaries established. 




Fig. 41. 



Elevated vertices define croun . 
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achieving uniformity, the modeling uias performed without 
use of the CLIP function, allowing direct comparisons to 
be made between foreground and background lines. Notice 
the depth cueing in the right rear quadrant of figure 41 
which gives the viewer a better perception of the display 
orientation . 




Fig. 42. Elevating the brim. 



The final step in forming the control graph was to 
elevate the first and last rows of the control graph to 
achieve the raised brim (figure 42} . It was passible to 
da this very easily in the ROW PICK mode, even though 
points in the rows had been moved around to accomplish the 
rounding of the edges of the brim. The entire row of 
vertices could be moved simultaneously without affecting 
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the relative position of the individual vertices in the 
row to each other . 




Fig. 43. Stetson hat with control graph. 



Figure 43 presents the final control graph of the hat 
with the resulting B-spline surface. Notice that uhile 
the vertices form a relatively rounded brim (figure 40), 
the resulting surface has distinct corners Cfigure 43). 
This effect is due to the fact that only the corner 
vertices of the control graph are interpolated by the 
surface. Since the adjacent vertices are not 
interpolated, the result is a sharp discontinuity at each 
of the four corners of the control graph. In order to 
present the desired smooth boundary, these four vertices 
must be relocated to lie on the intended edge of the 
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surface, rather than on a smooth curve connecting the 
adjacent vertices. 

While the hat perhaps belongs more correctly on 
Huckleberry Finn’s head than on a cowboy’s, it 
nevertheless illustrates well the great potential for 
freeform design with surface B-splines. Figure 44 
presents a section view of the same hat and control graph. 




Fig. 44. Section view of hat with control polygon. 




CHAPTER 5 



RECOMMENDATIONS AND CONCLUSIONS 

6.1. Hardware and Firmware Modifications 

Several modifications to the hardware and firmware of 
the system upon which SPLINE was developed would 
significantly enhance the program. At this time the 
University of Washington Mechanical Engineering 
Department’s PDP-11 computer is being replaced by a UAX- 
11/750. This replacement makes possible most of the 
enhancements discussed below. 

6.1.1. High-speed interface. It was initially 
planned that SPLINE would make use of rubber sheeting, the 
real-time modification of the B-spline surface display in 
response to control vertex manipulations. However, it 
soon became apparent that communication between the two 
computers, at 3600 baud, was far too slow to effect real- 
time response. 

Rogers and Satterfield C171 found calculation of the 
full B-spline algorithm too slow to allow surface dragging 
or rubber sheeting, so developed an abbreviated algorithm 
which took into effect the localized effect of a single 
control vertex movement. Using the present system, even 
the time for a full B-spline surface calculation was found 
to be insignificant compared to the time required to 
communicate a new vector list, however calculated, to the 
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PS 300. It takes approximately ten seconds for the host 
computer to send the vector list for the B-spline surface. 

The PS 300 manual C20D states, ’’This interface is 
best suited for infrequent communication of small amounts 
of data.” The interface simply is not designed to support 
a program such as SPLINE. Upgrading the system with Evans 
and Sutherland’s 5E kilobaud high-speed interface would 
provide a distinct improvement in the program operation. 

5.1.2. UAX host computer . The greatest potential 
benefit of interfacing the PS 300 with the newly installed 
UAX computer is avoidance of the PDP-11 operating system 
discrepancies. It is expected that SPLINE would require 
few modifications to compile on the UAX system. 

5.1.3. Graphics Support Routines. A totally 
different host-resident communications software package is 
available for the PS 300 on UAX systems than is used on 
the PDP-11. The Graphics Support Routines provide 
significantly more reliable communications between the 
host computer and the PS 300. Additionally, these 
routines perform some preprocessing of data in the host 
CPU, resulting in substantially faster communications. 
Upgrading the system with the Graphics Support Routines 
would require some revision of the manner in which SPLINE 
presents data for communication to the PS 300, but would 
provide a worthwhile enhancement. 
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6.S. Program Expansions 

SPLINE was designed to support many features which 
are not yet Fully implemented. For many of these 
features, the PS 300-resident portion of the program 
already includes the necessary statements to handle 
display and selection of the additional options. The 
programs are carefully documented to facilitate expansion. 

E.S.l. Input/output menu selections. The first 
program expansion which should be investigated is the 
implementation of all input/output menu selections. At 
the present time the function key presses are reported to 
the FORTRAN subroutine, 10. However only a shell of the 
subroutine exists and is totally nonfunctional. 

The most needed output ability is that to print the 
final control vertex positions on disk or paper. The 
primary purpose of the program is actually to obtain such 
a listing, the control vertex positions being an integral 
part of the B-spline algorithm which defines the modeled 
surface. Such a listing should include the order of the 
surface, the homogeneous coordinates, and the knot 
vectors. These values are readily available and 
identified within SPLINE. FTN Cthe FORTRAN portion of 
SPLINE) and interfacing with the printer should be 
relatively straightforward. 

The variable, KEY, in subroutine 10 is equal to the 
number of the function key which was pressed. 
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The next input/output menu functions which need ta be 
implemented are the saving to disk and loading from disk 
of control vertex positions, allowing work in progress to 
be saved and restarted later at the point from which it 
was saved. This task will be very similar to that of 
implementing the print function, except that some effort 
may be required to avoid having control vertex 
manipulations for the currently displayed model apply to a 
new set of vertices when loaded from a disk file. 

The most difficult input/output menu option to 
implement may be the plot function, intended as a screen 
dump. While the host program is aware of current vertex 
locations relative to one another, it is not informed of 
rotations, scaling, global translations, and clipping 
which are performed by the PS 300. Even if this 
information were communicated to the FORTRAN program, the 
programming required to accomplish the matrix 
transformations in the host-resident program would be 
sobering. The PS 300 programming language has a function 
called XFDRNDATA, however, which retrieves the transformed 
data for an object as ASCII data for transmission to the 
host computer. Using this function, it should be possible 
to implement a screen dump function without excessive 
difficulty . 

S.5.5. Inversion algorithm. The time required to 
model a desired surface could be greatly reduce with the 
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implementation of an inversion algorithm. Such an 
algorithm would calculate an initial set of control vertex 
positions, the required surface order, the homogeneous 
coordinates and the knot vectors which would yield a B- 
spline surface approximating the surface represented by 
data in the known-surface data file. The program would 
branch to such a subroutine almost immediately after 
program execution, from subroutine INITCP. 

Because the B-spline surface does not interpolate the 
control vertices, developing and implementing a suitable 
inversion procedure is not a trivial task. Two- 
dimensional approaches to the problem are presented by 
Yamaguchi C55J and Wu, Abel, and Breenberg C211. A more 
complex Cbut more immediately applicable) three- 
dimensional algorithm is presented by Barsky and 
Breenberg 1112,561. Rogers, Satterfield, and Rodriguez 
offer yet another three-dimensional approach C5,171, which 
they describe as ’’conceptually simpler, but 
computationally less efficient” than Barsky and 
Breenberg ’s. They offer as Justification a reminder that 
the inversion algorithm is used only once during the 
program, to obtain an initial approximation. Tiller C18J 
presents some thoughts relevant to approximating surfaces 
with nonuniform rational B-splines. 

It should not be assumed that an inversion algorithm 
will render interactive modeling obsolete. The algorithm 
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offered by Rogers, Satterfield, and Rodriguez, for 
example, does not automatically create hard chines, which 
must then be added interactively . Even if an inversion 
algorithm was capable of yielding a B-spline surface which 
exactly interpolated the known surface data, some 
interactive manipulation would be required to make the 
surface suitable. Digitized data from a sculpted surface 
is almost never fair, having small oscillations due to the 
digitization process. These oscillations must be removed 
either by extensive and difficult preprocessing of the 
data or through interactive manipulation. 

5.5.3. Circular arcs. Consideration of Tiller’s work 
with the representation of circular arcs by nonuniform 
rational B-splines suggests an additional feature useful 
during interactive modeling. With a suitable inversion 
algorithm, the user could place several control vertices 
under program control, identify three points on a desired 
circular arc, and allow the program to relocate the 
vertices to yield a circular arc interpolating the three 
points. Since the manipulation of homogeneous coordinates 
and of nonuniform knot spacing is not highly intuitive, 
such an approach is necessary in order to make use of the 
power of the nonuniform rational B-spline. 

5.3.4. Orthogonal Plane intersection mappino. The 
parametric representation of the surface presented by 
SPLINE is adequate for accurate visualization of the B- 
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spline surface. For some applications, however , a more 
traditional representation of waterlines, buttocks lines, 
and transverse sections is desirable. In order to 
implement such a display, it is necessary to map the 
intersection of the B-spline surface with the three basic 
orthogonal planes. 

The PS 300-resident portion of SPLINE, SPLINE.DAT, 
already has the necessary structure to display the three 
intersection maps. SPLINE. FTN would have to present the 
vector lists to the PS 300 with a command of the format: 
XY: -UECT0R_LIST ITEtllZED N- . . . . 

The vector list would be followed by the sending of an 
appropriate value to the PS 300 program variable, XYMQDE, 
and the sending of a triggering message to input number 
one of NEW_M0DE . The communication of control graph 
vector lists in subroutine CPLINE provides a similar 
example . 

Calculation of the vector list should be implemented 
in subroutine PLANES, which is currently called from 
subroutine WAIT when the appropriate function key is 
pressed . 

Mapping of the intersection between the planes and 
the B-spline surface can be difficult, especially when 
trying to obtain sufficient accuracy for an adequate 
display. Multiple disconnected contours, such as contours 
of adjacent peaks, present a particular problem. 
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Heap CE71 presents three algorithms For the 
production of contour maps of a Function deFined at the 
vertices oF an irregular triangular mesh. Such a mesh is 
oFten used in Finite element techniques, particularly For 
the solution oF partial diFFerential equations. Heap’s 
algorithm First creates the triangular mesh, with vertices 
oF the mesh lying on the surface to be mapped. The 
algorithm then Follows the contour From element to element 
oF the mesh. Since straight lines are generated through 
each mesh element, a very Fine mesh is required For an 
accurate contour, consuming large amounts oF memory. 
SPLINE already approaches the memory limit oF the PDP- 
11 . 

SatterField and Rogers C28J suggest using Heap’s 
algorithm to obtain a triangular mesh along the path oF 
the contour, but then present a variation. They use the 
triangular mesh as a basis For piecewise B-spline 
calculation oF the surFace in the region oF the mesh. 
Pairs oF surFace points separated by a desired tolerance 
are generated on the B-spline surFace and checked For 
spanning oF the contour. When they are Found to span the 
contour, a more exact contour location is determined by 
interpolating between the two spanning paints. This 
procedure handles properly the problems oF correctly 
ordering the points oF a contour and oF multiple 
disconnected contours in a single contour plane. 
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6.5.5. Reflection about a plane. Like the orthogonal 
plane intersection mapping, display of a reflection of the 
model about a plane has been fully implemented in 
SPLINE.DAT, but remains to be implemented in SPLINE. FTN. 
Calculation of the vector list would be trivial, requiring 
only that an offset be added to or deducted from the 
appropriate vector information already held in memory . 
Communication of the vector information would be similar 
to that detailed in section 6.5.4, above. Pressing the 
appropriate function key already activates a call to 
subroutine REFLECT. 

6.5.6. Calculation of form parameters. The 
mathematical nature of the surface representation could be 
made even more useful with automatic calculation of form 
parameters, such as center of gravity, section areas, or 
waterplane area. No provision has been made in SPLINE at 
this time for such a feature. Cruetz and Schubert C593 
take this approach one step further yet. They propose 
allowing the designer to specify initial form parameters, 
such as a section area curve, block coefficient, and beam 
to draft ratio, and letting the computer program generate 
from that an initial body plan. 

6.3. Conclusions 

SPLINE has demonstrated the tremendous potential 
which exists for mathematical definition of sculptured 



surfaces 



coupled with highly interactive computer 
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graphics. Dramatic changes in the efficiency with which 
automobile bodies, aircraft fuselages, and ship hulls are 
designed and manufactured can be expected in the very near 
future. Using B-splines, or similar mathematical 
representations, the definition of a complex sculptured 
surface will take hours, not days. Modifications to a 
developed surface will be accomplished in mere minutes, no 
longer requiring the designer to essentially start over 
with calculations and fairing of lines. 

Not only will surface definitions be achieved more 
accurately and efficiently, but those definitions will be 
directly interfaced with computer programs which will 
evaluate the performance characteristics of the form and 
with the machinery which will manufacture it. Further 
effort in the development of accurate inversion 
algorithms, in user/computer interfacing, and in 
generation of form parameters from surface definitions 
will be particularly important. 
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APPENDIX A 



FLOW CHARTS AND FUNCTION NETWORKS 

The Figures contained in this appendix present flow 
charts and Function networks For the major program 
Functions oF both SPLINE. FTN Cthe host-resident program} 
and SPLINE.DAT Cthe PS 300-resident program}. The program 
listings are given in appendices B and C. The comments 
contained in those listings will be particularly useFul in 
understanding the Function networks. Those who are 
unFamiliar with the PS 300 programming language will want 
to consult the computer manual C203 For an explanation oF 
display trees and Function networking and For details 
concerning the individual Function blocks. 
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Fig. 45. SPLINE. FTN Flaw chart. 
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KL09 



Fig. 45. SPLINE.DAT display tree. 
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ROTAT I ON 



FETCHMODE 




Fig. 47 



Rotation function network. 
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HOST 

nODE 



SCALING 



FETCHMODE 



SW ] T CH_SCALE 




Fig. 4B . Scaling function network. 
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Fig. 43. Typical function network for function key 
mode selection and labeling. 
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DIAL1 

] NPUTS_ 
CHOOSE! 5 ) 



’GLOBAL X'- 
’ GLOBAL X’- 
’ VERTEX X'- 



FETCHMODE 






HOST 

MODE 




< J 

< 2 

< 3 

< 4 

< 5 



\7 



-< 4 ) 



< ] > 




FALSE- 
FALSE- 
TRUE - 
FALSE- 



BL ] NKER 
] NPUTS_ 
CHOOSE! 5 ) 



< 1 ) 
< 2 ) 

< 3 ) 

< 4 ) 

< 5 ) 



< 1 ) 






Fxg. 50. Function network for blinking of 
translation dials when in interactive modeling mode. 
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FETCHCL I PMODE 



HOST 
CL I PMODE 



F : FETCH 
< 1 > < 1 > 
< 2 >C 



CL IP- TRUE 
] NPUTS- 
CHOOSE ( 3 ) 
FALSE - < J )C <1 
TRUE — < 2 >C 

-ft <JL±_ 




4 



CLIP-DIAL 
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< 1 >C ( 1 ) 

-fr |< 2 > < 2 ) H 
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DISCARD 
VARIABLE 



CLIP-TOTAL 



0 
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< J ) < I 
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— m ? 




Fig. 51. Clipping function network. 
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DISCARD 
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Fig. 5 2. Picking function network . 
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CONVERT PICK ID 
TO INTEGER 



P] CK_ASC ] 1 




Fig. 53. PICKID conversion Function network . 
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Fig. 54. Function network for global translation and 
translation of selected vertex during interactive 
modeling . 
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FETCHMODE 



SV ] TCH_TRAN 




Fig. 55. Function network for reporting of 
movement to FORTRAN program . 



vertex 
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FIG. 

53 
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SET-P I CKED2 
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< 1 > ( 1 > 
< 2 >C 
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SV ]TCH_P]CK 
F CROUTEl 90 ) 




INITIALIZED 
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Fig. SB. Function network to reset vertex movement 
accumulator when new vertex is picked. 



APPENDIX B 



SPLINE. FTN PROGRAM LISTING 



SPLINE. FTN is the host-resident portion of 
program, SPLINE. 



the 
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mm-m mm-m mm m mm * ***** mm mmmmmmmmmmmmmmm mmm mm mm mmm ***** *•* 

mmmmmmmmmmmmmmmmmmm program SPLINE. FTN ■*-*mmmmm-*mmmmmmmmmmmm 
mmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmm 



mmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmrn 



m 

♦Input 

♦Output 

♦Calls 

♦Alters 

♦Description 

m 

m 

m 

* 



m 



m 

♦Author : : 

♦Created On : 



* 



m 

INITPS, PSEXIT, INITCP,WAIT * 

* 

This program facilitates three-dimensional * 
surface modeling or design using a surface * 
B-spline algorithm. The surface is displayed* 
and interactively manipulated and modified * 
in real time on the PS-300 vector refresh * 
display. The file, ’’SPLINE.DAT,” contains * 
the PS-300 command language program which * 
produces and controls the graphics display * 
using data supplied by SPLINE. FTN. * 

□. Michael Bryant * 

5/21/B7 Modified On: * 

* 



***f*****itii*****iii«»*****«.(.(««tt«*»*..««**.«**«*«*ii*ii«..»* 



♦Uariable 

■at 


Type 


Explanation 


•* 


* ISTAT C 1 3 


Integer 


Reports error status from PS-300 


m 



*************************************************«********* 



PROGRAM SPLINE 

COMMON /SPLINE/ MOP, NCP, CPC10,S,43, MORDER, NORDER 
COMMON /BASIS/ RC13,4,313, SC11,4,313, SURF C 3 , 31 , 31 3 
COMMON /UECS/ UECSC 10003 
DIMENSION ISTATC13 

***** Send PS-300 command language program to PS-300 and 
***** execute. 

CALL INITPS 

***** Establish initial control point locations. 

CALL INI TCP 

***** p^ait and act upon signals from PS-300 program. 

CALL WAIT 

***** Oe-establish link between PS-300 and PDP-11, using 
***** library subroutine. 

ISTATC13-0 

CALL PSEX I T C 1 , I STAT 3 

STOP 

END 
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••******•****••****••**•*****»**•****••***•*••••**••**•*•** 
******** Subroutines Follow in Alphabetical Order «****««»* 

mmmmmmmmmmmmmmmmmmm mmmmmmmmmmmmmmmmmmm mmmmmmmmmmmmmmmmmmmrn m 



mmmmmmmmmmmmmmmmmmmm Subroutine BLANK ********************* 



♦Input 

♦Output 

♦Calls 

♦Alters 

♦Description 



* 



* 



* 

^Author : 

^Created On : 



I, BUFFER * 

BUFFER * 

* 

BUFFER * 

This program makes extensive use of the * 

’ENCODE’ command to convert integers to * 

characters in a LQGICAL*1 array C ’BUFFER’D. * 
The buffer is then passes to the PS-300, * 

with the encoded integers forming portions * 
of the names of nodes, I.E. ’CP23’. When a * 
single-digit integer is encoded with an ’12’* 
specification, however, the result is not * 
’03’, for example, but ’ 3’. The blank * 

would not be allowed as part of a node name.* 
This subroutine replaces the unwanted blanks* 
with ”0’s”. * 

D. Michael Bryant * 

11/17/86 Modified On: * 



* * 

*Uariable Type Explanation * 

4t *t 

*BUFFERC80) Logical Contains the string to be sent to * 

* the PS-300. * 

*1 Integer Identifies the array element which * 

* contains the blank to be replaced. * 

* * 

*********************************************************** 

SUBROUTINE BLANKC I , BUFFER) 

L0GICAL*1 BUFFERC80) 

IF CBUFFERCID .ED. ’ ’ D BUFFERC I D- ’ 0 ’ 



RETURN 

END 
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* ***************************************************** ***** 
*********************************************************** 



******************** Subroutine BSPLIN ******************** 
*********************************************************** 



*********************************************************** 



* 



* 



* Input 

^Output 

* 

♦Calls : 

♦Alters : 

♦Description : 

* 

* 

* 

♦Author: 
♦Created On : 

* 



MCP , NCP , CPC 10,9,11, MORDER , NORDER ♦ 
R(13, 4, 313, 5(11,4,313, UKNOT C 153 , UJKNOT C 131 , * 
UC313 , UJC313 ♦ 
PSSEND, SURFAC * 

* 

This subroutine calculates a full B-spline ♦ 
surface. Subroutine SURFAC is then called ♦ 
to extract and send a vector list to the ♦ 
PS300 for display of the calculated surface.* 
D. Michael Bryant * 
2/22/07 Modified On: ♦ 

* 



*********************************************************** 



* * 

*Uariable Type Explanation * 

m — ___________ * 



♦A 


Real 


* 




*B 


Real 


*C 


Real 


* 





*CPC 10,8,4) 

* 

* 

* 

* 


Real 


*D 


Real 


*ncp 

* 


Integer 


*M0RDER 

* 


Integer 


*NCP 

* 


Integer 


*N0RDER 

* 


Integer 


*RC 13 , 4 , 31 ) 

* 

* 

* 

* 

* 

* 


Real 


*S( 11,4,31) 

* 

* 

* 


Real 



Denominator of the first term of the ♦ 
basis function. ♦ 

First term of the basis function. ♦ 

Denominator of the second term of the* 
basis function. * 

Holds the X,Y,Z coordinates for each * 
of the control points (vertices! . ♦ 

Subscripts denote row, column of * 

vertex and X,Y,Z or homogeneous * 

component . * 

Second term of the basis function. * 

Number of rows of vertices in the * 

B-spline control polygon. * 

The order of the B-spline surface in * 
the longitudinal Cu3 direction. * 

Number of columns of vertices in the * 
B-spline control polygon. 

The order of the B-spline surface in 
the transverse Cuj 3 direction. 

The basis functions in the longitud- 
inal direction. The first subscript 
indicates the associated vertex row, 
the second indicates the order, and 
the third indicates the number of 



parametric increments along the curve* 
(corresponding to the ”u” value! . * 

The basis functions in the tranverse* 
direction. The first subscript indi-* 
cates the associated vertex column, * 
the second indicates the order, and * 
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m 

m 

0 

*SURFC3, 31,313 


the third indicates the number of * 
parametric increments along the curve 1 * 
([corresponding to the ” uj ” value). * 
The calculated points on the B-spline* 


m 


Real 


surface, in cartesian coordinates, as* 


0 

m 

*UC31 ) 


Real 


a function of the parametric * 
coordinates. * 
The parametric value in the longitud-* 


m 

0 

*LIINC 


Real 


inal Cu) direction at each of 31 * 
increments along the curve. * 
Parametric spacing of the 31 ”u” * 


m 

♦UKNDTC1S5 


Real 


lines. * 

The knot vector in the longitudinal * 


m 

*umx 


Real 


Cu) direction. * 

The maximum value in the ”u” knot * 


* 

*ldC 31 ) 


Real 


vector . * 

The parametric value in the trans- * 


* 

m 

*UINC 


Real 


verse Cuj) direction at each of 31 * 
increments along the curve. * 
Parametric spacing of the 31 ” uj ” * 


m 

*UIKNDT C 13 D 


Real 


lines. * 

The knot vector in the transverse Cuj)* 


m 

*umx 


Real 


direction. * 

The maximum value in the ” uj ” knot * 


* vector . * 

*0 00000*0000000000*0 + 0 0000000000000 0 00* 0*0000 00* *0000000 0**0 



SUBROUTINE BSPLIN 

connoN /spline/ ncp, ncp, cpcio,b,4), morder, norder 
COMMON /BASIS/ RC13,4,31), sell, 4, 31), SURFC3.31.31) 
DIMENSION UKN0TC15), WKN0TC13), UC31), UJC31) 

***** Blink ’’B-SPLINE” label on keyboard while calculating. 
CALL PSSEND C 24, ’SEND TRUE TO <2>FLABEL2 ; ’ ) 

CALL PSSEND C22,’SEND 1 TO < 1>NEUJ_M0DE; ’ ) 

***** Determine knot vectors Cdefault * uniform). 

UMAX - REALCMCP-MORDER+1) 

UMAX - REALCNCP-NORDER+1) 

DO 150 K-l .MCP+MORDER 

IF CK.LE. MORDER) UKNOTCJO-O.O 
IF C K . GT . MORDER . AND . K . LE . MCP ) UKNOTCK) - 
+ UKNOTCK-D + l .0 

IF CK . GT .MCP) UKNOTCK) - UMAX 
150 CONTINUE 
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DO ISO K - 1 , NCP+NORDER 

IF CK . LE . NORDER3 UIKNOTCK3 - 0.0 
IF CK . GT . NORDER . AND . K . LE . NCP 3 UJKN0TCK3 - 
+ WKNOTCK-13+1 .0 

IF CK.GT.NCP3 UKN0TCK3 * UIMAX 
160 CONTINUE 

***** Determine parametric spacing cf displayed B-spline 
***** surface net lines. 

UINC - UHAX/30 . 0 
UINC - UiriAX/30 .0 

DO 200 K - 1,31 

UCK3 - CREALCK-133*UINC 
UJ C K 0 - CREALCK-13 3*UJINC 
200 CONTINUE 

***** Determine first order basis functions. 

DO 250 M - l.nCP+nORDER-l 
DO 240 K - 1,31 

IFCUKNOTCm .LE.UCK3 . AND . UCK 3 . LT . UKNOT C f1+ 1 3 3 
+ THEN 

RCn, 1,10 -1.0 

ELSE IF CUKN0TCn3.EQ.UCKn THEN 
RCfl, 1,K3 -1.0 

ELSEIF CK.EQ.31 . AND . M . EQ . f1CP3 THEN 
RCH, 1,K3 -1.0 

ELSE 

RCN, 1 ,K3 -0.0 

END IF 

240 CONTINUE 

250 CONTINUE 

DO 270 N - 1 , NCP+NQRDER-1 
DO 260 K - 1,31 

IF CUJKN0TCN3 .LE.UJCK3 .AND.UJCK3 .LT . WKNOTCN+1 3 3 
+ THEN 

SCN, 1 ,K3 - 1.0 

ELSEIF CWKN0TCN3 .EQ.WCK33 THEN 
SCN, 1 , K 3 - 1.0 

ELSEIF CK.EQ.31 .AND. N.EQ.NCP3 THEN 
SCN, 1,K3 -1.0 
ELSE 

SCN, 1,K3 - 0.0 

END IF 

260 CONTINUE 

270 CONTINUE 
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290 

330 

340 

350 



390 



390 

430 

440 

450 



Calculate higher order basis functions. 

DO 350 K - 2 , flORDER 

DO 340 M - l.MCP+MORDER-K 
DO 330 I - 1,31 

A - UKNOTCtl+K-1 3 - UKN0TCN3 
IF CA. ED. 0.03 THEN 
B - 0.0 
EL5E 

B - CUC n-UKNOTCm 3*RCf1, K-l , I3/A 

ENDIF 

C - UKNOTCfl+lO - UKNOTCM+13 
IF CC.EQ.0.03 THEN 
D " 0.0 
ELSE 

D-CUKN0TCM+K3-UC I 3 3*RCh+l , K-l , I 3/C 
ENDIF 

RCM.K, 13 - B + D 
CONTINUE 
CONTINUE 
CONTINUE 

DO 450 K - 2 , NORDER 

DO 440 N - 1 , NCP+NORDER-K 
DO 430 I - 1,31 

A - UKNOT C N+K-l 3 - UKN0TCN3 
IF CA.EQ.0.03 THEN 
B - 0.0 
ELSE 

B - CUIC I 3-UKNOTCN3 3*SCN, K-l, 13/A 

ENDIF 

C - UKNOTCN+K 3 - UKNOTCN+13 
IF CC.EQ.0.03 THEN 
D “ 0.0 
ELSE 

D- C UKNOTCN+K 3 -UC I 3 3+SCN+l , K-l , 13/C 
ENDIF 

SCN, K, 13 - B + D 
CONTINUE 
CONTINUE 
CONTINUE 
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***** Cal CU late the surface coordinates . 

DO BOO I - l.MCP 

DO 730 J - 1 , NCP 

DO 700 L - 1,31 

IF CRC I , MORDER, L3 . EQ .0 . 03 GOTO 7B0 
DO 770 M - 1,31 

IFCSC J , NQRDER , M3 .EQ.0.Q3 GOTO 770 
DO 7B0 K - 1,3 

SURFCK , L , M3 - SURFCK,L,M3 + 

+ CPC I , J , K3 * 

+ CPC I , J , 43 * 

+ RC I .MORDER, L3* 

+ SC J. NORDER, M3 

CONTINUE 
CONTINUE 
CONTINUE 
CONTINUE 
CONTINUE 

***** Send the surface coordinates to the PS-300. 

CALL SURFAC 

***** Stop blinking ’’B-SPLINE” label on keyboard. 

CALL PSSEND C 25,’ SEND FALSE TO <2>FLABEL2 ; ’ 3 
CALL PSSEND C 22,’ SEND 1 TO < 1>NEU1_M0DE ; ' 3 

RETURN 

END 



7B0 

770 

7B0 

730 

BOO 
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*********************************************************** 

******************** Subroutine CPLINE *************** : * : **** 
*********************************************************** 

*********************************************************** 



* 

♦Input : 

♦Output : 

* 

♦Calls : 

♦Alters : 

♦Description : 
* 



* 

♦Created On : 
* 



* 

mcp.ncp.cpcio.b, 4) ♦ 

UECSC500) , LPC333) , C vector list defining ♦ 

lines connecting vertices) ♦ 

PSSEND , PSUECS * 

SPLINE . DAT structure ”CP_LINES” * 

This subroutine prepares in the UECS array a* 
vector list defining the lines connecting * 

the B-spline vertices Ccontrol points) . The* 
vector list is then passed to the PS-300. * 

13/0S/S6 Modified On: * 



*********************************************************** 

* 
* 



♦Uariable 



Type 



Explanation 



♦BUFFERC80) Logical 

* 

* 

*CPC10,e,4) Real 



* 

* 

♦I 



Integer 



♦ILP 

m 

* IOSTAT 

♦IUEC 

* 

♦LPC333) 

* 

* 

* 



Integer 

Integer 

Integer 

Logical 



* 

m 

*M 



♦MCP 

* 

♦N 

♦NCP 



Integer 

Integer 

Integer 

Integer 



♦UECSC500) Real 



Holds a message for transmital to the 
PS-300 program with the PSSEND 
library subroutine. 

Holds the X,Y,Z coordinates for each 
of the control points Cvertices) . 
Subscripts denote row, column of 
vertex and X,Y,Z or homogeneous 
component . 

Indicates the X,Y, or Z component of 
the vertex vector is under 
consideration . 

Counter indicates the current array 
element of LP . 

Input/output error status. 

Counter indicates the current array 
element of UECS. 

.TRUE, indicates corresponding 
elements of UECS represent a line 
vector (draw to this point with pen 
down). .FALSE, indicates 
corresponding elements of UECS 
represent a point vector Cmove to 
this point with pen up) . 

Row of vertices under consideration. 
Number of rows of vertices in the 
B-spline control polygon. 

Column of vertices. 

Number of columns of vertices in the 
B-spline control polygon. 

Uector list defining the lines 
connecting the B-spline vertices. 



* 

* 

* 



* 



* 

* 

* 

* 

* 

* 

* 



* 



* 

* 

* 

* 

* 

* 

* 



*********************************************************** 



135 



**»»* 



* 

* 



80 



90 

100 



180 



190 

200 



SUBROUTINE CPLINE 

connoN /uecs/ uecsc ioood 

LOGICAL LPC333D 

CONNON /SPLINE/ flCP , NCP , CP C 10 , 8 , 4 5 , flQROER , NORDER 
L0GICAL*1 BUFFERC80D 

IUEC-0 

ILP-0 

Uector list far lines along raws of vertices. 

do ioo n-i,ncp 

DO 90 N-l.NCP 

Count the vectors created. 

ILP-ILP+1 
DO 80 1-1,3 

Count the vector components CX,Y,ZD 
created . 

IUEC - IUEC + 1 
UECSC IUECD - CPCn,N,ID 
CONTINUE 



mmmm indicate whether vector is point or line vector. 

LPCILPD - .TRUE. 

IF CN.EQ.1D LPCILPD - .FALSE. 

CONTINUE 

CONTINUE 

Uector list for lines along columns of vertices. 

DO 200 N-l.NCP 

do 190 n-i,ncp 

ILP-ILP+1 
DO 180 1-1,3 

IUEC - IUEC + 1 
UECSC IUECD - CPCfl.N, ID 
CONTINUE 

**»« Indicate whether vector is point or line vector. 

LPCILPD - .TRUE. 

IF CN.EQ.1D LPCILPD - .FALSE. 

CONTINUE 

CONTINUE 
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***** 



210 



Send vector list to PS-300 program to define 
substructure CPL of structure CP_LINES. 

ENCODE C31 , 210, BUFFER!) ILP 

FORMAT C ’CPL: -UECTOR_LIST ITEMIZED N-’,I3) 
CALL PSSEND C 31, BUFFER 3 

CALL PSUECS C4, ILP.UECS.LP, .TRUE. , IOSTAT) 

RETURN 

END 
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*********************************************************** 

*********************************************************** 

******************** Subroutine INITCP ******************** 
*********************************************************** 
*********************************************************** 
* * 



♦Input 

♦Output 

♦Calls 

♦Alters 

♦Description 

* 



♦Created On : 

* 



RATIO, MCP, NCP 
CPC10, 0,41 

PSSEND , INUERT , BLANK , PREPS 
SPLINE . DAT 

This subroutine determines an initial set 
of control points Cverticesl and passes 
commands to the PS-300 which create the 
associated display tree branches. 

11/14/86 Modified On: 



m 

m 

m 

« 

* 

* 

* 

* 



* 



************************************************************ 



♦Uariable Type 



♦ANSWER C 1 ) Char 
♦CPC 10,8,43 Real 

* 

* 

* 

* 



♦DM 

* 


Real 


♦DN 

* 


Real 


♦ncp 

* 

* 


Integer 


♦flQRDER 

* 


Integer 


♦NCP 

* 

* 


Integer 


♦NQRDER 

* 


Integer 


♦Nuncp 

* 

* 

* 


Integer 


♦RATIO 

* 


Real 


♦YNIN 

* 

* 


Real 



Explanation 



Retrieves CYles and CNlo answers. 
Holds the X,Y,Z coordinates for each 
of the control points Cverticesl. 
Subscripts denote row, column of 
vertex and X,Y,Z or homogeneous 
component . 

The X increment between initial 
control point positions. 

The Y increment between initial 
control point positions. 

The number of control points in the 
X direction in the initial control 
point matrix. 

The order of the B-spline surface in 
the longitudinal Cul direction. 

The number of control points in the 
Y direction in the initial control 
point matrix. 

The order of the B-spline surface in 
the transverse Cwl direction. 

The identifying number of the control 
point under consideration, 1 — 80, 
where point Cl, 11 is 1, Cl, El is E, 

. . . , CIO, 81 is 80 . 

The length to width ratio of the 
initial control point matrix. 

The minimum Y coordinate of the 
initial control point matrix. 



*********************************************************** 
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SUBROUTINE INITCP 

connoN /spline/ ncp, ncp, cpc 10 , 8 , 4 ), morder, norder 

CHARACTER* 1 ANSWER 
LOGICAL*l BUFFERC80) 

***** Determine whether it is desired to display a known 
***** sur face . 

50 TYPE SO 

60 FORMATC///, ’ Do you wish to display a known ’ 

+ ’surface from a PREPS-Format file?’,/, 

+ ’ Please answer CY) or CN).’,/) 

READC5 , 70 , ERR”50) ANSWER 
70 FORMAT CA1) 

IF C ANSWER . EQ . ’ Y ’ 5 SOTO 100 
IF C ANSWER. EO. ’N’ ) SOTO 200 
SOTO 50 

***** fliiou input of coordinates of points on a known 
***** surface from a PREPS-Format file. 

100 CALL PREPS 

***** Establish size and length/width ratio of initial 
***** control point matrix. 

200 TYPE 201 

201 FORMATC//’ A matrix of vertices up to CIO X 81 ’ , 

+ ’ may be selected.’,/, 

+ ’ ENTER THE DESIRED MATRIX SIZE.’,/, 

+ ’ (Number of transverse, longitudinal ’, 

+ ’ stations . 1 ’ , / , 

+ ’ Example: 7,5’,/) 

READ C5,*,ERR-200) MCP,NCP 

IF CMCP.GT.10.0R.NCP.GT.8) THEN 
TYPE 220 

FORMAT C/,’ The maximum value for the’, 
+ ’ matrix is 10,8.’/, 

+ ’ ENTER BOTH UALUES AGAIN.’) 

GOTO 210 

ELSEIF CMCP.LT. 5. OR.NCP.LT. 5) THEN 
TYPE 230 

FORMAT C/, ’ The minimum value for the’, 
+ ’ matrix is 5,5.’/, 

+ ’ ENTER BOTH UALUES AGAIN.’) 

GOTO 210 



210 

220 

230 



END IF 
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235 TYPE 240 

240 FORNATC///, ’ Enter the desired length/width 

+ ’ ratio for the vertices matrix.’,/! 

READ C5, *, ERR-2351 RATIO 
IF C RATIO .LE .0! THEN 
TYPE 250 

250 FORNAT C / , ’ THE LENGTH/lilIDTH RATIO NUST ’ , 

+ ’ BE POSIT IUE . ’ 1 

GOTO 235 

END IF 

***** Determine the order of the B-spline surface equation. 
TYPE 255 

255 FORflAT C//, ’ Enter the order of the B-spline’,/, 

+ ’ surface in the longitudinal direction.’! 

270 TYPE 275 

275 FORNAT C’ Use integer values of 3 or 4. ’! 

READ C5, *, ERR-270! NORDER 

IF CnORDER.LT. 3. OR. NORDER. GT. 4! GOTO 270 

TYPE 280 

280 FORNAT C//, ’ Enter the order of the B-spline’,/, 

+ ’ surface in the transverse direction.’! 

285 TYPE 275 

READ C5, *, ERR-285! NORDER 

IF CNORDER . LT . 3 . OR .NORDER .GT .4! GOTO 285 

***** Calculate the coordinates of the initial vertices 
***** matrix. 

»«*«•*»*** Determine the X and Y spacing of the vertices. 

DN“1 . 0/ FLOAT CNCP-1 ! 

IF CNCP.EQ.il THEN 
DN-0.0 

ELSE 

DN=*1 . 0/ C RAT IQ*FL0AT C NCP-1 1 1 
END IF 

*»****«*«» Calculate the coordinates. 

YNIN - 0.0-1 .0/CRATI0*2. 01 

DO 310 N-l.NCP 

DO 300 N-l.NCP 

CPCN , N , 1 1 - -0.5 + FLOAT CN-11 * DN 
CPCN, N, 21 - YNIN + FLOATCN-11 * DN 
CPCN, N, 31 - 0.0 
CPCN, N, 41 - 1.0 

* Default homogeneous coordinate - 1 for 

* non-rational B-spline. 



300 CONTINUE 

310 CONTINUE 
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***** Send the necessary program lines to the PS-300 to 
***** display the vertices correctly positioned. 

500 CONTINUE 



DO 700 M - 1,MCP 

DO 630 N - 1 , NCP 

NUMCP - M + CN-1 )*10 

ENCODE C50 , 510 , BUFFER} NUMCP , NUMCP 
510 FORMAT C ’CP ’, 12, ’: -BEGIN .STRUCTURE SET’, 

+ ’ PICKING IDENTIFIER-CP’ , 12, ’ ; ’ ) 

CALL BLANK C 3, BUFFER) 

CALL BLANKC48, BUFFER) 

CALL PSSENDC50, BUFFER) 

CALL PSSEND C 25 ,’ TRAN : -TRANSLATE BY 0,0,0;’) 
CALL PSSEND C 17, ’ INSTANCE OF CUBE;’) 

CALL PSSEND C 14 , ’ END ..STRUCTURE ; ’ ) 

ENCODE C28, 520, BUFFER) NUMCP 
520 FORMAT C ’ TRAN_TOTAL_’ , 12 , ’ : -F : ’ , 

+ ’ACCUMULATE;’) 

CALL BLANK C 12, BUFFER) 

CALL PSSEND C2B, BUFFER) 

ENCODE C3B, 530, BUFFER) NUMCP , NUMCP 
530 FORMAT C ’CONNECT TRAN_TOTAL_’ , 12, 

+ ’<1> :<1>CP’ , 12, ’ .TRAN; ’ ) 

CALL BLANKC20, BUFFER) 

CALL BLANKC31 , BUFFER) 

CALL PSSEND C3B, BUFFER) 

ENCODE C52, 540, BUFFER) CPCM.N.l), CPCM,N,2), 
+ CP CM , N , 3 ) , NUMCP 

540 FORMAT C ’ SEND UC ’ , F7 . 4 , ’ , ’ , F7 . 4 , ’ , ’ , 

+ F7 . 4 , ’ ) TO < 2> TRAN_TOTAL_ ’ , 12, ’ ; ’ ) 

CALL BLANKC50, BUFFER) 

CALL PSSEND C52, BUFFER) 

ENCODE C34, 545, BUFFER) NUMCP 
545 FORMAT C ’ SEND UC0,0,0) TO’, 

+ ’ <1>TRAN_T0TAL„’ , 12, ’ ; ’ ) 

CALL BLANKC32, BUFFER) 

CALL PSSENDC34, BUFFER) 
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ENCODE CBB, 550, BUFFER) NUMCP 

550 FORMAT C ’SEND .1 TO <4>TRAN_T0TAL_’ , IB, 

* ’ 5 ’5 

CALL BLANK CBS, BUFFER) 

CALL PSSEND CBB, BUFFER) 

ENCODE C 41 , 5S0, BUFFER) NUMCP , NUMCP 
560 FORMAT C ’ CONNECT SUITCH_PICK< ’ , IB, 

+ ’ > : <1>TRAN_T0TAL_’ , IB, ’ ; ’ ) 

CALL BLANK CB1, BUFFER) 

CALL BLANK C 33, BUFFER) 

CALL PSSEND C 41, BUFFER) 



***** Connect the dial signal to the switch outputs for row 
***** anc j C oi um n movement. 



ENCODE C41, 560, BUFFER) M+BO, NUMCP 
CALL BLANK C39, BUFFER) 

CALL PSSEND C41, BUFFER) 



ENCODE C41 , 560, BUFFER) N+SO, NUMCP 
CALL BLANK C39, BUFFER) 

CALL PSSEND C41, BUFFER) 



690 CONTINUE 

700 CONTINUE 



TYPE 701 

701 FORMAT C / , ’ 

+ ’ 

- /, * 

+ /, ’ 

+ /, ’ 



************************************** > 

** PRESS < SH I FT-LI NE/LOCAL> ** ’ 

** 

** 



ON THE PS-300 KEYBOARD ** ’ 

FOLLOWED BY <TERM> . ** ’ 



************************************** 



//////) 



RETURN 

END 
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A*-*-****-**-*-**-*-*-**-******-***-*****-*-**-**-*-*****#*-*-*-***-*-***-*-*-*-**-*-* 
****************«*****«****•*•*•*-*•*-«*•*•*•*•***•*•**•***•*-*•*•*-*-•*•**•*-*•**•* 
*****-**•*•**•**•*•*-*•***-*•* Subroutine INITPS ***■**•*•*■*•*•*•*•*■*’*■*•*•*■*■** 

**-«-*-*-*-*-*-*-*•*-* m-mmmm + + mmm*-mmm-*-*m-m-mmmm m*km**mmm*-mmm**m*mm***k*-mm* 






m 



* 



“Input : PS-300 command language program file, 

* ’’SPLINE . DAT” 

“Output : PS-300 command language file output to 

* PS-300. 

“Calls : PSETUP, PSSEND, PSEXIT 

“Alters : 

“Description: This subroutine sends the PS-300 command 

* language program stored in the PDP-11 disk 

* file, ’’SPLINE.DAT, ” to the PS-300. 

“Created On : 11/3/06 Modified On: 

* 



* 

* 

* 



* 



* 

* 



*********************************************************** 



“Uariable 

m 


Type 


Explanation 




* 

m 


“BUFFERCBO) 


Logical 


Used to pass each program line to 


-* 


* 




PS-300. 




* 


“ I COUNT 


Integer 


Number of characters 


in program line 


* 


* 




currently being read 


and passed . 


* 


“ISTATC1D 

* 


Integer 


Reports error status 


from PS-300. 


m 

m 



*•»***•••«*****»•**•*•**•*«••***••**•«•*»*«•**«***•*•«*•**« 

SUBROUTINE INITPS 
DIMENSION ISTATCID 
LOGI CAL*1 BUFFERCBOD 
ISTATC1D-0 

*“““* Initialize PS-300 and set up link with PDP-11 using 
***** library subroutine. The first two parameters set up 
***** logical unit numbers for input from and output to the 
***** ps-300. The third parameter indicates whether of not 
***** the PS-300 is to issue an ”INIT” command. The last 
***** two parameters indicate the length and name of the 
***** array used to report error status from the PS-300. 

TYPE 100 

200 FORMAT C ’ a**-*-*-*-*-*-*-*-*-*-**-*-*-**-*-*-*-*-*'*-*-*'*-*-***-*'*'*'**- ’ ^ 

+ / i ’ ******** SETTING UP PS-300 ********* ’ , 

+ / » ************************************’ / } 
CALL PSETUPC7.7, .TRUE. ,1, ISTATD 
CALL PSSENDC11, ’ INITIALIZE; ’ 5 
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***** Send File to PS-300 one line at a time, using library 
***** subroutine. 

OPENCUNIT-3, NAME” ’SPLINE .DAT ’ , TYPE- ’OLD ’ , ERR-210) 

LINE - 0 

120 READC3, 130, END-140, ERR-220) I COUNT, CBUFFERC I ) , 

+ I -1,1 COUNT) 

130 FORMATCQ, 132A1) 

IFC I COUNT . GT . 0 )CALL PSSENDC I COUNT , BUFFER ) 

LINE - LINE+1 

MARK - LINE-LINE/ 50*50 

IF CMARK.EQ.O) TYPE 135 

135 FORMAT C’ Working ’ ) 

IF CLINE. EQ. 1000) TYPE 135 
135 FORMAT C ’ Don ’ ’ t give up on me . . . . ’ ) 

GOTO 120 

140 CL0SEC3) 

RETURN 

***** Error Routines. 

210 TYPE 211 

211 FORMAT C ’ ERROR WHILE OPENING DISK FILE ’, 

+ ’ SPLINE . DAT . ’ ) 

CALL PSEXIT Cl , I ST AT) 

GOTO 140 

220 TYPE 221 

221 FORMAT C ’ READ ERROR WHILE TRANSFERRING ’, 

+ ’DATA FILE. ’ ) 

CALL PSEXITC1 , ISTAT) 

GOTO 140 



END 

*W***M«W*W*******W***W«**«***«**««««*'»'*'»'*'l«*'*'*'*'»«'*'*'*'l»'**'»'*'»'l»'* 

*********************************-************************** 

********************** Subroutine 10 ********************** 
*********************************************************** 

SUBROUTINE IOCKEY) 

***** This subroutine directs all input/output functions, 
***** based on the function key press reported in the 
***** variable ’’KEY”. 



RETURN 

END 
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*********************************************************** 

*************************************-**********•*•****•**•*•**•** 

******************** Subroutine MQUING ******************** 
*********************************************************** 



*********************************************************** 
* * 



♦Input : 

* 

♦Output : 

♦Calls : 

♦Alters 

♦Description : 
* 

* 

* 

* 

* 

* 

* 

* 

* 

♦Created On : 
* 



I CP , IBS , flSG , flCP , NCP , CPC 10 , 8 , 4 ) , vector from * 
PS-300 program indicating vertex movement. * 
Updated CPC10,a,4) * 
PSSEND , PSREAD , CPL I NE , PSPOLL , BSPL I N * 
CPC10,a,41 * 
This subroutine receives vectors from the * 
PS-300 program which define vertex movement.* 
The array, CP, is updated to reflect that * 
movement. CPLINE is called to recalculate * 
the vector list for the lines connecting the* 
vertices if the lines are being displayed * 
and the B-spline surface is recalculated * 
using an abbreviated algorithm which takes * 
into account the localized influence of the * 
relocated vertex. * 
12/17/86 Modified On: * 



* 



*********************************************************** 



* 



* 



*Uariable Type Explanation 

* 



* 

* 



*BUFFER( 80 1 Logical 



* 

*cpcio, a, 

* 

* 

* 

*DX 

*DY 

*DZ 

* I BS 

* 

* 

*ICP 

* 

* 

* 

* I UEC 

* 

♦LENGTH 

* 

*n 

* 

*ncp 



) Real 

Real 
Real 
Real 
I nteger 

I nteger 

Integer 

Integer 

Integer 

Integer 



Used to pass text messages to and * 
from the PS-300 program. * 
The coordinates of the vertices, with* 
the subscripts indicating row, * 
column, and X,Y,2 or homogeneous * 
component . * 
The X movement of a selected vertex. * 
The Y movement of a selected vertex . * 
The Z movement of a selected vertex. * 
Indicates the display status of the * 
B-spline surface. 0 - not displayed, * 
l-4“varying combinations of spacing. * 
Indicates the display status of the * 
vertices. ICP-2 indicates the lines * 
connecting the vertices are * 
displayed. * 
The position of the picked vector in * 
the vector list. * 
Byte length of a message returned * 
through a call to PSPOLL or PSREAD. * 
The row of vertices under * 
consideration. * 
The number of rows of vertices. * 



145 



*nsg 

* 

* 



* 

* 

* 

* 

m 

m 

m 

m 

m 

m 

*N 

* 

*NCP 

*nun 

m 



m 

* 

* 



Integer 



Integer 

Integer 

Integer 



This variable is passed from the 
PS-300 program. A value of 1-BO 
indicates that the vector to be sent 
next represents movement of the 
correspondingly numbered vertex. 
Ualues of Bl-30 indicate the vector 
represents movement of a row of 
vertices, and values of 91-BB 
indicate movement of a column of 
vertices. Ualues greater than 100 
are used to report the picking of a 
roui or column of vertices, with the 
value being the position of the 
vector in the vector list, offset by 
100 . 

The column of vertices under 
consideration . 

The number of columns of vertices. 
Indicates which switch position the 
PS-300 program should send the 
translation dial signals to when a 
row or column of vertices has been 
picked for movement. 



* 

* 

* 

* 

m 

m 

m 

* 

* 

m 

m 

m 



m 

« 

m 



m 

m 



*•*•**»••»••*««*•••**••#•**•**********•*••**•••***«•*»*•••* 



SUBROUTINE NOUINGCICP, IBS , MSG , LENGTH , BUFFER) 

COfinON /SPLINE/ NCP , NCP , CPC 10 , B , 4 ) , flORDER , N0RDER 
L0GICAL*1 BUFFERCB0) 

***** If MSG>100, row or column of vertices has been picked 
***** for movement. Compute the appropriate identifying 
***** number C81-90 “ rows 1-10; Sl-SB - columns 1-B) and 
***** send it to the PS-300 program to use in routing the 
***** dial signals. 

IF Cf1SG.GT.100) THEN 
I UEC - MSG - 100 
IF C IUEC . LE . f1CP*NCP) THEN 

* a vector in a row of vertices was picked 

NUfl - C IUEC+NCP-15/NCP+B0 

ELSE 

* a vector in a column of vertices was picked 

NUN - C IUEC+NCP-1-NCP*NCP)/NCP+S0 

ENDIF 

ENCODE CBG, 100, BUFFER) NUN 

100 F0RNAT C ’SEND FIXC’,12,’) TO <1>PICKED;’) 

CALL PSSEND CB6, BUFFER) 

RETURN 



ENDIF 
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mmmmm 


A vertex has been moved . 


m 


IF CMSG.LT.100) THEN 

wait for the vector defining the movement 
CALL PSREAD CBO , BUFFER , LENGTH) 



********** Decode the message, which is in the form of 
»***»*«*«* »i_ x,Y,Z”, where X, Y, and Z are in an E10 
********** format if positive, and an Ell format if 

nSQ3tlVS 



200 


IF C LENGTH. EQ. 34) THEN 

DECODE C34, 200, BUFFER) DX,DY,DZ 
FORMAT C1X,3C1X,E10.4)) 


201 


ELSEIF C LENGTH . EQ . 37) THEN 

DECODE C37, 201 , BUFFER) DX , DY , DZ 
FORMAT C1X, 3C1X, Ell .4) ) 


202 


ELSEIF CLENGTH.EQ .35) THEN 

IF CBUFFERC3) . EQ . ’ - ’ ) THEN 

DECODE C35, 202, BUFFER) DX,DY,DZ 
FORMAT C2X.E11 .4, 2C IX, E10 . 4) ) 


203 


ELSEIF C BUFFER C 14 ) . EQ . ’ - ’ ) THEN 

DECODE C35, 203, BUFFER) DX , DY , DZ 
FORMAT C2X, E10.4, IX, Ell .4, 1X.E10.4) 


204 


ELSE 

DECODE C35, 204, BUFFER) DX , DY , DZ 
FORMAT C1X,2C1X,E10.4) , IX, Ell. 4) 

END IF 


-t 

205 


ELSEIF CLENGTH.EQ. 36) THEN 

IF CBUFFERC3) . EQ . . AND . BUFFER C 15 ) . EQ . ) 

THEN 

DECODE C 36, 205, BUFFER) DX , DY , DZ 
FORMAT C IX, 2 C IX, Ell .4) , IX, E10.4) 


206 


ELSEIF CBUFFERC3) .EQ . ’- ’ ) THEN 

DECODE C36, 206, BUFFER) DX,DY,DZ 
FORMAT C2X, Ell .4, IX, E10.4, IX, Ell .4) 


207 


ELSE 

DECODE C 36, 207, BUFFER) DX.DY.DZ 
FORMAT C 2X ,E10.4,2C1X,E11.4)) 

END IF 


208 


ELSE 

TYPE 208 

FORMAT C//, IX, ’UECTOR LENGTH OUT OF RANGE’) 
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RETURN 

END IF 

**«******« An individual vertex has been moved. 

*****»***»«» Filter out occasional erroneous transmissions 
IF C ABS C DX ) . GE .0.1) DX-0 . 0 
IF CABSCDY) .GE.0.1) DY-0.0 
IF CABSCD2) .GE.0.1) D2-0.0 

IF CnSG.LE.80) THEN 
N - CnSG+3)/10 
n - nsG-ccnsG-i)/iO)*io 



cpcn, n, i) 

CPCN , N , 8) 
CPCN,N,3) 



CPCN.N, 1) + DX 
CPCn,N,B) + DY 
CPCN,N,3) + D2 



END IF 



**»*«**»«** ^ roU) Qr C Qi um n of vertices has been moved 
IF CnSG . GT . 80) THEN 



810 



880 



IF CnSG.LE.30) THEN 
a row has been 

n - nsG - 80 

DO 810 N - 1 , NCP 
CPCn.N, 1) = 
CPCn.N, 8) - 
CPCN.N, 3) - 
CONTINUE 



moved 



cp c n , n , i ) 

CPCN.N, 8) 
CPCN, N, 3) 



END IF 



IF CnSG.GT.SO) THEN 

a column has been moved 
N - NSG - 30 
DO 880 n - i,ncp 

cpcn.N, i) - cpcn,N, i) 

CPCn,N,8) - CPCn,N I 8) 
CPCn,N,3) - CPCn,N,3) 
CONTINUE 

END IF 



DX 

DY 

D2 



DX 

DY 

D2 



END IF 



***** Report the new B-spline surface and vertices’ 
***** vectors to the PS-300 only if vertex movement 
***** is completed, to minimize transmission time. 

CALL PSPOLLCBO, BUFFER, LENGTH) 

IF C LENGTH. E0.0) THEN 

IF CICP.EQ.8) CALL CPLINE 
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IF CIBS.BT.05 CALL BSPLIN 

END IF 
RETURN 

ENDIF 

STOP 

END 

*********************************************************** 

*********************************************************** 

******************** Subroutine PLANES ******************** 
•******«***••****•********••****«•*•••*******«*••**•******* 

••••*•«•»•*****«•*•««•*•***•*******•*•***•**•**•*****•*•»*• 

SUBROUTINE PLANESC IPLANE1 

***** When implemented, this subroutine will calulate and 
***** transmit to the PS-300 the intersection of the 
***** B-spline surface with the orthogonal plane specified 

***** by i plane. 

RETURN 

END 
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*********************************************************** 

*********************************************************** 

******************* SUBROUTINE PREPS ********-*-*-*-*-*'*-*-*-*-*-*-*-*'* 
************ *********************************** ************ 
*********************************************************** 



* Input 

* 

♦Output 

* 

♦Calls 

♦Alters : 

♦Description : 



* 

m 

* 

* 

* 

* 



♦Created By 
♦Created On 



PREPS-format data file defining 
structure to be displayed. 

PS-300 vector list for display of given 
structure . 

PSSEND , PSUECS , BLANK 

PS-300 program structure KLINES. 

This subroutine reads a data file of points 
known to be on the surface or limiting the 
surface to be modeled with the B-spline. 

The data is interpreted in accordance with 
the format of the program, PREPS. This 
format is utilized at the University of 
Washington as a standard format for all 
digitized data to be input to the PS-300. 
The data is then passed to the PS-300 in th 
form of a vector list for display. 

D. flichael Bryant 

4/20/87 Modified On: 



* 

•* 



* 

* 

* 

* 



e* 

* 



*********** 
* 

*Uariable 



»****#*«******#«***#«»***##*****»****«**»#***»#* 



Explanation 



♦BUFFERC80D 

* 

♦FILNAfl 



♦IC00RD 

•* 

* 



IGNORE 



♦IOSTAT 
♦I REF 



♦ISTAR 

* 

* 

m 



Type 

Logical Contains strings being received from 
or sent to the PS-300. 

Char Name of the PREPS-format data file 

containing the data for the known 
points . 

Integer Coordinate to be held constant for 
each repetition of the substructure 
directed by NPTS2 . 0=none; 1=X; 2=Y ; 
3=2. 

Integer Dummy variable used to read data from 
PREPS-format file which is not needed 
by this program. 

Integer Reports input/output error status. 

Integer Plane about which substructure is to 
be reflected. 0 = none; 1 = YZ; 

2 = X2; 3 = XY. 

Integer Number of vector components required 
to define the position of stars used 
to represent data points from data 
file on the PS-300 screen. A 
maximum of 133 points can be so 
labeled, although the structure can 



* 

•* 

* 

* 

* 

* 

* 

* 

* 
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* 

* I UEC 


I nteger 


contain many more points. * 

Identifies the element of UECS which * 


m 

m 

m 

*LP C 1000 ) 


Logical 


contains the Z vector component of * 
the final vector in the previously * 
processed substructure . * 
Contains ’’point” (.FALSE. 3 or ’’line” * 


* 

* 

*F1AXUEC 


I nteger 


(.TRUE. 3 designations for the * 
itemized vector list defined by UECS.* 
Number of vector components which the* 


* 

*NDONE 


Integer 


current substructure adds to UECS . * 

Ualue greater than zero indicates end* 


* 

*NPTS 


I nteger 


of data file has been reached. * 

Number of paints defining the current* 


* 

*NPTS2 


Integer 


substructure . * 

Number of repetitions of the sub- * 


m 

*NSTRUC 


Integer 


structure to be constructed. * 

Number identifying the current * 


* 

* 

m 

m 

m 

^REPEAT CIO) 


Real 


portion of the known data being * 
processed. To minimize memory * 
consumption, the known data is split * 
into a number of vector lists rather * 
than a single long one. * 
Contains the values of the coordinate* 


* 

* 

*UECS C 3000 D 


Real 


specified by ICDDRD to be used with * 
repetitions of a substructure. * 
Contains the vectors for the vector * 


m 

*XBDOY 


Real 


list to be transmitted to the PS-300.* 
Global translation along the X axis * 


* 

* 

*XNDRf1 


Real 


to be applied to the structure. * 
Applied before normalization. * 
Normalization factor by which all * 


m 

* 

*XTRAN 


Real 


vectors in data file are to be * 
divided. * 
Translation of substructure along the* 


* 

*YBQOY 


Real 


X axis. Applied before normalization.* 
Global translation along the Y axis * 


* 

* 

*YTRAN 


Real 


to be applied to the structure. * 
Applied before normalization. * 
Translation of substructure along the* 


m 

*ZBOOY 


Real 


Y axis. Applied before normalization.* 
Global translation along the Z axis * 


* 

* 

*ZTRAN 


Real 


to be applied to the structure. * 
Applied before normalization. * 
Translation of substructure along the* 


* Z axis. Applied before normalization.* 

*******»********************************>*****•******•****-*-•*-** 



SUBROUTINE PREPS 



connoN /uecs/ uecsc 10003 

DIMENSION REPEAT C 10 3 
LOGICAL LPC3333 
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LOGICAL 1 * 1 BUFFERC 801 
CHARACTER*80 FILNAN 
NSTRUC - 1 
NDONE = 0 
IUEC = 0 

TYPE », ’ ENTER FILE NAME.’ 
READ C5 , 51 FILNAH 
5 FORNAT CAB01 



OPEN CUNIT-B.NAHE-FILNAn, TYPE- ’OLD’ , ERR-2001 



»**•** 



Read global header. 



READ C3, *, ERR-300 1 XNORH 

READ C 3 , * , ERR-300 1 IGNORE 

READ C 3 , * , ERR-300 1 XBODY , YBODY , 2B0DY 



***** Read substructure header and data . 



10 



READ C3, *, END-500, ERR-3001 NPTS , NPTS2 , I COORD 
READ C3, *, ERR-3001 XTRAN , YTRAN , 2TRAN , I REF 



****«•«**• Ensure the UECS array dimension will not be 
********** exceeded by reading the substructure data. 



NAXUEC - NPTS*NPTS2*3 
IF CIREF.NE.01 NAXUEC - NAXUEC*2 
IF C IUEC+NAXUEC.GT. 10001 THEN 
GOTO 400 

END IF 



********** Read substructure data . 



15 CONTINUE 

IF CICOORD.NE.il THEN 

READ C3, *, ERR-3001 C UECS C I UEC+ I *3-21 , I -1 , NPTS 1 
ELSE 

READ C3, *, ERR-3001 C REPEAT C I 1 , I -1 , NPTS2 1 

END IF 



IF C I COORD . NE . 21 THEN 

READ C3,*, ERR-3001 CUECS C IUEC+ I *3-1 1 , I -1 , NPTS1 
ELSE 

READ C3, *, ERR-3001 CREPEATC II, 1-1 , NPTS21 

ENDIF 



IF C I COORD. NE. 31 THEN 

READ C3, *, ERR-3001 CUECSC IUEC+I *31 , 1-1 , NPTS1 
ELSE 

READ C3, *, ERR-3001 CREPEATC II, 1-1 , NPTS21 
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END IF 

TYPE *, ’ Working . ’ 

«**«»***»* Construct substructure repetitions and perform 
********** translations and normalization. 

DO 30 I - 1 , NPTS2 

DO 20 J - l.NPTS 

L - C I-11»NPTS*3 

IF C ICOORD . EQ . 1 1 THEN 

UECSC IUEC+J*3-2+Ll - 

+ C REPEAT C I l+XTRAN+XBODYl/XNORfl-O .5 

ELSE 

UECSC IUEC+J*3-2+L!) - 

+ CUECSC IUEC+J*3-2D+XTRPN+XB0DYJ/ 

+ XN0Rn-0.5 

END IF 

IF C ICOORD .EQ .23 THEN 

UECSC IUEC + J*3 _ 1+L) - 

+ CREPEAT C I 1+YTRAN+YB0DY1 /XNORfl 

ELSE 

UECSC IUEC+J»3-1+LD - 

+ CUECSC IUEC+ J*3-l l+YTRAN+YBODY 1 /XNORN 

END IF 

IF C ICOORD . EQ . 33 THEN 

UECSC IUEC+J*3+L1 - 

+ CREPEATC I 1+ZTRAN+ZBODYl/XNORfl 

ELSE 

UECSC iuec+j*3+li - 

+ CUECSC I UEC+ J*3 1 +ZTRAN+ZBQDY 1 /XNORM 

END IF 

LPC Cl — 1 1 *NPTS+J+ I UEC/31 - .TRUE. 

IFCJ.EQ.il LPC CI-11*NPTS+J+IUEC/31- .FALSE . 

20 CONTINUE 

30 CONTINUE 

********** Construct substructure reflections. 

IF CIREF.GT.O) THEN 

DO 40 I - IUEC, IUEC+flAXUEC/2-3, 3 
DO 35 J - 1,3 

UECSC I+J+NAXUEC/21-UECSC I+Jl 
35 CONTINUE 

UECSC I + IREF+NAXUEC/21 - UECS C I + 1 REF 1 * C -1 . 0 1 
LP C C I +nAXUEC/2 1 /3+1 1 = LPCI/3+11 
CONTINUE 



40 
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END IF 

IUEC - IUEC + NAXUEC 
GOTO 10 

***** Error statements. 

E00 WRITE C5 , 210} FILNAN 

E10 FORNAT (.//,’ Error opening disk file’,/, 

+ 1X.AE0) 

STOP 

300 WRITE C5 , 310) FILNAN 

310 FORNAT (.//, ’ Error reading data from ’ , AEO , / , 

+ ’ Ensure File is Formatted correctly //) 

CL0SEC3) 

STOP 

***** OeFine the position oF stars to represent up to 1SS 
***** data points. 

400 IF CNSTRUC.EO. 13 THEN 
ISTAR - 199*3 

IF C IUEC. LT. ISTAR) ISTAR - IUEC 
DO 450 I - 3, ISTAR, 3 

IF C I- I / 10*10 . EQ . 0) TYPE *,’ Working....’ 
ENCODE C61, 410, BUFFER) 1/3, UECSC 1-3) , 

+ UECSC 1-1 ) , UECSC I ) 

410 FORNAT C ’ KPT ’, 1 3 -TRANSLATE BY ’, 

+ F7.4, ’ , ’ , F7.4, ’ , ’ , F7.4, 

+ ’ APPLIED TO STAR; ’ ) 

CALL BLANKC4, BUFFER) 

CALL BLANK C 5, BUFFER) 

CALL PSSENDCG1, BUFFER) 

450 CONTINUE 

END IF 

***** Send vector list For partial structure to PS-300. 

ENCODE C33,4G0, BUFFER) NSTRUC, IUEC/3 
4G0 FORNAT C ’ KL ’ , I E , ’ : -UECTOR_L I ST ITENIZED N-’,I4) 

CALL BLANK C 3, BUFFER) 

CALL PSSEND C 33, BUFFER) 

CALL PSUECS C4, IUEC/3, UECS.LP, .TRUE. , IOSTAT) 

NSTRUC - NSTRUC + 1 
IUEC - 0 

IF CNDONE.EQ.l) GOTO GOO 
GOTO 15 
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***** The end of the data file has been reached. 

500 CLOSE C3) 

IF CIUEC.GT.03 THEN 
NDONE - 1 
GOTO 400 

END IF 

500 CONTINUE 

RETURN 
END 

4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft 
4ft 4 *4ft4ft4ft4ft*4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft*4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft 

Subroutine REFLCT ***■*■*■*■*■*■*■*■*■*■*•*■*■*■*■*■*■* 

**m*m*m****mm**mi*m**i*m*i****m*m**i***i***i*i*m**m*m**m**mmmm-*m-** 



SUBROUTINE REFLCT 

***** This subroutine will calculate and display a 
***** reflection of the displayed B-spline surface about 
***** the body’s XZ plane. 

RETURN 

END 



4ft*4ft*4ft*4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft-*4t4ft4ft4ft-*-*4ft4ft4ft4ft-*-*4ft*t-*4ft-*4ft4ft-*4ft4fc-*4fc4ft-*4ft4*-*4*4ft4ft-*4ft4ft4ft 

*4**m***m*mm**4*mm***mm*mmi**mm*m*4**i*i****4***mm*4**i**m*****m*i**m 
4ft 4ft 4ft 4ft 4ft 4ft 4ft 4ft 4ft 4ft 4ft 4ft 4ft 4ft 4ft 4ft 4lt 4ft 4ft SUBROUTINE RESET 

4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4K4ft4ft4ft4ft4ft4ft4ft 

4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft 



*Input 

•Output 

•* 

•Calls 

•Alters 

4ft 

•Description : 
* 

* 

m 



•Created By 
•Created On 



IRESET * 
PS-300 program lines resetting rotation, * 
translation, and scaling nodes. * 
PSSEND * 
PS-300 program translation, rotation, and * 
scaling accumulators. * 
When triggered by the appropriate function * 
key, this subroutine resets the translation * 
accumulator to CO, 0,03, the scaling * 
accumulator to 1, and the rotation * 
accumulator to the value indicated by * 
IRESET: 1=C0,0,03, 2=C30,0,03, 3=CS0, SO, 0) . * 
D. Michael Bryant * 
12/17/86 Modified On: * 



4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft4ft 

4ft 4ft 

*Uariable Type Explanation * 

4ft * 



155 



*BUFFERCB03 Logical Contains strings being received from * 

* or sent to the PS-300 . * 

*IRESET Integer Indicates the current value of the * 

* three-way rotation resetting toggle. * 

* 1=C0,0,03, 2=C90 ,0,03 , 3=C90 , 90 , 0} . * 

* * 
*************«********«*********««*****«******************* 



SUBROUTINE RESET C IRESET 3 
LOGICAL*! BUFFERC803 



CALL PSSENDC31 , ’SEND UC0,i 
CALL PSSENDC29, ’SEND UCO, 

CALL PSSENDC 19 , ’SEND 1 TO 
CALL PSSENDC 19, ’SEND 0 TO 

IF C IRESET. GT. 13 THEN 

CALL PSSENDC 22, ’SEND 
ELSE 

CALL PSSENDC21 , ’SEND 

END IF 

IF C IRESET . EO . 33 THEN 

CALL PSSENDC22 , ’SEND 
ELSE 

CALL PSSENDC21 , ’SEND 

END IF 



1,03 TO <2>TRAN_T0TAL; ’ 3 
1,03 TO < 1> ALL . TRAN ; ’ 3 

<2> SCALE; ’ 3 
< 1>SCALE ; ’ 3 

90 TO <1> RESET _X; ’ 3 
0 TO <1>RESET_X; ’ 3 

90 TO < 1 > RESET_Y ; ’ 3 
0 TO < 1 > RESET_Y ; ’ 3 



RETURN 

END 
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A***********************-***-******-************-*******-****-*-** 

Subroutine SURFAC * *■*»***■**»■* ■***■* »»■**■* 
•*••«*•***•*••*«»**•*•»**••**»*********»••**•***•••••*•»«** 



« 

♦Input : 

♦Output : 

♦Calls : 

♦Alters : 

♦Description : 

m 

•* 

* 

♦Author: : 

♦Created On : 



SURFC3,31,315 ♦ 
UEC5C5005 , LPC3335 ♦ 
PSSEND , PSUECS ♦ 
Structure ”BS” in PS-300 program. ♦ 
This subroutine extracts a vector list from * 
the array variable SURF and passes it to ♦ 
the PS-300 For display of the calculated * 
B-spline surface. ♦ 
D. Michael Bryant ♦ 
E/55/87 Modified On: ♦ 



***************************•*•*•*•*•*•*•*•*•*•*•*•*•*•*•*•*-*•*•*•*•*•*•*-*•*•*•*•*•*•**•* 



♦Uariable Type Explanation 

-* 



*BUFFERCBO) 

* 

* 


Logical 


* IBS 

m 

•* 


Integer 


* I LP 

* 


Integer 


*IQSTAT 


Integer 


* IUECS 

m 


I nteger 


*LP( 170} 

* 

* 

•* 

* 

m 

•* 


Logical 


*NU 


I nteger 



♦SURF (3,31 , 315 



•* 


Real 


m 




* 




*UECS C 500 D 


Real 



•* 



Holds a message for transmital to the 
PS-300 program with the PSSEND 
library subroutine. 

Indicates status of display of 
B-spline surface: 0 - not displayed, 

1 = displayed. 

Counter indicates the current array 
element of LP . 

Input/output error status. 

Counter indicates the current array 
element of UECS . 

.TRUE, indicates corresponding 
elements of UECS represent a line 
vector (draw to this point with pen 
down! . .FALSE, indicates 
corresponding elements of UECS 
represent a point vector Cmove to 
this point with pen up5 . 

Total number of vectors in the vector 
list for the B-spline surface. 

The calculated points on the B-spline 
surface, in cartesian coordinates, as 
a function of the parametric 
coordinates . 

Uectar list defining the parametric 
lines representing the B-spline 
surf ace . 



■m 

•* 

* 



m 



* 



* 

* 

* 

* 

* 

•* 

* 

* 

* 



•* 



* 

* 



•* 



I*******************************************-***-*-*-*-*-***-**-*-*-*-* 
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SUBROUTINE SURFAC 

LOGICAL LP(333) 

COMMON /UECS/ UECS(IOOO) 

COnnON /BASIS/ RC13, 4,313, 5(11,4,31), SURF ( 3 , 31 , 31 ) 
LOG I CAL*1 BUFFER (80) 

NSTRUC - 0 


***** 


Parametric lines in u (longitudinal) direction. 


10 


CONTINUE 

NSTRUC - NSTRUC + 1 
NU - 310 
IUECS - 0 
ILP - 0 

NBEGIN - ( NSTRUC- 1 ) *10+1 
NEND - NSTRUC*10 
IF (NSTRUC. EQ. 4) THEN 
NEND - NBEGIN 
NU - 31 

END IF 


* 


DO 100 I - NBEGIN, NEND 
DO SO J - 1,31 

Count the vectors created . 
ILP - ILP + 1 


* 


DO BO K = 1,3 

Count the vector components (x,y,z). 
IUECS - IUECS + 1 


BO 

* 

* 


UECS( IUECS) = SURF ( K , I , J ) 

CONTINUE 

Indicate whether vector is point or line 
vector . 

LP(ILP) - .TRUE. 

IF (J.EO.l) LP(ILP) = .FALSE. 


30 

100 


CONTINUE 

CONTINUE 

GOTO 300 


***** 


Parametric lines in u (transverse) direction. 


110 


CONTINUE 

NSTRUC - NSTRUC + 1 
NU - 310 
IUECS = 0 
ILP = 0 

NBEGIN = (NSTRUC-5)*10+1 
NEND « (NSTRUC-4)*10 
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IF CNSTRUC . EQ . 8) THEN 
NEND - NBEGIN 
NU - 31 

END IF 

DO 800 J - NBEGIN, NEND 
DO ISO I - 1,31 

ILP = ILP + 1 
DO 180 K - 1,3 

IUECS = IUECS + 1 

UECSC IUECS) = SURF C K , I , J ) 

180 CONTINUE 

LP C I LP ) = .TRUE. 

IF CI.EGJ.l) LP C ILP ) = .FALSE. 

130 CONTINUE 

800 CONTINUE 

300 CONTINUE 

ENCODE C31 , 310, BUFFER) NSTRUC,NU 
310 FORNATC ’ BS ’ , II, ’ : **UECTOR_L 1ST ITEMIZED N- ’, 13) 

CALL PSSEND C 31 , BUFFER) 

CALL PSUECS C 4 , ILP , UECS , LP , . TRUE . , IOSTAT) 

TYPE *, ’ Working ’ 

IF CNSTRUC . LT . 4) GOTO 10 
IF CNSTRUC.LT. 8 5 GOTO 110 

RETURN 

END 



****•**»•**•*•*•**»•*•*•*•**•*•*-*-* «**•*•*•*•*•***•*•*•*•*•*•*•*■*■*•*•* *•*•*•*•*•**•*-**»•*•**•*-*•*•* 
****-*•**-****•*•*•**•*•*•*•*•** Subroutine U1AIT 



♦Input : 

♦Output : 

* 

♦Calls : 

m 

“Alters : 

“Description : 

m 

m 

m 



m 



* 



BUFFER, LENGTH 

ICLIP, ICP, I KNOWN , IPICK, IRFLCT , IRESET, 
IXY,IXZ,IYZ, and PS-300 equivalents 
BSPLIN, CPLINE, 10 , HOUING , PLANES , PSREAD , 
PSSEND , REFLCT , RESET 



Awaits signals under the control of the “ 
PS-300 program which report Function key “ 
presses or the picking of vertices For “ 
movement. Determines the mode oF operation “ 
dictated by the key presses, updates the “ 
related variables and inForms the PS-300 * 
program oF the changes, and invokes the “ 
appropriate FORTRAN subroutines For any “ 
resultant actions. * 
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•Created On : 11/25/95 



Modified On: 



* 

* 



*********************************** 



************************ 

* 



•Uariable 

* 


Type 


•BUFFERCE30) 

* 


Logical 


•IBS 

* 


I nteger 


* I CL IP 

*ICP 

* 


I nteger 
Integer 






* I KNOWN 

* 


I nteger 


* 




*IPICK 

* 


I nteger 


* 




* I RESET 

* 


I nteger 


* I RFLCT 

* 


Integer 


* 




*IXY 


INTEGER 



Explanation 



*IXZ 

* I YZ 

^LENGTH 

* 

♦MODE 

* 

* 

* 

* 

* 

* 

•msg 

* 

* 

* 

* 

* 

* 

* 

* 

* 



I nteger 
I nteger 
I nteger 

I nteger 



I nteger 



Contains strings being received from 
or sent to the PS-300 . 

Status of B-spline surface display: 

0 - not displayed, 1 = displayed. 
2-clipping status: l=off, 2=on . 

Control polygon display status: 0-not 
displayed, 1-points displayed, 
2-connecting lines displayed. 

Display status of known surface: 

O-not displayed, 1-points displayed, 
2-connecting lines displayed. 

Picking mode indicator: 1-picking 

individual vertices, 2-picking rows 
or columns of vertices . 

Reset toggle indicator: 1-CXO, YO, 20) , 

2- CX90, Y0, 20) , 3-CX90, Y90.20) . 

Status of display of reflection of 
B-spline surface about an X2 plane: 

0- not displayed, 1-displayed. 

Status of display of intersection of 
B-spline surface with XY orthogonal 
plane Cwaterlines) : 0-not displayed, 

1- coarse spacing, 2-fine spacing. 

See I XY . 

See I XY . 

Length of incoming message from 
PS-300, in bytes. 

Function key availability status: 
1-standard mode, most functions 
available, 2-display of reflection or 
intersection with planes selected, 
interactive mode not available, 

3- interacti ve mode selected, 
4=input/output menu selected. 

Contains a numeric message from the 
PS-300, as decoded from BUFFER. 

1-B0 indicate next message will be a 
new vector for the vertex of the same 
number. 101-242 indicate next 
message will be new vector of the 
vertex in the row or column of the 
control polygon line which has the 
position in the vector list CN-100) . 
301-312 indicate a keypress of 
function keys 1-12. 
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***«*************«*******«**«l«**«|HM«*|Ht«««***« < « << | ll H l «l l< l < |* lt 

SUBROUTINE WAIT 





L0GICAL*1 BUFFERC80) 


mm-*** 


Initialize variables For default operating modes. 

LENGTH = 0 

H0DE-1 

I CP-1 

IBS-0 

IKNOWN-O 

IPICK-1 

ICLIP-1 

I RFLCT-0 

IXY-0 

1X2=0 

IY2-0 

IRESET-3 


***** 


Receive and decode integer message from PS-300 


***** 


program . 


100 

* 


IF CLENGTH . GT . 0} GOTO 130 

LENGTH > 0 indicates message already received in 


* 


another subroutine. 


110 


CALL PSP0LLC80, BUFFER, LENGTH) 
IF CLENGTH . EQ . 0) GOTO 110 


130 

E00 


DECODECLENGTH, 200, BUFFER) MSG 
FORMAT C 13 ) 

LENGTH - 0 


***** 


If message indicates new vertex vector is coming, 


***** 


go to subroutine MOUING to receive and act upon the 


***** 


vector . 

IF CMODE.EQ. 3. AND.MSG.LT. 300) THEN 

CALL MOUING C ICP , I BS , MSG , LENGTH , BUFFER ) 

ENDIF 

IF CMODE.EQ. 3. AND.MSG.LT. 300) GOTO 100 


***** 


Uertex display toggle. 


300 


IF CM0DE.LE.3) THEN 

IF CMSG.EQ.301) THEN 
ICP-ICP+1 

IF CICP.GT.2) I CP-0 
IF CICP.EQ.2) CALL CPLINE 
ENCODE C25, 300, BUFFER) ICP 

FORMAT C ’ SEND FIXC’Il’) TO <1>CPM0DE;’) 
CALL PSSEND CBS, BUFFER) 



ENDIF 
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***** B-spline surface display toggle. 



310 


IF CnSB.EQ.302) THEN 
IBS-IBS+1 

IF CIBS.GT.l) IBS-0 
IF CIBS.ET.O) CALL BSPLIN 
ENCODE C25, 310, BUFFER) IBS 

FORMAT C ’ SEND FIXC’Il’) TO < 1> BSMODE ; ’ ) 
CALL PSSEND C25, BUFFER) 

END IF 


***** 


Known surface display toggle. 


320 


IF C MSB . EQ . 303 ) THEN 
I KNOWN- I KNQUN+1 
IF C IKNOUN .BT .2) IKNOUN-O 
ENCODE C2B, 320, BUFFER) IKNOUN 

FORMAT C ’ SEND FIXC’Il’) TO < 1> KNOUNMODE ; ’ ) 
CALL PSSEND C2B, BUFFER) 

END IF 


***** 


Interactive mode toggle. 


330 


IF CMSB . EQ . BOH) THEN 

IF CMODE.EQ.l) THEN 
MODE-3 

ELSEIF CM0DE.EQ.3) THEN 
MODE-1 

END IF 

ENCODE C23, 330, BUFFER) MODE 

FORMAT C ' SEND FIXC’Il’) TO <1>M0DE;’) 
CALL PSSEND C23, BUFFER) 

END IF 


***** 


2-clipping toggle. 


340 


IF CMSB. EQ. 305) THEN 
ICLIP-ICLIP+1 
IF CICLIP.GT.2) ICLIP-1 
ENCODE C27, 340, BUFFER) ICLIP 

FORMAT C ’ SEND FIXC’Il’) TO < 1>CL IPMODE ;’ ) 
CALL PSSEND C27, BUFFER) 

END IF 


***** 


Reflection display toggle. 


n n n n 


IF CMSB. EQ. 306) THEN 

IF CMODE.EQ.l) THEN 
MODE-2 
IRFLCT-1 



nnnnn * * nnnnnnnnnn 



1B2 



C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 



+ 



CALL PSSENDC23, 

’SEND FIXC23 TO <l>nODE; ’3 
CALL PSSENDC22, 

’SEND 1 TO <l>NEU_f10DE; ’ 3 
CALL REFLCT 
CALL PSSENDC 30 , 

’SEND FIXC 1 3 TO < 1> REFLECTflODE ;’ 3 
ELSEIF CflODE . ED . 2 . AND . IRFLCT . ED . 03 THEN 
I RFLCT-1 
CALL REFLCT 
CALL PSSEND C30, 

’SEND FIXC13 TO < 1> REFLECTflODE ; ’ 3 
ELSEIF CflODE . EQ . 2 . AND . IRFLCT . EQ . 1 3 THEN 
IRFLCT=0 

CALL PSSEND C30, 

’SEND FIXC03 TO < 1> REFLECTflODE ;’ 3 
GOTO 400 

ENDIF 

END IF 



***** XY intersection display toggle. 



IF CfISG .EQ .3073 THEN 

IF CflODE. LT. 33 THEN 
IXY-IXY+1 

IF CIXY.GT.23 IXY-0 
IF CIXY.EQ.13 CALL PLANESC 1 3 
ENCODE C25, 360, BUFFER3 IXY 

FORflAT C ’ SEND FIXC’,11, 
’ 3 TO < 1> XYflODE ; ’ 3 
CALL PSSEND C 25, BUFFER 3 
IF CIXY.EQ.03 GOTO 400 



ENDIF 



*** Point-to-row picking toggle. 

IF CflODE. EQ. 3. AND. IPICK.EQ. 13 THEN 
I P I CK=2 

CALL PSSEND C27, 

+ ’SEND FIXC23 TO < l>PICKf10DE ; ' 3 

ENDIF 

ENDIF 



**** XZ intersection display toggle. 

IF CfISG. EQ. 3003 THEN 

IF CflODE. LT. 33 THEN 
IXZ-IXZ+1 

IF CIXZ.GT.23 IXY-0 
IF CIXZ.EQ.13 CALL PLANESC 2 3 
ENCODE C 25, 370, BUFFER 3 IXZ 
C370 FORNATC ’SEND FIXC ’, II, 

C + ’ 3 TO < 1 > X2f10DE ; ’ 3 
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CALL PSSEND CEB , BUFFER) 
IF CIXY.E0.01 GOTO 400 



C 
C 

C END IF 

***** Row-to-point picking toggle 



IF CNODE . EQ . 3 . AND . IPICK . EQ . S) THEN 
IPICK-1 

CALL PSSEND CE7, 

’SEND FIXC 1 ) TO < l>PICKnODE ; ’ ) 

END IF 



END IF 

***** Y2 intersection display toggle, 



C 

C 

c 

c 

c 

C380 

c 

c 

c 



IF CF1SG . EQ . 309 . AND . MODE . LT . 3 3 THEN 
IYZ-IYZ+1 

IF C I Y2 . GT . E) IYZ»0 
IF CIYZ.EQ.l) CALL PLANESC 3) 

ENCODE CSS , 380 , BUFFER ) IYZ 

FORHAT C ’ SEND FIXC’Il’) TO < 1> YZNODE ; ’ ) 
CALL PSSEND CE5, BUFFER) 

IF C I YZ . EQ . 0 ) GOTO 400 

END IF 



***** Select input/output menu. 

IF CHSG .EQ .310) THEN 
LAST - NODE 
NODE - 4 

CALL PSSEND CE3,’SEND FIXC4) TO <l>nODE;’) 

END IF 

***** Reset function toggle. 

IF C tISG . EQ .311) THEN 
IRESET” IRESET+1 
IF C I RESET . GT . 3) IRESET=1 
CALL RESET C IRESET) 

END IF 



***** Quit . 



IF CnSG.EQ.31S) THEN 

CALL PSSEND Cll, ’INITIALIZE;’) 
RETURN 
END IF 



END IF 



***** i/q menL | selections. 

IF CH0DE.EQ.4) THEN 

IF CnSG . LT . 31E) CALL I0CMSG-300) 
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IF CnSG.EQ.31E) THEN 
NODE = LAST 

ENCODE CE3, 330, BUFFER) NODE 
CALL PSSEND CE3, BUFFER) 

END IF 

END IF 

***** Trigger PS-300 program to send new variable values to 
***** all ’FETCH’ nodes. 

CALL PSSEND CEE, ’SEND 1 TO < 1>NEWJ10DE ; ’ ) 

GOTO 100 

***** IF toggling a Function ’oFF’ while in mode E, ensure 
***** all mode-E Functions are oFF beFore returning to 
***** mode 1 . 

400 IFCIRFLCT.EQ.O.AND. IXY . EQ . 0 . AND . IXZ . ED . 0 . AND . IYZ.EQ.O 
+ ) THEN 

NODE - 1 

CALL PSSEND CE3,’SEN0 FIXC1) TO <l>nODE;’) 

CALL PSSEND CEE, ’SEND 1 TO < 1>NEUM10DE ; ’ ) 

ENDIF 
GOTO 100 



END 



APPENDIX C 



SPLINE . DAT 
program, SPLINE 



SPLINE.DAT PROGRAM LISTING 
is the PS 300-resident portion of 



the 



IBB 



C 

SPLINE . DAT ****+***-*-**-*******-** : ** ; * : #- 
***■**■-*-*-**■****•**-*-**■ Prograrn Organization *********’********’** 



I . Uariables . 

II. flain display tree. 

A. Display vertices. 

B. Display B-spline surface. 

C. Display points on known surface. 

D. Display reflection of all displayed structures 

about a plane . 

E. Display intersection of XY planes with B-spline 
surf ace . 

F. Display intersection of XZ planes with B-spline 
surface . 

B. Display intersection of YZ planes with B-spline 
surf ace . 

III. Connect dials. 

A. Translation. 

B. Rotation. 

C. Scaling. 

IU. Uertex picking network. 

A. Enable picking. 

B. Convert cursor position to picklist. 

C. Convert picklist to integer identifying vertex. 

D. Report pick of row or column to hast computer. 

E. Route translation dial signals to picked vertices 

U. Report vertex movement to host. 

UI . Clipping . 

UI I . Function keys and labels. 

A. Report key press to host computer. 

B. Label function keys and dials. 
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*<t**<***»*****<**«*«*«***<t*«*«****«)H|l***«****«*<**********l| 

•**<**»*•*•*********«*•• Variables ************************ 
•*******§•*******•**»***•***•********•*****•**«•****•***•*• 
*«#»««*»«»«»««««««»#<IHt4»«»«»***«*«*»*»*««»*«*«W*«**«******«* 

} 

SEND 

’SPLINE, A B-SPLINE SURFACE NODELER BY □. NICHAEL BRYANT’ 
TO <1>FLABEL0; CDisplay title message . 3 

O 

TO_HOST : -F : PRINT ; CConvert messages going 

to host to strings} 
CONNECT T0_H0ST<1> : < 1>H0ST J1ESSAGE ; 

SEND TRUE TO <2>T0_H0ST; 



O 

NEUJ_f10DE: “F: NOP; 



O 

UARI ABLE B5M0DE; 



SEND FI X Cl 1 TO <1>BSN0DE; 

FETCHBSMODE : -F : FETCH ; 

SEND 1 BSflODE ’ TO < E>FETCHBSNODE ; 

CONNECT NEU_nODE< 1> : < 1>FETCHBSN0DE ; 

O 

UARIABLE CLIPNODE; CZ-clipping mode: l=off, 

E=on} 

SEND FIXC 1 ) TO <1>CLIPN0DE; CInitial default value} 
FETCHCL I PNODE : -F : FETCH ; 

SEND ’CLIPNODE’ TO < S> FETCHCL I PNODE ; 

CONNECT NEW_N0DE<1> :< 1> FETCHCL I PNODE ; 

O 

UARIABLE CPNODE; CControl point display 

mode: O-not displayed, 
1-points displayed, 
E-connecting lines} 

SEND FIXCE) TO <1>CPN0DE; CInitial default value} 

FETCHCPNODE : -F : FETCH ; 

SEND ’CPNODE’ TO < E> FETCHCPNODE ; 

CONNECT NEW_NQDE< 1> : < 1>FETCHCPN0DE ; 



CConvert floating-point 
numbers to exponential 
format . } 

CNo-operation function 
used as a tie point for 
the triggers of FETCH 
commands for all 
variables, which are 
triggered by the host 
computer sending any 
value to NEUi__NODE} 

CB-spline surface 
display mode: O=not 
displayed, 1-displayed . } 
CInitial default value} 



C} 
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CProvides connection for 
function outputs left 
intentionally unused, 
without receiving a ”no 
connection made to . . . . ” 
error message on screen} 

CKnouin surface display 
mode: O-not displayed, 

1- points displayed, 

2- connecting lines} 

SEND FIXC05 TO < 1 > KNOUJNNODE ; CInitial default value} 
FETCHKNOWNNODE : -F : FETCH ; 

SEND ’ KNOUJNNODE ’ TO < 2> FETCHKNOWNNODE ; 

CONNECT NEW _NODE< 1> :< 1> FETCHKNOWNNODE ; 

O 

UARIABLE MODE; 



SEND FIXC 1 ) TO <1>N0DE; 

FETCHNODE : =F : FETCH ; 

SEND ’NODE’ TO < 2> FETCHNODE ; 

CONNECT NEWJ10DE<1> :< 1> FETCHNODE ; 

O 

UARIABLE PICKNODE; CUertex picking mode: 

1- individual points, 

2- rouis/columns of 
points} 

SEND F I X C 1 } TO <1>PICKN0DE; CInitial default value} 
FETCHP I CKNODE : -F : FETCH ; 

SEND ’PICKNODE’ TO < 2> FETCHP I CKNODE ; 

CONNECT NEUI_NODE< 1 > : < 1> FETCHP ICKNODE ; 

C} 

UARIABLE REFLECTNODE; CRef lection-about-plane 

display status: O-off, 
1-on} 

SEND F I X C 0 } TO < 1> REFLECTNODE ; C Initial default value} 
FETCHREFLECTNODE : -F : FETCH ; 

SEND ’REFLECTNODE’ TO < 2> FETCHREFLECTNODE ; 
CONNECT NEWJ10DE<1> :< 1 > FETCHREFLECTNODE ; 



CFunction key menu mode: 
l=main menu, most 
functions available 

2- display of reflection 
or intersecting planes 
selected, interactive 
mode not available 

3- interactive mode 
selected, mode 2 
displays not available 

4- input/output menu 
selected} 

•CInitial default value} 



UARIABLE DISCARD; 



O 

UARIABLE KNOWNNODE; 



O 
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UAR I ABLE XYHODE, XZNODE , YZNODE ; CDisplay of inter- 
section of B-spline 
surface with orthogonal 
planes C i . e . 
waterlines}: O-not 



displayed, 1-coarsely 
spaced, B-finely spaced! 

SEND FIXCOD TO <1>XYN0DE; CInitial default value! 

FETCHXYNODE : -F : FETCH ; 

SEND ’ XYNODE ’ TO < 2> FETCHXYNODE ; 

CONNECT NEU_NODE< 1> : < 1>FETCHXYN0DE ; 

SEND FIXCOD TO <1>XZN0DE; CInitial default value! 

FETCHXZnODE : -F : FETCH ; 

SEND ’XZnODE’ TO < E> FETCHXZNODE ; 

CONNECT NEUi_NODE<l> : < 1> FETCHXZNODE ; 

SEND FIXCO} TO <1>YZN0DE; CInitial default value! 

FETCHYZNODE : -F : FETCH ; 

SEND ’ YZNODE ’ TO < 2> FETCHYZNODE ; 

CONNECT NEUJ_NODE< 1> : < 1> FETCHYZNODE ; 

C 



m*-**-**-*-m-**m-*-*-**#*-*-*-*-**#***t-m»»»-m-*-*-»»**t**mmm*-*m-»n***-***nn»nnn*-*-m 

****«••**•••*•***«•*•*«•*•<***•****•§•***•••***••**•***•*** 

* * **-*♦-*-* * Nain Display Tree <***»«'***'*'*'***■*■*»'*•*** 

***«««****«*«**f *•*••••***••****••***•*•**•****•««••*•***** 
•*•**•••»•**••*••«••*•******••*«•••**••****•*••*«******•••• 

1 

CL I PON : — SET DEPTH_CLIPPINS OFF APPLIED TO CLIPPER; 
CLIPPER:- Id I NDOUJ X — 1:1 Y — 1:1 FRONT-O 

BACK-O.l APPLIED TO CLIP_INTENSITY; 

CL I P_I NTENS I TY : -SET INTENSITY OFF 1:1 APPLIED TO ALL; 

CDisables depth cueing 
when clipping is on, to 
allow viewing within 
narrow clipping windows! 

DISPLAY CLIPON; 

ALL : -BEGIN_STRUCTURE 

U I EWPORT HOR— 1:1 UERT— 1:1 I NTENS I TY- . 75 : 1 ; 

TRAN: -TRANSLATE BY 0,0,0; 

ROT : -ROTATE 0; 

SCALE: -SCALE BY 1; 

SET .-“SET COND I T I ONAL _B I T 1 ON; 

INSTANCE OF CP_SET , BS_SET , KNOUJN_ SET , PLANE_SET ; 
END_STRUCTURE; 

C 



f 
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************** Display Control Points (Uert ices )«*****»**»* 

A**********************-*-**************'**-****'***-**'**-*-***'*-*.** 



} 

CP_SET : “SET LEUEL_0F_DETA I L TO 1 THEN CP_IF; 

{Node to connect toggle 
to> 

CONNECT FETCHCPf10DE< 1> : < 1>CP_SET ; 

CP_I F : ” I F LEUEL_OF_DETAIL >0 THEN CP; 

{Display vertices unless 
toggled off.} 

CP: -INSTANCE OF CP_.PT , CP_L INES ; {Structure consists of 

two branches: points and 
connecting lines.} 

CP_L I NES : =BEGI N_STRUCTURE 

IF LEUEL_OF_DETAIL =2; 

{Display if toggled on.} 
SETP I CK : “SET PICKING OFF; 

{Create picking node to 
allow vertices to be 



moved a row or column at 
a time by picking the 
connecting line} 

SET PICKING I DENT I F I ER“CPSS ; 

INSTANCE OF CPL ; {CPL is a vector list 

which will be provided 
dynamically by the host 
computer . } 

CPL: “UECT0R..LIST ITEMIZED N“200 

0.0,0. 0,0.0; {Reserve memory for 

vector list.} 



END_STRUCTLJRE ; 

{} 

CP_.PT : “BEGIN _STRUCTURE 

SETPICK:“SET PICKING OFF ; {Create picking node to 

allow vertices to be 



INSTANCE OF 

CP01 , CP02 , CPOG , 
CP0S,CP10,CP11 
CP17,CPia,CP19, 
CP25 , CP2S , CP27 
CP33 , CP34 , CP3S 
CP41 , CP42, CP43 
CP49, CP50, CP51 
CP57 , CP50 , CP59 
CP65 , CPGG , CP67 
CP73 , CP74 , CP75 
END_STRUCTURE; 



selected in interactive 
mode by picking with 
tablet and stylus} 

CP04 , CP05 , CPOB , CP07 , CPOB , 
CP12 , CP13 , CP14 , CP15 , CPIS , 
CP20 , CP21 , CP22 , CP23 , CP24 , 
CP2S , CP29 , CP30 , CP31 , CP32 , 
CP36 , CP37 , CP3B , CP39 , CP40 , 
CP44 , CP45 , CP4S , CP47 , CP4B , 
CP52 , CPS3 , CP54 , CP55 , CP56 , 
CPGO , CPG1 , CPG2 , CP63 , CP64 , 
CP6B , CPG9 , CP70 , CP71 , CP72 , 
CP7S , CP77 , CP7B , CP79 , CPBO ; 
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C 

Definitions of CP01 through CP80 will be supplied by the 
host computer, with the necessary nodes and connections to 
allow the vertices to be individually manipulated in the 
interactive mode. The program lines supplied dynamically 
from the host computer will have the following form Cwith 
the 01 replaced with the appropriate number for the nodes 
associated with Bach vertex): 



*CP01 : =BEGIN_STRUCTURE 

* SET PICKING IDENTIFIER-CP01 

« 

* 

* TRAN: -TRANSLATE BY 0,0,0 

* 

m 

m 

m 



m 



* 

* INSTANCE OF CUBE 

* END STRUCTURE 



m 

CAllows identif icai ton * 
of which node was * 

picked for interaction.* 
CThis node will connect* 
through an accumulator * 
and switching network * 
to the dials for * 

manipulation of the * 
picked vertex. Initial* 
placement of the vertex* 
will also be specified * 
through the network * 
connected to this * 

node . ) * 

* 

* 



*TRAN_T0TAL_01 -F : ACCUMULATE 



(Create translation 
accumulator for each 



vertex . ) 

CONNECT TRAN TOTAL 01< 1> : < 1> CP01 . TRAN 



* 



* 



* 

* 



CUBE 



SEND U ( X , Y , Z ) TO < 2> TRAN__T0TAL_Q1 

(Initial vertex 
position . ) 

SEND U(0,0,0) TO < 1> TRAN_T0TAL_01 

(Give accumulator an 
initial value to act 
upon before input from 
dial is received) . 

SEND 0.1 TO < 4> TRAN_T0TAL _01 (Multiplier for dial 

sensitivity . ) 

-INSTANCE OF FRONT , BACK, SIDES ; CEach vertex is 

represented on the 
screen by a cube . ) 

BACK: -TRANSLATE BY 0,0,0.01 APPLIED TO FRONT; 

FRONT : -UECT0R LIST Z— 0.005 N-5 

-0.005,-0.005 0.005,-0.005 0.005,0.005 

-0.005, 0.005 -0.005,-0.005; 



* 



* 

* 

* 

* 

* 

* 

* 

m 
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C 



SIDES: =UECTOR__LI ST SEP N-S 

-0.005, 0.005,0.005 
-0.005,-0.005,0.005 
0.005, 0.005,0.005 
0.005,-0.005,0.005 



-0.005, 0.005,-0.005 
-0.005,-0.005,-0.005 
0.005, 0.005,-0.005 
0.005,-0.005,-0.005; 



**•••**••**«•******••••************•«**•**<•••**•**••••***•* 

**************** Display B-Splme Surface +*+**++■*****«*+** 
******************************************* 1 **************** 



} 

BS_SET : —SET LEUEL_0F_DETA I L TO 1 THEN BS._IF; 

CDisplay B-spline 
surface if toggled on . 1 
CONNECT FETCHBSt10DE< 1> : <1>BS_SET; 

BS_IF:=IF LEUEL_OF_DETAIL >0 THEN BS ; 

BS : “INSTANCE OF BS1 , BS2 , BS3 , BS4 , BS5 , BS6 , BS7 , BS8 ; 

C 

The host computer will dynamically provide program lines 
defining BS . These commands will be of the form: 



BS1 : -UECTQR LIST ITENIZED N-XXX 



} 

C 

*********************************************************** 



************* Display Points on Known Surface ************* 
*********************************************************** 



KNOWN _SET : -SET LEUEL_0F _DETA I L TO 1 THEN KN0WN.._IF; 

CDisplay known points 
if toggled on.} 

CONNECT FETCHKNQWNf1QDE< 1> : < 1> KN0WN_SET ; 

KNOWN_IF:=IF LEUEL_OF_DETA I L >0 THEN KNOWN; 

KNOWN: -INSTANCE OF KN0WN_PTS , KN0WN_L INES ; 

C} 

KNOWN _PTS:- INSTANCE OF 

KPT001 , KPT002 , KPT003, KPT004 , KPT005, KPT005, 

KPT007 , KPT008, KPTOOS, KPT010, KPT011, KPT012, 

KPT013 , KPT014, KPT015, KPT016, KPT017, KPT01S, 

KPT013 , KPT020 , KPT021, KPT022, KPT023, KPT024, 

KPT025 , KPT025 , KPT027, KPT028, KPT02S , KPTOES, 

KPT030 , KPT031 , KPT032, KPT033, KPT034, KPT035, 

KPT03G , KPT037 , KPT03S , KPT039, KPT040, KPT041 , 

KPT042 , KPT043 , KPT044, KPT045, KPT045, KPT047, 

KPTO40, KPT04S , KPT050, KPT051 , KPT052, KPT053, 

KPT054 , KPT055 , KPT055, KPT057, KPT058, KPT05S , 

KPT050 , KPT0E1 , KPT052, KPT053, KPT054, KPT055, 

KPT06S , KPT057, KPT0G8 , KPT0S3, KPT070, KPT071 , 

KPT072 , KPT073 , KPT074, KPT075, KPT076, KPT077, 

KPT078 , KPT078, KPTOBO, KPT091, KPT0B2 , KPT0B3 , 

KPT084 , KPT085 , KPT086, KPT087, KPT088, KPT088 , 
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C 



KPT090 , 
KPT09S , 
KPT102 , 
KPT108, 
KPT114 , 
KPT120 , 
KPT12S , 
KPT132, 
KPT138, 
KPT144 , 
KPT150 , 
KPT15B , 
KPT162, 
KPT1BB, 
KPT174, 
KPT180, 
KPT1BS, 
KPT132, 
KPT198 , 



KPT091 , 
KPT097 , 
KPT103 , 
KPT109 , 
KPT115, 
KPT121 , 
KPT127 , 
KPT133, 
KPT139, 
KPT145 , 
KPT151 , 
KPT157 , 
KPT1S3 , 
KPT169 , 
KPT175 , 
KPT181 , 
KPT1B7 , 
KPT193 , 
KPT199; 



KPT092, 
KPT098 , 
KPT104 , 
KPT 110, 
KPT11E, 
KPT122, 
KPT12B , 
KPT134 , 
KPT140 , 
KPT14B, 
KPT152, 
KPT15B , 
KPT1B4 , 
KPT170 , 
KPT17B , 
KPT1B2 , 
KPT1BB , 
KPT194 , 



KPT093 , 
KPT099 , 
KPT105 , 
KPT111 , 
KPT 117, 
KPT123 , 
KPT129 , 
KPT135 , 
KPT141 , 
KPT147 , 
KPT1B3 , 
KPT159, 
KPT1B5 , 
KPT171 , 
KPT177 , 
KPT1B3 , 
KPT1B9 , 
KPT19B , 



KPT094 , KPT095 , 
KPT100 , KPT101 , 
KPT10B , KPT107 , 
KPT 112, KPT113, 
KPT11B , KPT119, 
KPT124 , KPT125 , 
KPT130 , KPT131 , 
KPT13B , KPT137 , 
KPT142 , KPT143 , 
KPT14B , KPT149, 
KPT154 , KPT15S , 
KPT160 , KPT1B1 , 
KPT1BB , KPT1B7 , 
KPT172 , KPT173, 
KPT17B , KPT179 , 
KPT1B4 , KPT185 , 
KPT190 , KPT191 , 
KPT19B , KPT197 , 



Structures KPT001 through KPT199 mill he provided 
dynamically by the host computer. The commands provided 
will be of the form: 



* 


KPT037 : 


-TRANSLATE BY 


x , y , z APPLIED TO 


STAR 


J 

STAR: 


-UECTOR 


..LIST SEP N-14 










0.000, 


0.005, 0.000 


0.000, 


-0.005, 


0.000 




-0.005, 


0.000, 0.000 


0.005, 


0.000, 


0.000 




-0.004, 


0.004, 0.000 


0.004, 


-0.004, 


0.000 




-0.004, 


-0.004, 0.000 


0.004, 


0.004, 


0.000 




0.000, 


-0.004,-0.004 


0.000, 


0.004, 


0.004 




0.000, 


0.004,-0.004 


0.000, 


-0.004, 


0.004 




0.000, 


0.000,-0.005 


0.000, 


0.000, 


0.005; 


O 












KNOWN 


l_L I NES : 


-IF LEUEL_OF_DETAIL -2 


THEN KLINES; 



CDisplay lines 
connecting known points 
if toggled on . 1 

KLINES: -INSTANCE OF KL01 , KL02 , KL03 , KL04 , KL05 , KLOB , KLG7 , 

KLOB , KL09 , KL10 , KL1 1 , KL12 , KL13 , KL14 ; 

C 

KLINES, the vector list defining the lines connecting the 
known points, will be provided dynamically by the host 
computer. The command provided by the host will be of the 
form: 



* KL01 : -UECTDR_LIST ITEMIZED N-xxx 
C 
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****************«********«**********»-!»*-»**.»W«'*-*W.»W*-*-*-»-»-«-i»4t-* 

************ Display Reflection About a Plane ************* 

*iiai«****i>»****i>******»*<iHiil*<«*#iif««****it**«**«***»ili**f«ji««* 

REFLECT_SET : “SET LEUEL_OF._DETA I L TO 0 THEN REFLECT_IF ; 

{Display reflection 
if toggled on . } 

CONNECT FETCHREFLECTnODE< 1> : < 1> REFLECT_SET ; 

REFLECT_IF: -IF LEUEL_OF_DETAIL >0 THEN REFLECT; 

{ 

REFLECT, the vector list defining the reflection of all 
displayed structures about a plane, will be provided 
dynamically by the host computer. The command provided by 
the host will be of the form: 

* REFLECT: =UECT0R_LIST ITEMIZED N-xxx * 

1 

{ 

*********************************************************** 

*********** Display Intersection With XY Planes »******»«*» 
*********************************************************** 

XY_SET : -SET LEUEL_OF_DETA I L TO 0 THEN XY._IF; 

{Display intersection 
if toggled on . 1 

CONNECT FETCHXYf10DE< 1> : <1>XY_SET; 

XY__I F : — I F LEUEL_OF_DETAIL >0 THEN XY; 

O 

XY: -INSTANCE OF XYC,XYF; {The sturucture has two 

levels of display 
detail: coarse and fine} 
XYF : — I F LEUEL_OF_DETAIL -2 THEN XYFINE ; 

{ 

The structures XYC and XYFINE will be provided dynamically 
by the host computer. The commands provided will be of the 
form : 

* XYC: -UECTOR_LIST ITEMIZED N-xxx * 

* XYFINE:- UECTOR LIST ITEMIZED N-xxx * 

} 

{ 

*********************************************************** 
*********** Display Intersection With XZ Planes *********** 

* *»«**«****««*»l«*«*K*«««*«**lll*»«*««<l>««*«***«<li«»****»»**»»** 

} 

XZ_ SET : —SET LEU£L_0F _DETA I L TO 0 THEN XZ„IF; 

{Display intersection 
if toggled an.} 

CONNECT FETCHXZMODE< 1> : <1>XZ_SET; 

XZ_ I F : = I F LEUEL_OF,_DETAI L >0 THEN XZ; 

{} 
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X2: “INSTANCE OF XZC.XZF; CThe sturucture has two 

levels of display 
detail: coarse and Fine} 
XZF : — I F LEUEL_OF_DETA I L -5 THEN XZFINE; 

C 

The structures XZC and XZFINE mill be provided dynamically 
by the hast computer. The commands provided will be of the 
f orm : 



* XZC: “UECTOR_LIST ITEMIZED N-xxx * 

* XZFINE:- UECTOR_LIST ITEMIZED N=xxx * 

} 

C 

**•***••****••*****•••****•****•**»**•*•****•****••*•**•••* 

*********** Display Intersection With YZ Planes »*»»****»»* 
*********************************************************** 

} 

YZ_SET : “SET LEUEL_OF _DETA I L TO 0 THEN YZ_IF; 

CDisplay intersection 
if toggled on.} 

CONNECT FETCHYZMODE< 1> : < 1> YZ_SET ; 

YZ_I F : = I F LEUEL_OF_DETAIL >0 THEN YZ; 

O 

YZ: “INSTANCE OF YZC,YZF; CThe sturucture has two 

levels of display 
detail: coarse and fine} 
YZF : = IF LEUEL_OF_DETAIL -E THEN YZFINE; 

C 

The structures YZC and YZFINE will be provided dynamically 
by the host computer. The commands provided mill be of the 
form: 

* YZC: “UECTOR _LIST ITEMIZED N-xxx * 

* YZFINE:- UECTOR_LIST ITEMIZED N-xxx * 

} 

C 
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*********************************************************** 

••******«f«*******«*«»*i**««*«**«<«***«««i*f*«***»«*«*»»«i 

********************** Connect Dials ********************** 
*********************************************************** 

*********************************************************** 

! 

c 

(•«***f***»t«**t»t******»**«t**»***l«***t*i***************ii> 

******************** Translation Dials ******************** 
*********************************************************** 



************ Convert dial signals to vectors. ************* 
i 

TRAN_X : -F: XUEC; 

TRAN_Y: =F: YUEC; 

TRAN_2:=F:2UEC; 

CONNECT D I ALS< 1 > : < 1 > TRAN_X ; 

CONNECT DIALS<2> : <1>TRAN_Y; 

CONNECT DIALS<3> : < 1>TRAN_2; 

C 

***** Create switching network which applies ***** 

***** translation dials to global translation Cif in ***** 
***** modes 1 or 2D, vertex translation Cif in mode s,***** 
***** interactive mode} , or no translation Cmode 4, ***** 

***** I/O menu} ***** 

! 

SU I TCH_TRAN : -F : CROUTE Cl}; 

CONNECT TRAN_X<1> : < 2> 50 I TCH_TRAN ; 

CONNECT TRAN__Y< 1> : < 2> SU I TCH __TRAN ; 

CONNECT TRAN_2<1> : < 2> SU I TCH_TRAN ; 

CONNECT FETCHt10DE<l> : < 1 > SU I TCH_TRAN ; 

C 

***** Connect switch to global translation node if in ***** 
***** mode 1 or 2. ***** 

! 

TRAN_TOTAL : -F : ACCUMULATE ; C Accumulator for global 

translation . ! 

CONNECT SUITCH TRAN< 1> : < 1> TRAN_TOTAL ; 

CONNECT SUITCH_TRAN< 2> : < 1>TRAN_T0TAL ; 

CSUITCH_TRAN<3> will be connected later! 

CONNECT SU I TCH_TRAN< 4> : < 1 > D I SCARD ; 

CONNECT TRAN_TOTAL< 1> : < 1> ALL . TRAN ; 

O 

SEND UC0,0,0} TO < 2> TRAN.JTOTAL ; 

CReset value for 
translation . ! 

SEND 1 TO < 4> TRAN_ TOTAL ; CSignal multiplier for 

dial sensitivity! 

SEND FIXC 1 } TO < 1 > SU I TCH_TRAN ; ^Default initially to 

mode 1! 
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Switch will be connected to control point translation 
nodes at a later point in the program.! 



*********************************************************** 
**•*********-•**•••-•**• Rotation Dials *******-*********-*-*-»■** 

*»-*»*»**-*-*j***#'*Wj*i*i*i****j*i***»*j*i**j#*»'l#i***'ilHIHIHIH**H«HIHIHIH**HimMHIH**H<H* 

! 

XMLIL : ”F : NULC ; 



YnUL : =F : 

ZnUL : “F : 

CONNECT 

CONNECT 

CONNECT 

O 

ROTX : ”F : 



CNultipl iers for dial 
sensitivity! 



MULC; 

MIILC; 

DIALS<5> 

DIALS<6> 

DIALS<7> 

XROTATE; 



< 1> XMLIL ; 
< 1 > YMUL ; 

< 1> ZMUL ; 



CConvert dial signal to 
rotation! 



ROTY : "F : YROTATE ; 
ROTZ: ”F : ZROTATE ; 

CONNECT XnUL<l> 
CONNECT YNUL<1> 
CONNECT ZNUL<1> 



< 1> ROTX ; 
< 1> ROTY ; 
< 1> ROTZ ; 



O 

ROT 



C ! 



ACCUM : »F: CNUL; 



CONNECT R0T_ACCUf1< 1> : < 1> ALL 
CONNECT ROT ACCUN< 1> : < 1 > ROT 



C Accumulator 
rotations! 

. ROT ; 

_ACCUn ; 



f or 



RESET _X : =F : XROTATE ; 



CUpon signal from host 
computer, resets 
accumulator with given X 
rotation and 0 For Y and 
Z.T 

CONNECT RESET_X<1> : < 1> ROT_ACCUM ; 

RESET_Y: =F: YROTATE; CUpon signal From host 

computer, adds Y 
rotation reset value to 
accumulator and triggers 
accumulator to Fire.! 

CONNECT RESET_Y<1> :<2>R0T_ACCUN; 

C 

***** Create switching network For ROTATE Function ***** 
***** which will route dials to global rotation node ***** 
***** in all modes Cl-33, except when in interactive ***** 
***** moc j e £43 > i n which case no transformations are 
***** allowed. 

1 






'*<*<**-* 
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SWITCH_ROT: -F: CR0UTEC4! ; 

CONNECT R0TX<1> :<2>SWITCH_R0T; 
CONNECT ROTY< 1> : <2>SWI TCH_ROT ; 
CONNECT ROTZ< 1> : <2> SWITCH_ROT ; 
CONNECT SUJ I TCH_ROT < 1> : <2> ROT_ACCUH ; 
CONNECT SWITCH_R0T<2> : < 2> ROT _ACCUf1 ; 
CONNECT SUJ I TCH_ROT< 3> : < 2> ROT_ACCUN ; 
CONNECT SWITCH _R0T<4> :<1>DISCARD; 
CONNECT FETCHNODE< 1> : < 1> SUJI TCH_ROT ; 
C 



********** send initial values 

SEND 200 TO <2>Xf1UL; 

SEND 200 TO <2>Yt1UL; 

SEND 200 TO <2>2NUL; 

SEND FIXC1! TO < 1> SUI I TCH_ROT ; 

SEND P13DC 1,0,0 0,1,0 0,0,11 TO 

C 



to rotation network. ******* 

{Multiplication factor 
for dial sensitivity! 



CDefault initially to 
mode 11 
< 1> R0T_ACCUf1 ; 

CPrime accumulator! 



*********************************************************** 



*********************** Scaling Dial **************-****-**** 
*********************************************************** 



1 

SUIITCH_SCALE: -F: CR0UTEC4) ; CSee SUIITCH__ROT above 

for description! 

SCALE : “F : DSCALE ; 

CONNECT D I ALS< 4> : < 2> SUJ I TCH_.SCALE ; 

CONNECT SUJ I TCH_SCALE< 1> : < 1>SCALE; 

CONNECT SUJ I TCH_SCALE< 2> : < 1 > SCALE ; 

CONNECT SUJ I TCH_SCALE< 3> : <1>SCALE; 

CONNECT SUJITCH_SCALE<4> : <1>DISCARD; 

CONNECT SCALE<1> : <1> ALL .SCALE; 

Cl 



SEND FIXC11 TO <1>SUJITCH SCALE: 



SEND 1 TO <2> SCALE; 
SEND 1 TO <3> SCALE; 

SEND 20 TO <4> SCALE; 

SEND 0.1 TO <5> SCALE; 



CDefault initially to 
mode 11 
CReset value! 

CMul tipi ication factor 
for dial sensitivity! 
Criaximum scaling 
allowed! 

CMinimum scaling 
allowed! 



Cl 

CONNECT FETCHMODE< 1> :< 1> SWITCH.. SCALE ; 



C 
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^^^^^^^^^^^^^^^^^^^-tt*-*-**-***-*-**-***-*-*-*-*-****-*-*-****-***-**-*-*-****-** 

*******•*•*-*-**-*•*•*-*-* Uertex Picking N s t uj □ r k **■*-*-*■*-*'*•*-*■**-*-*-*•*-*-#• 



**********«****«*******•»•*•*•*•*•*•*•*•*•*•*•**•*•*•**-*-*•***•»•*•*•*•**•*•*-*•»*•»** 

> 

c 

***** Enable picking when stylus tipswitch is pressed.***** 
1 

PICK_IF_N0DE3: =F:CEOUTEC43 ; 

CONNECT FETCHM0DE<1> : < 1> PICK_IF J1QDE3 ; 

CONNECT TABLET I N<4> : < 2> P I CK_I F J10DE3 ; 

CSends ’TRUE’ when 



tipswitch pressed! 

CONNECT PICK_IF_J10DE3< 1> : <1>DISCARD; 

COnly regard stylus 
in interactive made 



if 

C3)! 



CONNECT PICK_IFJ10DE3<2> :< 1> DISCARD ; 

CONNECT P I CK_I F_N0DE3< L t> : <1>DISCARD; 

O 

PICKN0DE_T0GGLE1:=F.-CR0UTEC2! ; CDirects ’true’ tc 

proper set_picking node, 
depending on whether in 
single point or raw pick 
modes! 



C 



CONNECT FETCHP I CKf10DE< 1> : < 1>PICKN0DE_T0GGLE1 ; 

CSee variables! 

CONNECT P I CK_I F_f10DE3< 3> : <2> P I CKN0DE_T0GGLE1 ; 

CPasses ’TRUE’ to SET 
PICKING node only if in 
mode 3! 

CONNECT PICKN0DE_T0GGLE1< 1> : < 1> CP_PT . SETP ICK ; 

CONNECT PICKN0DE_T0GGLE1<B> : < 1>CP_LINES . SETP ICK ; 



***** Convert cursor 
•*•*•*** 

! 



position to picklist 
the component picked 



identifying 



•»•***•»* 

•*•*•*•** 



CONNECT TABLETIN<6> : < 1>PICK; CSend cursor coordinates 

to ’PICK’ for conversion 
to picklist! 

CONNECT PICK<2> : <2> PICKMODE_ TOGGLE 1 ; 

CSend ’FALSE’ to 
SET_P I CK I NG node to 
disable picking once 
pick occurs! 

CONNECT P I CK< 3> : < 2> P I CKMODE_TOGGLE 1 ; 

CDisable picking if 
timeout occurs without a 



pick! 



1B0 



SEND FALSE TO <2>PICK: 



SEND FIXC BO ! TO <3>PICK; 



CCoordinates of picked 
point not required to be 
reported! 

CSet timeout limit for 
when no pick occurs! 



C 

**-***- 

*--*-*•** 



Convert picklist info to integer 
the vertex picked. 



identifying 



*<*•*•*■* 



PICK _ ID: ”F : PICKINFO; 

CONNECT P I CK< 1> : < 1 > P I CK_I D 
C! 

PICKNODE T0GGLE2: -F: CR0UTEC2! : 



CExtracts pick ID Ci.e, 
’CP32’! from picklist! 
CSend picklist! 



CONNECT 

CONNECT 

CONNECT 



CSends pick index if 
picking a row of 
vertices, discards it if 
picking individual 
vertex . 1 

FETCHP I CKNODE< 1 > : < 1> P I CKNODE _T0GGLE2 ; 

PICK_ID< 1> : < 2> P I CKN0DE_T0GGLE2 ; 

P I CKNQDE_TQGGLE2< 1 > : < 1 > D I SCARD ; 



O 

PICKNODE 



CONNECT 

CONNECT 

CONNECT 

CONNECT 

CONNECT 

CONNECT 

CONNECT 

CONNECT 

CONNECT 

C! 

PICK ASCII 



T0GGLE3 : ”F : CROUTE C 2 ! ; CSends pick ID if 

picking individual 
vertex, discards it 
picking a row of 
vertices . 1 

FETCHP I CKMODE< 1> : < 1> PI CKN0DE_T0GGLE3 ; 
PICK_ID<2> : <2>PICKn0DE_T0GGLE3; 

P I CKM0DE_T0GGLE3< 2> :<1> DISCARD; 

P I CK_I D< 3> : < 1 > D I SCARD ; 

P I CK_I D< 4> : < 1> D I SCARD ; 

PICK_ID<5> : <1> DISCARD; 

P I CK_ I D<G> : < 1> D I SCARD ; 

PICK_ID<7> : <1> DISCARD; 

PICK_ID<8> : <1> DISCARD; 



if 



■F: CHARCONUERT; 



CONNECT PICKNODE T0GGLE3< 1> 



SEND TRUE TO < 2> P I CK _ASC I I 
C! 

DIGITS: =F: SUBC; 



CConverts characters of 
pick ID to stream of 
integers representing 
ASCII equivalents! 

< 1>PICK__ASCI I ; 

CProvides pick ID for 
conversion to ASCII! 



CConverts ASCII codes 
pick ID characters to 
numeric digits.! 



of 



lai 



CONNECT PICK_ASCI I < 1> : <1>DIGITS; 

{Provides ASCII-coded 
pick ID For conversion} 
SEND FIXC485 TO <2>DIGITS; {Subtract 4B From ASCII 

code} 

O 

SUJITCH_DIGIT : "F : ROUTE C 4} ; {This switch will send 

each oF the Four 
numerals From the 
converted pick ID For 
separate processing} 
CONNECT DIGITS< 1> : < 2> SW I TCH_D I G I T ; 

{Provides the pick ID 
digits to the switch} 

{ 

The Following is a toggling network which provides an input 
oF integers which progress From 1 to 4, then start again, 
to the switch. This causes each oF the Four successive 
digits From the converted pick ID to be routed to the 
corresponding output port For processing. 

} 

TE I GGER_TOGGLE : -F : CONSTANT ; 

{For each digit 
received oF the pick ID, 
sends a 1 to ADD_ONE} 
SEND FI X C 1 } TO < 2> TR I GGER_TOGGLE ; 

CONNECT PICK_ASCI I<1> : < 1>TRIGGER_T0GGLE ; 

{Provide a signal to 
trigger the toggle} 

ADD_ONE : =F : ADD ; {Provides an integer 

which increases by 1 
each time a digit is 
sent oF the pick ID} 
CONNECT TRIGGER_TOGGLE< 1> : <1>ADD_0NE; 

{Supply the 1 to be 
added each time} 



CONNECT ADD_ONE< 1> : < 2> ADD_QNE ; 

{Supply the sum back to 
the input For 
incrementing next time} 

SEND FIXC3) TO <2>ADD_0NE; 

{Prime the summing 
Function} 

{} 

CYCLE0_3: -F:MODC; {Converts the constantly 

increasing output From 
ADD_ONE to an integer 
cycling From 0 through 
3} 

CONNECT ADD_0NE<1> : <1>CYCLE0_3; 

SEND FIXC45 TO < 2> CYCLE0_3 ; 
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O 



O 



o 



o 



o 



CYCLE1 _4 : =F : ADDC ; {Converts 0-3 cyclic 

output from previous 
function to 1-4 cyclic 
output! 

CONNECT CYCLE0_3< 1> : < 1>CYCLE1_4 ; 

SEND FIXC 1 ) TO <2>CYCLE1_4; 

CONNECT CYCLE1_4< 1> : < 1>SUIITCH_DIGIT ; 

{Provide integers 1-4 to 
control switch! 

CONNECT SWITCH _DIGIT< 1> : < 1>DISCARD; 

{Discard the ”CP” 
pref ix! 

CONNECT SUITCH_DIGIT<2> : <1>DISCARD; 

TENS_DIGIT : -F: MULC; {Multiplies tens digit 

of ID number by 10! 

SEND FIXC 103 TO < 2> TENS_D I G I T ; 

CONNECT SUIITCH_DIGIT<3> : < 1>TENS_DIGIT ; 

{Provides first of two 
digits of ID number for 
conversion to tens.! 

PICKED_CP: =F: ADD; {Reconstructs the ID 

number of the picked 
vertex from the separate 
ones and tens digits! 

CONNECT TENS_DI G I T< 1> :< 1> PICKED CP; 

CONNECT SU1ITCH_DIGIT< 4> : <2>PICKED_CP ; 

PICKED: “F: NOP ; {Provides a node for 

other other portions of 
the program to obtain 
the picking ID from, and 
for the host computer to 
send it to in some 
cases . ! 

CONNECT PICKED._CP<1> :<1> PICKED; 



{ 
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***** 

***** 

***** 

***** 

***** 



If picking rouus and columns of vertices, pass 
the pick index to the host computer for 
conversion to a number identifying the row or 
column to be connected to the translation 

dials . 



***** 

***** 

***** 

***** 

***** 



1 

□FFSET_PICK: -F : ADDC; CAdds offset to pick 

index to put it into 
range which the hast 
program recognizes as 
identifying a pick 
index . 1 



CONNECT P I CKNDDE _T0GGLE2< 2> : < 1 > DFFSET_ PICK; 
SEND FIXC 1005 TO <2>0FFSET_PICK ; 

CONNECT 0FFSET_PICK< 1> : <1>TD_H0ST; 



After conversion by the host to a number identifying 
the row or column picked, the identifier is returned 
from the host in the form: 



1 

C 

***** 

***** 

***** 

***** 



SEND FIXC D TO <1>PICKED 



Establish a switch which directs the output 
of the translation dials, when in mode 3, 
to the translation node of the selected 

vertex . 



* 



***** 

***** 

***** 



! 

SLJ I TCH _P I CK : =F : CROUTE C 9B! ; CA switch with 98 

positions. The first 80 
are for the BO vertices, 
the next ten for the ten 
raws of vertices, and 
the last eight are for 
the columns of vertices! 



CONNECT SUITCH_TRAN<3> : < 2> SU I TCH _P I CK ; 

CProvide the dial output 
to the switch if in 
interactive mode C3!.! 

CIn the program section which reports vertex movement 
to the host computer, the value in PICKED will be 
reported to < 1>SUJITCHPICK to provide the identifying 
number of the selected vertex to the switch.! 



The FORTRAN program on the host computer will provide 
PS-300 commands which will connect the outputs of the 
switch to the inputs of the translation accumulators 
of the vertices. These commands will be of the form: 



* 

1 



CONNECT SUIITCH_PICK<32> : < 1> TRAN_T0TAL _32 
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( 

********* Report Uertex flovement to Host Computer ********* 

»m**m*#m***»****m***»-m*»-m-m*m*-m-m*m*m*-mi*mm****m»*it’mmi*m»»m-#-m-mm 



*********************************************************** 



CP_MOUE: -F: ACCUMULATE; {Accumulates translation 

dial signals if in mode 
3 (interactive model. ! 
CONNECT SW I TCH _TRAN< 3> : < 1> CP_MOUE ; 

SEND UCO.O.Ol TO 



SEND 0.001 



SEND 0.1 



TO 



TO 



RESET CP MOUE :-F: CONSTANT: 



<2>CP_M0UE; 

<3>CP_M0UE; 

(Report accumulation at 
this interval . 1 
<4>CP_M0UE; 

(Scale signal to match 
scaling at individual 
vertex accumulators.! 

_ _ (Resets accumulator to ( 

each time the 
accumulator Fires, so 
only delta values are 
sent to host computer . ! 

CONNECT CP_M0UE< 1> : < 1> RESET_CP J10UE ; 

SEND UCO.O.Ol TO <2>RESET_CP_M0UE 

CONNECT RESET_CP M0UE< 1> : <2>CP MOUE; 

O 

MOUE _SYNC : -F : SYNC (21; (Ensures that currently 

picked vertex identifier 
is reported to the host 
computer immediately 
before the vertex 
position delta is 
reported . 1 

M0UE_TH I S_0NE : -F : CONSTANT ; 

(Sends the current 
picked vertex identifier 
to M0UE_SYNC each time 
the CP_M0UE accumulator 
fires . 1 

CONNECT CP_M0UE< 1 > : < 1> M0UE_TH I S_0NE ; 

(PICKED2<1> will be connected to 

<2>M0UE_THIS__0NE below. 1 
CONNECT MOUE THIS_0NE< 1> : < 1>M0UE SYNC; 

CONNECT CP_M0UE< 1 > : < 2> M0UE_SYNC ; 



O 



CONNECT M0UE_SYNC< 1> : < 1> TO HOST; 
CONNECT M0UE__SYNC< 2> : < 1> T0__H0ST ; 
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«***• Before reusing the accumulator, CP_MOUE, far 
***** another vertex, report any remnant in the 
***** accumulator to the host computer. 

1 



CProvides a zero to be 
sent to <3>CP_nDUE, 
allowing accumulator to 
fire without waiting to 
reach normal firing 
increment . 1 

CTrigger the accumulator 
firing increment to be 
reset to zero whenever a 
new vertex is picked for 
movement . } 

CSends a zero vector to 
<1>CPJ10UE, causing the 
accumulator to fire.} 
CSend the zero vector 
when triggered by output 
from ZERO.} 

(Ensures <3>CP_f10UE is 
set to zero before 
<1>CP.J10UE is provided 
with signal to fire 
accumulator . > 

CONNECT ZER0<1> : < 1> ZERO _F I RST ; 

CONNECT UECZERO< 1> : <2>ZEIR0_FIRST ; 

CONNECT ZERO_FIRST< 1> :<3>CP_N0UE; 

CONNECT ZERO_F I RST< 2> : < 1 > CP J10UE ; 

STEP__S I ZE : -F : CONSTANT; CAfter CP_NOUE fires, 

resets CP_f10UE firing 
increment . } 

CONNECT CP__f10UE< 1> : < 1> STEP SIZE; 

SEND 0.005 TO <2>STEP_SIZE ; 

CONNECT STEP_SIZE< 1> :<3>CPJ10UE; 

SET_PICKED2: =F: CONSTANT; CSends new value of 

PICKED to PICKED2, after 
any remnant in the 
CP_f10UE accumulator has 
been dumped to the host 
for addition to the 
position of the 
previously picked 
vertex . } 



ZERO: -F: CONSTANT; 

CONNECT P I CKED< 1> : < 1 > ZERO ; 

SEND 0 TO < 2> ZERO ; 

UECZERO : -F : CONSTANT ; 

CONNECT ZER0<1> :<1> UECZERO; 

SEND U C 0 , 0 , 0 3 TO <2>UECZER0; 
ZER0_F I RST : “F : SYNC C 2 } ; 



CONNECT PICKED< 1> : < 2> SET_P I CKED2 ; 
CONNECT STEP_SIZE< 1> : < 1>SET_ PICKED2 ; 
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PICKED2: =F : NOP; CContains identifier of 

currently selected 
vertex, but is not 
updated until CP_J10UE 
accumulator remnant has 
been reported in 
accordance with previous 
value of PICKED2 . ! 

CONNECT SET_P I CKED2< 1> : <1>PICKED2; 

CONNECT P I CKED2< 1 > : < 1 > SO I TCH„P I CK ; 

CSee SWITCH_PICK near 
end of Uertex Picking 
Network . ! 

CONNECT PICKED2<1> : < 2> NOUE_TH I S_ONE ; 

CSee MOUE._SYNC above . 3 

SEND F I X C 1 3 TO <1>PICKED2; 

C 

*-*•*•****•*•*****■*-**•*•*•*•*•**•*•* Clipping ***** ■* -m * + ■*»*■*** -m 



CL I P _TRUE : -F: INPUTS_CH00SE(3! ; 



CSend a FALSE if clip- 
ping is off, TRUE if 
clipping is on . 3 
< 1> CL I P TRUE: 



CL I P_D I AL : =F : CROUTE C 2 ! ; 



SEND FALSE TO 

SEND TRUE TO <2>CLIP_TRUE 

CONNECT FETCHCL I PNODE< 1> : < 3> CL I P_TRUE ; 

CONNECT CL I P_TRUE< 1> : < 1> CL I PON ; C Enable clipping.! 
CONNECT CL I P_TRUE< 1> : < 1> CL I P_INTENSI TY ; 

CDisable depth cueing! 
CRoutes signal from dial 
8 to control clipping 
window width if clipping 
enabled . 1 

CONNECT FETCHCL I PMODE< 1> :<1>CLIP DIAL; 

CONNECT DIALS<8> : <2>CLIP_DI AL; 

CONNECT CLI P_D IAL<1>:<1>DI SCARD ; 

CDiscard signal from 
dial 8 if clipping not 
enabled . 1 

CL I P_TOTAL : -F : ACCUMULATE ; CAccumulate dial signal! 

CONNECT CL I P_D I AL< 2> : < 1 > CL I P_TOTAL ; 

CRoutes signal from dial 
to accumulator only if 
clipping is on CCLIPflODE 
= 2 ! 

CReset value.! 

CScale dial signal and 
convert to real number.! 



SEND 0.5 
SEND 0.1 



TO <2> CLIP, TOTAL; 
TO < 4> CL I P TOTAL; 
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SEND 1 TO < 5> CLI P_TOTAL ; {Maximum clipping width} 

SEND 0.002 TO <6>CLIP_T0TAL ; Cflinimum clipping width} 
CL I P_SYNC : =F : SYNC C 2 } ; (Ensures clipping window 

limit is reset before 
the new window matrix is 



CONNECT CLIP_T0TAL<1> : <1>CLIP 
CONNECT CL I P_TQTAL< 1> : < 2> CL I P. 
Ill I NDOW J1ATR I X : -F : W I NDOlil ; 



SEND -1 
SEND 1 
SEND -1 
SEND 1 
SEND 0 



TO 

TO 

TO 

TO 

TO 



< 2> UJ I ND0U_ 
<3> WINDOW." 

< 4> W I NDOW 

< 5> W I NDOW. 
<S> WINDOW 



NATRIX 
NATRI X 
M1ATRIX 
M1ATRIX 
NATRIX 



sent to CLIPPER . 1 
_SYNC ; 

.SYNC ; 

^Provides a new matrix 
defining the window with 
the new clipping width. 1 
CWindow X minimum.} 
CWindow X maximum.} 
CWindow Y minimum.} 
{Window Y maximum . } 
{Front 2-clipping plane} 



CONNECT CL I P_SYNC< 1> : < 7> W I ND0W_T1ATR I X ; 

{Back 2-clipping plane 
defined by dial signal} 
CONNECT CL I P_SYNC< 2> : < 1> W I ND0W_NATR I X ; 

CTrigger WINDOW_nATRIX 
to send new matrix.} 



CONNECT W INDOW_MATR I X< 1> : < 1> CL I PPER ; 
C 



**************************************** ‘ft-*'*'*--*-'*-*-*--*'*'*--*-**-*-'*'*-** 
*********************************************************** 



************** Function Keys Operational ************* 

*********************************************************** 



*********************************************************** 



KEY_OFFSET : =F : ADDC ; CAdd offset to function 

key signals to put them 
into range which host 
computer program 
recognizes as indicating 
a pressed function key} 
SEND FIXC 300} TO < 2> KEY_OFFSET ; 

CONNECT FKEYS< 1> : < 1> KEY_OFFSET ; 

CONNECT KEY_OFFSET< 1> :<1>T0_H0ST; 

C 



laa 



■******-»***»****j_abel the Function Keys and Dials -*■******■*■*■* 
1 

KEY1 : =F: INPUTS CHOOSEC53 ; CStores labels for 

FKEY 11 

SEND ’ UERTICES ’ TO <1>KEY1; 

SEND ’UERTICES’ TO <2>KEY1; 

SEND ’UERTICES’ TO <3>KEY1; 

SEND ’ PRT UERT ’ TO <4>KEY1; 

CONNECT FETCHNODE< 1> : <S>KEY1 ; CNotifies KEY1 of which 

label to use, depending 
on mode . 1 

CONNECT KEY1< 1> : < 1>FLABEL1 ; CSends label to 

keyboard . 1 

KEYS : “F : I NPUTS JZHOOSE C 5 ) ; 

SEND ’ B-SPL I NE ’ TO <1>KEY2; 

SEND ’B-SPLINE’ TO <2>KEY2; 

SEND ’B-SPLINE’ TO <3>KEY2; 

SEND ’ ’ TO <4>KEY2; 

CONNECT FETCHNODE< 1> : <5>KEY2 ; 

CONNECT KEY2< 1> : < 1> FL ABELS ; 

KEY3: -F: I NPUTS_CHOOSE C S ) ; 

SEND ’ KNOWN PT ’ TO <1>KEY3; 

SEND ’ KNOWN PT ’ TO <2>KEY3; 

SEND ’KNOWN PT ’ TO <3>KEY3; 

SEND ’ ’ TO < 4> KEY3 ; 

CONNECT FETCHNODE< 1> : <5>KEY3; 

CONNECT KEY3< 1> : < 1>FLABEL3; 

KEY4:-F: INPUTS_CH00SEC5 ) ; 

SEND ’INTERACT’ TO <1>KEY4; 

SEND ’ ’ TO < 2> KEY4 ; 

SEND ’INTERACT’ TO <3>KEY4; 

SEND ’ PLT SCRN ’ TO <4>KEY l i; 

CONNECT FETCHnODE< 1> : <5>KEY4 ; 

CONNECT KEY L f< 1> : < 1> FLABEL4 ; 

O 

KEYS : =F : I NPUTS,_CHOOSE C 5 ) ; 

SEND ’CLIP’ TO < 1>KEY5 ; 

SEND ’CLIP’ TO < 2> KEYS ; 

SEND ’CLIP’ TO < 3> KEYS ; 

SEND ’ ’ TO <4>KEYS; 

CONNECT FETCHNODE< 1> : < S> KEYS ; 

CONNECT KEY5<1>:< 1> FLAB ELS; 

O 

KEYS : -F : I NPUTS _CHOOSE C53 ; 

SEND ’REFLECT’ TO <1>KEYE; 

SEND ’REFLECT’ TO <2>KEYS; 

SEND ’ ’ TO < 3> KEYS ; 

SEND ’ ’ TO <4> KEYS ; 

CONNECT FETCHnODE<l> :<S>KEYS; 

CONNECT KEY5< 1> : < 1> FL ABELS ; 
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O 

KEY7 : -F : INPUTS _CH00SEC5) ; 

SEND ’XY’ TO <1>KEY7; 

SEND ’XY’ TO < 2> KEY7 ; 

KEY7_nODE3: -F : INPUTS _CH00SEC3) ; 

SEND ’PICK ROW’ TO < l>KEY7_nODE3 ; 

SEND ’ ’ TO <2>KEY7_nODE3; 

CONNECT FETCHPICKnODE< 1> : < 3> KEY7 J10DE3 ; 

KEY7_SYNC: “F : SYNCC2) ; (Ensures neu PICKflODE 

value is in place on 
KEY7J10DE3 before neu 
NODE value is sent to 
KEY7, triggering that 
function . } 

CONNECT KEY7_nODE3< 1> : < 1>KEY7_SYNC ; 
CONNECT FETCHnODE< 1 > : < 2> KEY7_SYNC ; 
CONNECT KEY7_SYNC< 1> : <3>KEY7 ; 

SEND ’ ’ TO < 4> KEY7 ; 

CONNECT KEY7_SYNC<2> : <5>KEY7; 

CONNECT KEY7< 1> : < 1>FLABEL7 ; 

O 

KEYS: =F: INPUTS_CHOOSEC53 ; 

SEND ’ XZ ’ TO < 1> KEY8 ; 

SEND ’ X2 ’ TO < 2> KEYS ; 

KEY8_nODE3: -F: INPUTS_CH00SEC3) ; 

SEND ’ ’ TO <1>KEY8 J10DE3; 

SEND ’PICK PT’ TO < 2>KEY8_J10DE3 ; 

CONNECT FETCHP ICKnODE< 1> : <3>KEY8 J10DE3; 
KEY8_SYNC : = F : SYNC C 2 ) ; 

CONNECT KEYB J10DE3<1> : < 1>KEY8._SYNC ; 
CONNECT FETCHf10DE< 1> : < 2> KEYB _SYNC ; 
CONNECT KEY8_SYNC<1> : <3>KEY8; 

SEND ’ SAU UERT ’ TO <4>KEY8; 

CONNECT KEY8_SYNC< 2> : < 5> KEY8 ; 

CONNECT KEY8< 1> : < 1> FLABEL8 ; 

O 

KEYS: -F: INPUTS_CH00SEC5D ; 

SEND ’ Y2 ’ TO <1>KEY9; 

SEND ’ Y2 ’ TO <2>KEY9; 

SEND ’ ’ TO < 3> KEY9 ; 

SEND ’LD UERT’ TO <4>KEYS; 

CONNECT FETCHnODE< 1> : <5>KEYS ; 

CONNECT KEYS<1> :<1>FL8BEL3; 

O 

KEY 10 : =F : INPUTS..CH00SEC5D ; 

SEND ’I/O’ TO <1>KEY10; 

SEND ’I/O’ TO < 2> KEY10 ; 

SEND ’I/O’ TO < 3> KEY10 ; 

SEND ’ ’ TO < 4> KEY10 ; 

CONNECT FETCHNODE< 1> : < S> KEY10 ; 

CONNECT KEY10< 1> : < 1>FLABEL10 ; 
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KEY 11 : -F: INPUTS_CHOOSEC53 ; 

SEND ’RESET’ TO <1>KEY11; 

SEND ’RESET’ TO <2>KEY11; 

SEND ’RESET’ TO <3>KEY11; 

SEND ’ ’ TO < 4> KEY 1 1 ; 

CONNECT FETCHMODE< 1> : <5>KEY11 ; 
CONNECT KEY11< 1> : < 1> FLAB EL 11 ; 
O 

KEY12 : "F : INPUTS_CH00SEC S3 ; 

SEND ’QUIT’ TO <1>KEY12; 

SEND ’QUIT’ TO <2>KEY12; 

SEND ’QUIT’ TO <3>KEY12; 

SEND ’EXIT I/O’ TO <4>KEY12; 
CONNECT FETCHnODE< 1> : < 5> KEY12 ; 
CONNECT KEY12< 1> : < 1> FLABEL12 ; 
O 

DIAL1 : -F: I NPUTS_CH00SE C 5 3 ; 

SEND ’GLOBAL X’ TO <1>DIAL1; 

SEND ’GLOBAL X’ TO <2>DIAL1; 

SEND ’UERTEX X’ TO <3>DIAL1; 

SEND ’ ’ TO < 4> D I AL1 ; 

CONNECT FETCHMODE< 1> : <5>DIAL1 ; 
CONNECT DI AL1< 1> : < 1>DLABEL1 ; 

O 

DI AL2: -F: I NPUTS_CH00SE C 5 3 ; 

SEND ’GLOBAL Y' TO <1>DIAL2; 

SEND ’GLOBAL Y’ TO <2>DIAL2; 

SEND ’UERTEX Y’ TO <3>DIAL2; 

SEND ’ ’ TO < 4> D I AL2 ; 

CONNECT FETCHMODE< 1> : < 5> D I AL2 ; 
CONNECT DI AL2< 1> : < 1>DLABEL2 ; 

O 

DIAL3: -F: INPUTS _CH00SEC S3 ; 

SEND ’GLOBAL 2’ TO <1>DIAL3; 

SEND ’GLOBAL 2’ TO <2>DIAL3; 

SEND ’UERTEX 2’ TO <3>DIAL3; 

SEND ’ ’ TO < 4> DI AL3 ; 

CONNECT FETCHMODE< 1> : <5>DIAL3; 
CONNECT D I AL3< 1 > : < 1> DLABEL3 ; 

O 



BLINKER: =F: INPUTS CHOOSEC53 



SEND 

SEND 

SEND 

SEND 



FALSE 

FALSE 

TRUE 

FALSE 



TO 

TO 

TO 

TO 



CCause translation dials 
to blink when in inter- 
active mode.! 

< 1> BLINKER : 



<2> BLINKER 
<3> BLINKER 
<4> BLINKER 
CONNECT FETCHNODE< 1> : <5> BLINKER 
CONNECT BLINKER< 1> : < 2> DLABEL1 
CONNECT BLINKER< 1> : <2>DLABEL2 
CONNECT BLINKER< 1> : <2>DLABEL3 
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O 

DIAL4: -F: I NPUTS _CHOOSE C 5 3 ; 

SEND ’ SCALING ’ TO <1>DIAL4; 

SEND ’ SCALING ’ TO <2>DIAL4; 

SEND ’SCALING’ TO <3>DIAL4; 

SEND ’ ’ TO < 4> D I AL4 ; 

CONNECT FETCHNODE< 1>:<5>DIAL4; 

CONNECT DIAL4<1>:< 1>DLABEL4 ; 

O 

DIALS: =F: INPUTS_CHOOSEC53 ; 

SEND ’ROTATE X’ TO <1>DIAL5; 

SEND 'ROTATE X’ TO <2>DIAL5; 

SEND ’ROTATE X’ TO <3>DIAL5; 

SEND ’ ’ TO <4>DIAL5; 

CONNECT FETCHNODE< 1> : <5> DIALS; 

CONNECT DI AL5< 1> : < 1>DLABEL5 ; 

O 

DIALS: -F: INPUTS_CHOOSEC53 ; 

SEND ’ROTATE Y’ TO <1> DIALS; 

SEND 'ROTATE Y’ TO <2>DIALS; 

SEND ’ROTATE Y’ TO <3>DIALS; 

SEND ’ ’ TO < 4> D I ALS ; 

CONNECT FETCH NODE < 1>:<5>DIALS; 

CONNECT DI ALS< 1> : < 1>DL ABELS; 

O 

D I AL7 : =F : INPUTS .CHOOSEC 53 ; 

SEND ’ROTATE 2’ TO <1>DIAL7; 

SEND ’ROTATE 2’ TO <2>DIAL7; 

SEND ’ROTATE 2’ TO <3>DIAL7; 

SEND ’ ’ TO < 4> D I AL7 ; 

CONNECT FETCHN0DE<1> :<5>DIAL7; 

CONNECT DI AL7< 1> : < 1> DLABEL7 ; 

O 

DIALS: =F: INPUTS_CHOOSEC33 ; 

SEND ’ ’ TO < 1> D I ALB ; 

SEND ’CLIPPING’ TO <2>DIAL8; 

CONNECT FETCHCLIPNODE< 1> : <3>DIALB; 

CONNECT DI AL8< 1> : < 1>DL ABELS ; 

O 

SEND 1 TO < l>NEW_t10DE; CTrigger initial sending 

□f default operating 
mode indicators. 

Whenever an operating 
mode indicator 
CvariableD changes, the 
host computer will send 
a similar triggering 
message . 1 



APPENDIX D 



LINK CONN AND FILE 



The following is the 
SPLINE. FTN with the various 



command file used 
host-resident PS 300 



to link 
library 



subroutines . 



SPL INE/FP/CP/ ID, “spline, bd 
Cl , 1DMELIB/LB 
C5, 200DENSNEU/LB 

/ 

LIBR-F4PRES: RO 
ACTFIL-4 
UNITS-B 
ASG-TT4B: 7 
ASG-SYO: 6 
PR I =50 
// 
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